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1 . o Sensitivity  o?  Liquid  Monoprop* Hants  to  Explosion  by  Rapid  compress! on 


1 . l origin  of  the  Problem  in  Practical  Liquid  Propellant  Gun  (LPG) Systems 

Monopropellant  liquid-propellant  gun  systems  he/e  baan  davalopad 
generally  along  two  lines,  the  pre-loaded  type  (ao-called  bulk-loaded 
LPG)  and  tha  direct-in jaction  tyre  (so-called  regenerative  type).  In 
each  of  thaaa  types,  the  Liquid  charge  ia  brought  rather  rapi.Lly  up 
to  gun  operating  pressure,  ca.  30  to  60  kpsi,  in  a rise  tima  on  tha 
order  of  0.1  to  1 maec.  Although  hundreds  of  successful  firings  have 
been  made  with  both  types  of  gun,  w.th  several  different  types  of  liquid 
monopropellant,  in  sixes  up  to  4-inch  barrels,  a few  destructive  explosions 
have  taken  place,  all  of  them  on  start-up  during  the  rapid  pressure  rise 
associated  with  ignition  of  the  gun.  This  almost-euccessful  test-record 
suggests  that  the  particular  monopropellants  selected  are  not  themselves 
at  fault,  or  explosions  would  have  occurred  with  far  greater  frequency. 

The  most  plausible  explanation  that  has  emerged  from  all  of  tha  evidence 
available  ia  that  air  bubbles  due  to  excessive  ullage  in  the  firing 
clumber,  brought  in  during  the  pre-firing  fill  process,  become  sudden  hot 
spots  of  unusual  severity  during  the  start-up  as  a result  of  rapid 
compression.  In  at  least  one  case,  the  evidence  suggested  that  the  bubbles 
might  not  have  been  air-filled  but  simply  vapor-filled,  formed  by 
cavitation  during  the  filling  process.  It  is  this  recent  LPG  development 
record  that  forms  the  background  and  the  rationale  for  tha  present 
investigation. 

On  the  experimental  side,  Princeton  Combustion  Research  Laboratories, 
.Tnc.  has  produced  on  a systematic  basis  rapid  pressurizations  of  liquid 
propellant  (LP)  columns  with  various  rates  of  pressure  rise,  with  various 
amounts  of  ullage, both  finely  distributed  and  undistributed,  and  with 
various  L?  monopropeiiantn,  and  has  found  operating  domains  of  safe 
start-up  for  NOS -type  and  Otto  Fuel  II5*®  propellants.  Generally  speaking,  safe 
start-up  without  any  sign  of  explosion  or  runaway  exothermic  reaction 
is  obtainable  by  (a)  elimination  of  all  aii  in  the  system,  a probably 
impractical  requirement,  (b)  avoidance  of  sharp  entry  ports  or  other 
channel  locations  where  cavitation  may  occur  in  the  gun  filling  system, 
probably  also  an  impractical  requirement,  (c)  keeping  th?  rate  of  pressure 
rise  in  start-up  to  less  than  ca.  20  kpai/rasec,  a requirement  that  can 
be  met  by  rational  design  of  the  ignition  system,  and  id)  pre-pressurizing 
the  charge  of  LP  before  the  onset,  of  the  rapid  pressure  rise  of  the 
start-up,  a fairly  practical  remedy  that  can  be  incorporated  in  an  I.PG 
control  system.  These  four  orescriotions  define  the  domain  of  explosion- 
free  start-up,  on  the  basis  of  the  LP  compression  experiments  performed 
by  PCRL. 

It  should  be  noted,  as  a raattsr  of  practical  interest,  tnat  in  defining 
this  domain  of  safe  start-up,  it  was  necessary  tc  provoke  explosions  in 
the  test  apparatus.  The  domain  of  safe  start-up  is  defined  by  the  "dividing 
surface"  between  the  conditions  that  always  produce  explosion  and  those 
that  never  produce  explosion.  The  apparatus  was  fitted  in  each  case  with 
a pressure-relief  shear  disc  assembly,  vfhich  always  functioned  properly. 

The  most  significant  point  of  interest  is  that  there  was  never  any  evidence 
of  a destructive  detonation,  and  that  the  apparatus  was  always  re-usable 
after  replacement  of  the  pressure-relief  shear  disc  plug.  This  experience 
supports  the  conclusions  drawn  from  the  evidence  after  the  few  destructive 
explosions  that  have  occurred  in  lecent  practical  LP  gun  development  programs , 
that  no  detonations  have  occurred. 


With  this  background, attention  in  tha  present  rasa arch  has  been  focused 
on  tha  mannar  of  initiation  of  run  iway  exothermic  raa  tion  in  an  LP  fiald 
by  a suddenly  compressed  hot  gas  bubble  imbedded  in  the  LP . In  passing, 
we  note  that  the  problem  of  ullage  in  an  explosive  medium  is  not  limited 
to  the  field  of  LP  monopropeliant  gun  systems.  It  can  arias  in  LP 
monopropall ant  (e.g.,  hydrazine)  power  generating  systems  designed  for 
aircraft  and  spacecraft  use,  and  it  can  arise  in  the  problem  of  premature 
warhead  explosions  due  to  setback  acceleration  of  plastic-type  warhead 
explosives  in  high  velocity  projectiles.  Defining  theoretically  and 
experimentally  the  safe  domain  of  operation  for  explosive  media  containing 
bubbles  is  a task  of  genaral  significance. 

1.2  sensitization  of  LP  by  'Jllaqs  or  Dispersed  Bubbles 

When  a liquid  monopropeliant  is  rapidly  compressed  to  about  50  kpsi, 
an  satimate  of  tha  tamperaturs  risa  shows  that,  if  the  liquid  is  free  of 
bubbles  and  if  the  process  is  taken  to  be  iaentropic,  it  can  increase  by 
about  10  to  20  deg  C,  with  some  uncertainty  in  the  calculation  due  to  lack 
of  preciza  data  on  compressibility  ar.d  thermal  expansion  coefficient, 
among  othar  properties.  A bubbls-free  monopropeliant  is  therefore  unlikely 
to  explode  upon  sudden  pressurization . 

However,  with  one  or  more  small  bubbles  in  the  liquid,  a "hot  spot" 
can  be  generated  by  the  adiabatic  compression  of  the  bubble,  either  air 
or  vapor,  and  if  tha  buhble  is  initially  large  enough,  the  sharp  local  rise 
in  temperature  can  cause  the  bubble  to  retain  its  heat  long  enough  to  initiate 
exothermic  chemical  reaction,  i.e.,  combustion.  The  resulting  rapid  gas 
generation  can  lead  to  rupture  of  the  chamber  or  to  a gun  breech  blow,  \f  it  is 
fast  enough  or  if  the  voIum  expansion  does  not  accommodate  the  generated  gas. 
There  is  soma  evidence  from  firing  test  pressure  records  of  re gone rati ve- 
type  direct-injection  LP  gun s45 ar.d  from  explosion  teals  conducted  at  Princeton 
Coafcuation  Research  Laboratories,  Inc.  with  nearly  bubble-free  LP  fillinqs 
that,  with  only  a few  bubblau  sparsely  distributed  in  the  liquid  field, 
the  rate  of  gas  generation  is  slow  enough  to  be  accommodated  by  an  LPG 
system  as  the  volume  expands  at  the  start  of  a firing.  In  the  present  tests 
the  severity  of  the  explosions  in  a confined  liquid  propellant  column 
was  proportional  to  tha  space  density  of  bubbles  — few  bubbles  led  to 
only  mild  explosions.  On  the  other  hand,  whan  in  the  PCRL  compression 
twits  the  number  of  injected  air  bubbles  was  deliberately  l*.rge,  the 
resulting  explosions  were  much  raore  severe.  A*  pointed  out  above,  these 
were  never  so  eevere  as  to  destroy  the  apparatus  (except  for  the  pressure  - 
reltel  shear  diac  device)  . The  conclusion  is  that  it  takas  only  a single 
buhbl  , .5  v.  l,*cge  enough,  to  start  a run -away  reaction,  but  it  requires 
a large  number  of  sucn  oubblei  acting  simultaneously  sa  hot  spots  to 
cause  the  hardware  damagxt'  associated  with  such  explosion.. 


2.0  flow  Visual izat ion  Studies  of  Dynamic  Loading  Characteristics  off  NOS-365 
Liquid  Mono pro pa  11 ant 

In  order  to  conduct  compression-ignition  ssnaitivity  studies  of 
a pre-compressed , multiple  bubble  liquid  monopropcllant  stadium  representing 
condition!  found  in  qun  simulation  studies  performed  as  part  of  a hazards 
analysis  by  Pulsopower  Systems,  Inc-,  Princeton  Combustion  Research 
Laboratories , Inc.,  efforts  have  been  directed  toward  the  design, 
fabrication,  and  development  of  an  apparatus  which  incorporates  features 
to  control  the  rate  of  dynastic  loading  of  the  liquid  monopropellant, 
the  volume  percent  of  injected  air  ullage,  and  the  .resulting  break-up 
phenomena  of  the  ullage  Into  some  mean  bubble  distribution  in  the  liquid 
propellant  charge.  The  objective  of  the  Flow  Visualization  phase  of  this 
study  wss  to  dsfine  thr>  physical  stats  of  tha  liquid  propallant  charge 
after  dynastic  loading,  at  the  system  "f lrs"condition . 

2.1  Description  of  Test  Equipment 

Figur#  1 shows  an  assembly  drawing  of  ths  Transparent  Visualization 
Chasiber  and  the  Pneumatic  Load  L.P.  Cylinder  utilized  in  tha  Flow 
Visualization  studies.  Ths  schematic  drawing  presented  in  Figure  2 
identifies  functional  componants  such  as  valves,  gas  linas,  liquid  fill 
and  discharge  lines,  and  electrical  connections.  The  Transparent 
visualization  Chamber  was  machined  from  a cast  acrylic  block  with  s 
0. 4 375-inch  diameter  bora.  The  chamber  is  provided  with  an  End  Plug  which 
carries  ths  contact  wire  for  actuation  of  the  electronic  strobe  light. 

In  the  actual  Liquid  Propellant  Compression  Chamber,  the  End  Plug 
carries  the  contact  wire  completing  the  circuit  for  ignition  of  an 
M52  electric  primer  which  ignites  the  smokeless  powder  starter  charge 
which , in  turn,  produces  the  desired  pressurization  rata  to  which  the 
liquid  propellant  column  in  the  Compression  Chamber  is  to  be  subjected. 

The  Transparent  Visualization  Chamber  also  contains  s Bleed  Valve  for 
flushing  residual  ullage  from  the  system  prior  to  dvnamic  injection  of  the 
liquid  propellant  charge.  A Kr rstal  Type  601A  pressure  transducer  is 
fitted  into  the  chamber  to  monitor  the  pressure-time  history  of  the 
dynamic  loading  process  end  is  located  0.70-inch  downstream  of  tha  centor- 
lina  of  the  injection  port.  The  injection  port  is  designed  so  as  to  accept 
various  flow  guide r.  to  alter  the  cavitation  dynamics  of  the  injected  liquid. 
The  clamber  bore  is  fitted  with  a Separator  Piston  which  seals  ths  bore  and 
simulates  ths  function  of  the  Separator  Piston  in  ths  Liquid  Propellant 
Compression  Chiimber,  i. . e . , to  isolate  the  liquid  propellant  charqe  in  the 
Compression  Chamber  bora  from  the  hot  combustion  gasos  in  the  starter 
Charge  Chaabar  (see  Assembly  Draw^ig  Figure  3).  A Projectile  Piston 
is  utilised  to  facilitate  chamber  filling  under  dynamic  loading  and 
provides  the  trigger  to  actuate  the  electronic  strobe  light  non  chamber 
filling  completion.  The  maximum  stroke  of  the  Projectile  i-  .con  is  2.7- 
inch,  providing  a maximum  volumetric  loading  of  liquid  propellant  of 
0.41  cu.in.  (6,65  cc) . 

The  Pneumatic  Load  L.P.  Cylinder  is  machined  from  416  S.S.  The  0.5- 
inch  diameter  bore  contains  a Pneumatic  Piston  which  transfers  the  liquid 
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proponent  charge  through  a Poppet  Valve  into  the  injection  port 
of  tii*  Transparent  Visualization  Chamber.  The  Pneumatic  Load  L P. 

Cylinder  acts  as  a liquid  syringe  to  control  the  dynamic  loading  process. 

A regulated  supply  cf  Nj  gas  fed  through  a Solenoid  Valve  provides  the 
driving  pressure  for  acceleration  of  the  Pneumatic  Piston.  The 
injection  time  is  controlled  by  this  driving  pressure,  for  a fired 
volumetric  loading  of  liquid  propellant.  The  Poppet  Valve  control*  the 
steady  state  (equilibrium)  pressure  in  the  Transparent  Visualisation 
Chamber  by  adjustment  of  the  Poppet  Valve  Spring  compression.  The 
Pneumatic  Load  L.P.  Cylinder  contains  a Siliccn  Rubber  Septum  through 
which  a Gas  Ullage  Syringe,  i.e. , hypodermic  needle  passes  to  introduce 
a precise  loading  of  air  ullagu  (volume  percent,  STV,  . Liquid  propellant 
ia  introduced  into  the  Pneumatic  Loading  L.P.  Cylinder  through  a fill 
port.  A photograph  of  the  Transparent  Visualisation  Chamber  assembled 
to  the  Pneumatic  Load  L.P.  Cylinder,  with  associated  fill  lines  and 
valves  is  shown  in  Figure  4. 

2.2  System  Operation 

The  Transparent  Visualisation  Chamber  is  first  flushed  clean  of 
air.  This  is  accomplished  by  first  positioning  the  Separator  Piston 
and  the  Projectile  Piston  with  ends  on  either  side  cf  the  Visualisation 
Chamber  fill  f.  ->rt , thereby  forming  a 5/32-inch  wide  passageway  for  flow 
of  liquid  from  l.'  11  port  to  Blned  Valve . The  Poppet  Valve  is  then 
manually  held  in  t*w  open  position  as  a gravity-fed  supply  'i  liquid 
propellant  flow*  into  the  Pneumatic  Load  L.P.  Cylinder,  through  the 
Poppet  Valve,  Flow  Guide,  fill  port,  the  5/32-inch  slice  of  Visualisation 
Chamber  bore,  and  out  the  Bleed  Valve.  In  this  manner  the  entire  system 
including  the  initial  be  re  volume  between  the  two  Visualisation  Chamber 
pistons  can  be  flushed.  With  the  Poppet  Valve  released , free  to  return 
to  its  normally  closed  position,  a prescribed  volumetric  .loading  of  gas 
ullage  is  introduced  into  the  Pneumatic  Load  L.P.  Cylinder  with  the 
gas  syringe.  The  displaced  liquid  propellant  volume  is  permitted  to  flow 
out  through  the  Isolation  Valve.  All  valves  are  then  closed. 

The  dynamic  loading  sequence  is  now  initiated  by  activating  the 
circuit  that  opens  the  Solenoid  Valve  which  ralaases  the  driving  N2  pressure 
into  the  Pneumatic  Load  L.P.  Cylinder  to  accelerate  the  Pneumatic  piston. 

The  notion  of  the  Pneumatic  Piston  forces  the  liquid  propellant  and  its 
gas  ullage  past  the  Poppet  Valve  once  the  force  balance  establishes  Poppet 
Valve  opening.  The  liquid  propellant  and  its  gas  ullage  flow  through  the 
Flow  Guide,  possibly  resulting  in  cavitation  dapending  on  the  nature  of 
the  orifice,  into  the  bore  of  the  Transparent  Viaualir*ti^n  Chamber.  The 
Projectile  Piston  is  then  driven  to  the  right  allowing  the  charge  to  fill 
the  bore  until  maximum  stroke  is  achieved.  The  contact  of  the  Projectile 
Platon  with  the  contact  wire  in  the  End  Plug  completes  a circuit  that  fires 
the  electronic  strobe  light,  since  the  camera  shutter  is  open  and  the  entire 
device  ia  placed  in  a darkened  enclosure,  a photograph  of  the  resulting 
bubble  distribution  in  the  liquid  propellant  charge  is  obtained.  The  complete 
Flow  Visualization  Tester  Assembly  Procedure,  Filling  Procedure,  Checklist, 
and  Electronics  and  Photographic  Equipment  Set-Up  is  presented  in 
Appendix  C. 


FIGUPE  4.  PHOTOGRAPH  OF  TRANSPARENT  VISUALIZATION  CHAMBER  MOUNTED  TO  THE  PNEUMATIC  LOAD  L.P. 
CYLINDER  WITH  ASSOCIATED  FILL  LINES  AND  VALVES. 
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2.3  Results  of  Flow  Visualisation  Studies 
2.3.1  Iwrt  Simulant  (w«t«r)  Tests 


Function*!  costa  of  ths  system  wrt  parfomd  witi;  water  to  determine 
that  all  aspects  of  system  oparation  checkad-out.  Thes*  functional  taata 
with  watar  also  provided  5.  baaclina  for  roaipariaon  of  dynamic  loading  test* 
performed  with  NUS-365  liquid  aonopropailant. 


Flow  visualization  taata  wars  conductad  with  watar  and  air  ullaga 
voluaaa  of  0%,  0.9%  (0.060  cc) , and  3.1%  (0.206  cc) . Tha  driving  praaaura 
in  all  thaaa  taata  waa  Lixmd  at  300  paig.  Figuraa  5a  and  5b  display 
oscillograph  racorda  of  thaaa  watar  taata.  A dual  tJjnabaa*  nod*  of  data 
recording  and  display  was  utilized.  Tha  time -scale  is  as  follows. 

Tha  first  centimeter  of  tha  oscillograph  racorrd  rapraaants  40  ns.  All 
subsaquant  racording  is  displayed  at  4 ns/cn.  Tha  lower  trace  rapraaants 
praaaura  in  tha  Transparent  Visualisation  Chamber  bora.  Tha  vertical 
seating  is  200  psi/cn.  Tha  upper  trace  represents  electronic  strobe  voltage, 
showing  discharge  characteristics  one*  tha  contact  wire  in  the  End  Plug 
completes  tha  circuit.  Tha  traces  in  Figure  5b  represent  a cism  expansion 
of  the  traces  in  Figure  5a  in  tha  neighborhood  of  the  "water  hamaar"  spike . 
Table  1 contains  information  extracted  frrra  the  oscillograph  racorda  on 
Visualisation  Chamber  fill  time,  t (nsec) , peak  pressure  associated 


with  the  "water  hamknr",p  ..  (psig) , and  equilibrium  pressure  in  the 
liquid  col  sen  in  the  Visualisation  Chamber, 


P il (psig) . Depending  on 

the  ullage  loading,  tha  fill  tima  varies  from  28.0  to  32.0  msec.  The 
fill  time  is  defined  as  that  tima  from  opening  of  the  the  Solenoid  Valve 
to  time  of  attainment  of  maximum  stroke  of  tha  Projectile  Platon  (coamen- 
sur«ts  with  discharge  of  the  strobe)  . The  spike  pressure  varies  from  420 
to  460  psig  and  the  equilibrium  pressure  in  tne  Visualization  Chamber  bora, 
a function  of  the  Poppet  Valve  Spring  compression,  is  approximately  240- 
260  psig.  Tha  Popper  valve  Spring  compression  in  responaibls  for  an 
observed  equilibrium  pressure  in  tha  liquid  column  less  than  that  of  the 
applied  driving  pressure  by  approximately  50  psi . The  more  the  initial 
spring  compression , the  higher  the  restoring  fores  exerted  cm  the  Poppet 
and  the  lower  ths  resulting  equilibrium  pressure  after  injection.  This  is 
of  no  concern.  We  deliberately  use  a stiff  spring  on  tha  Poppet  Valve 
in  order  to  assure  prompt  and  complete  closing  of  th*  poppet  Valve . Thus, 
if  compression-ignition  were  to  take  place  in  the  Sensitivity  tests  to 
he  par famed  with  live  propellant,  subsequent  runaway  reaction  in  the  liquid 
would  be  confined  to  the  Liquid  Compression  Chamber  and  would  not  propa- 
gate into  the  Pneumatic  Load  L.P.  Cylinder. 


Figures  6,  7,  and  0 are  photographs  of  the  dynamically- loaded  water 
column  illustrating  the  physical  condition  of  the  charge  at  the  instant 
of  complete  fill  of  tn«»  Visualisation  Chamber  bore.  These  correspond  to 
U%,  0.9%,  and  3.1%  volumetric  air  loading,  respectively.  It  i~  interesting 
to  nota  that  zero  percent  ullage,  Figure  6,  shows  a consider*.  amount  of 
bubbles  distributed  through  the  liquid  charge  and  intense  injection  activity 
at  the  injector  port.  It  should  be  realized  that  at  the  instant  the 
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DUVING  PRESSURE > 300  palq 


RUN 

NO. 

LIQUID 

TYPE 

ULLAGE 

(voliasaS,  STP) 

tfill 

(maac) 

?spika 

(psig) 

^aquil 

(psig) 

d 

ms  an 

(inch) 

24 

WATER 

0.0 

28.0 

460 

240 

0.007 

25 

WATER 

0.9 

32.0 

420 

260 

0.009 

23 

WATER 

3.1 

30.4 

460 

260 

0.009 

11 

NOS-365 

0.0 

21.0 

540 

210 

0.008 

12 

NOS-365 

0.9 

22.0 

480 

200 

< 0.001 

14 

NOS-365 

3.1 

20.0 

540 

200 

< 0.001 

Notw:  Tha  seal*  factor  anployad  in  tha  Flow  Visualisation  photographs 
is  obtainad  by  noting  that  tha  itriatal  Typa  601a  Praasura 
Transducar  port  thraads  am  machinad  3/8-24  , i .a. , a thraad 
spacing  of  0.0417  inch. 


TABLE 


TABULATION  OF  RESULTS  OF  VISUALISATION  S1UDIES  CONDUCTED 
WITH  WATER  AND  NOS-365  LIQUID  M0N0PR0PELLANT  AT  A DRIVING 
PRESSURE  OF  300  psig. 


KRISTAL  Tyoe  601-A 
Pressure  Transducer 


Figure  6.  WATER  TEST:  0*  ULLAGE,  300  psig  DRIVER  PRESSURE 


FIGURE  7.  WATER  TEST:  0.9%  ULLAGE,  300  psiq  DRIVER  PRESSURE 
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photograph  I a taken,  the  "water  htMMr"  pressure  wav*  generat'd  tha 
extreme  right  hard  alda  of  tha  liquid  column  by  tie  slashing  of  tha 
Prv'jcctila  Piston  against  the  Pr.d  Plug  haa  not  yat  propagated  to  tha 
injection  port  location.  Therefore  the  Poppa t Valve  la  still  in  tha 
open  position,  awaiting  the  pressure  signal  to  propagate  back  through 
tha  liquid  charge.  Near  tha  injection  port,  this  cloud  of  vary  fine 
bubbles  is  believed  to  be  due  to  cavitation  in  tha  flow  guide,  producing 
bubblea  of  vatei’  vapor,  which  than  rapidly  condsnae  and  disappear  further 
downstream.  'This  cavitation  phenomenon  is  not  saen  in  the  0.9%  ullage 
case  (Figure  7) . This  could  be  due  to  the  coincidence  that  the  electronic 
strobe  was  discharged  slightly  l&tar  in  time  compared  with  the  0% 
ullage  teat,  allowing  vapor  bubblea  in  tha  cavitation  region  to  condense. 

The  »ean  bubble  disaster  for  the  aero  percent  ullage  case  taken  over 
a large  enough  sample  of  tha  bubblea  dovnatrann  of  tha  injection  port, 
is  0.00;  inch. 

Figure  7 displays  the  bubble  distribution  for  tha  case  of  0.9% 
voluMtric  loading  of  air  ullage.  Tha  Man  bubble  disaster  for  thia 
case  is  0.009  inch.  By  comparison  of  this  case  with  tits  results  of 
the  zero  percent  ullage  case,  one  could  estimate  that  the  sero  ullage 
case  may  have  about  0.1%  of  dissolved  air  and  that  thia  dissolved  air 
tends  to  coma  out  of  tha  solution  at  the  Injector  port  where  cavitation 
would  taka  place.  Thia  ia  not  surprising  since  do-aerated  water 
is  not  used  for  the  tests. 

Figure  0 displays  tna  bubble  distribution  for  the  case  of  3.1% 
volumetric  loading  of  air  ullage.  The  mean  bubble  diameter  for  this 
case  is  0.009  inch.  Intense  turbulent  mixing  can  be  seen  in  the  vicinity 
of  tha  injector  port. 

• 2 ^OS-365  Liquid  Monopropellant  Taste 

The  first  series  of  liquid  monopropellant  dynamic  loading  rests  wars 
conducted  at  a driving  pressure  off  300  paig  with  ullage  volumes  off  0%, 

0.9%,  and  3.1%.  Figures  9a  and  9b  display  oscillograph  records  of  these 
propellant  loading  tests.  Tha  records  of  Figure  9b  represent  time 
expansions  of  the  traces  in  Figure  9a  in  the  neighborhood  of  the  pressure 
spike.  Propellant  fill  timea,  extracted  from  Figure  9a,  are  considerably 
less  than  corresponding  fill  timet  for  water  loading  tests,  i.e.,  ca.  20 
msac  for  NOS-365  compared  to  30  usee  for  water.  Also,  the  spike  pressure 
is  increased  to  ca.  540  psig.  A lower  equilibrium  pressure  is  attained 
in  the  propellant  column,  associated  w?.th  a slightly  increased  Poppet 
Valve  Spring  compression. 

Figures  10,  11,  and  12  are  photographs  of  the  dynamically-loaded 
propellant  column.  Thexe  corraspond,  respectively,  to  0%,  0.9%,  and 
3.1%  volumetric  loading  of  air  ullage.  The  zero  percent  ullage  case,  Figure 
10,  shows  results  very  similar  to  that  o t water,  with  intense  cavitation 
action  near  the  injector  port  and  a mean  bubble  diameter  in  the  column 
downstream  of  the  injection  port  of  approximately  0.008  inch. 
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(b)  0.9%  Ullage 


FIGURE  9a. 
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FIGURE  9b.  EXP  XSioN  OK  OSCILLOGRAPH  RECORDS 
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FIGURE  10.  NOS -36 5 TEST:  0%  ULLAGE,  300  psic;  DRIVER  PRESSURE 


FIGURE  11.  NOS-365  TEST:  0.9%  U1J.AGE , 300  psig  DRIVER  PRESSURE 


Figure  11,  depicting  the  physical  condition  of  the  dynamically-loaded 
propellant  charge  with  0.9%  volumatric  air  loading,  displays  a "fog-like" 
appearance.  The  swan  bubble  diameter  is  less  than  0.001  inch  for  this 

caae. 


Figure  12,  the  case  of  3.1%  volumetric  air  loading  produces  a "fog" 
that  is  nearly  opaque.  No  bubbles  of  measureable  diameter  can  be  seen. 

It  is  reasonable  to  assume  that  the  mean  bubble  diameter  for  this  case 
is  considerably  lec  ' than  0.001  inch. 

We  can  suemari ae  these  results  as  follows.  Under  identical  conditions 
of  ullage  load: Jig  and  driving  pressure,  the  bubbles  generated  in  NQS-365 
liquid  manopropellant  are  considerably  smaller  than  those  generated  in 
water,  by  a factor  on  the  order  or  10.  We  suspect  that  surface  tension 
plays  a aajor  contributing  role  in  altering  the  ultimate  size.  It  may 
be  worthwhile  for  further  experimentation  to  be  performed  with  water 
and  added  surfactants  and  perhaps  with  more  viscous  liquids  in  order  to 
understand  sore  accurately  the  factors  that  determine  the  bubble  size. 

This  could  be  of  immense  practical  importance  as  a means  of  ensuring  that 
all  bubbles  created  in  a practical  LPG  fill  procedure  are  smallest  possible 
size (safest  when  considering  coeipression-ignition  sensitivity) . 

We  anticipate  that  bubble  size  can  be  controlled  by  adjusting  the 
driving  pressure  and  hence  the  time  of  injection  associated  with  the 
Pneumatic  Load  L.P.  Cylinder.  We  expect  that  bubbles  of  finer  or  coarser 
diameter  can  be  established  by  choosing  either  faster  or  slower  fill  times, 
i.e.,  higher  or  lower  driving  pressure.  Two  additional  series  of  tests 
were  conducted  with  N05-365  liquid  monopropellant  with  driving  pressures  of 
5 00  psig  and  ISC  peig.  note  that  the  150  psig  driving  pressure  is  just 
sufficient  to  drive  the  pneumatic  piston  and  open  the  Poppet  Valve, 
thereby  injecting  liquid  into  the  Transparent  Visualization  Chamber. 

Figures  13c  and  13b  display  oscillograph  records  of  propellant  dynamic 
loading  teats  with  a driving  pressure  of  5C0  psig.  Table  2 contains 
information  extractsd  from  the  oscillograph  records  on  visualization  chamner 
fill  time , spike  pressure  , and  equilibrium  pressure  in  the  liquid  column. 
Figures  14,  15,  and  16  are  photographs  of  dynamically-loaded  propellant 
column  under  a driving  pressure  of  500  psig.  These  correspond,  respectively, 
to  0%,  0.9%,  and  3.1%  voliasetrlc  loading  of  air  ullage.  Figure  14, 
the  zero  percent  ullage  case,  shows  significant  departure  from  the 
equivalent  case  at  a driving  pressure  of  300  psig.  Finer  bubbles  can  now 
be  observed,  dispersed  throughout  the  column.  The  mean  bubble  diameter  is 
reduced  to  0.004  inch,  approximately  one-half  of  that  observed  for  the 
0%  ullage,  300  peig  driving  pressure  case.  Figure  15,  the  G.9%  ullage 
caae  at  500  psig  driving  prassure , shows  remarkable  similarity  to  the 
opaque  "fog”  produced  in  the  3.1%  ullage,  300  psig  driving  pressure  case. 

The  mean  bubble  diameter  is  estimated  to  be  considerably  less  than  0.001 
inch.  Increasing  the  ullage  to  3.1%  volumetric  loading  at  a driving 
pressure  of  500  psig  results  in  a photograph  (Figure  16)  much  li :e  that  of 
the  0.9%  ullage  case.  The  fineness  of  the  ullage  break-up  into  tiny  bubbles 
cannot  be  resolved  to  any  bettar  than  0.001  inch. 


(a)  0%  Ullage 


► m« 4 msec /cm ► 

I msec  > 


(b)  0.9%  Ullage 
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(c)  3.1%  Ullage 


FIGURE  13a.  OSCILLOGRAPH  RECORDS  OF  NOS-365  TESTS:  500  psig  DRIVER  PRESSURE 
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FIGURE  13b.  EXPANSION  OF  OSCILLOGRAPH  RECORDS 
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DEIVIMG  y««8PMi  500  p»iq 


RUN 

NO. 

LIQUID 

TYPE 

ULLAGE 

(VOLUME  »,  STP) 

tflll 

(MAC) 

P*pik« 

(psig) 

P«quil 

(psig) 

d 

wmm i 

(inch) 

20 

NOS-365 

0.0 

20.8 

1100 

350 

0.004 

23 

NOS-365 

0.9 

22.0 

800 

400 

< 0.001 

27 

NOS-365 

3.1 

16.0 

900 

.100 

<0.001 

TABU  2 


TABULATION  OP  RESULTS  OP  VISUALISATION  STUDIES  CONDUCTED  WITH 
WOE-365  LIQUID  MONOPROPELLANT  AT  A DRIVING  PRESSURE  OP  500  pwig 


FIGURE  16.  NOS-365  TEST:  3.1%  ULLAGE,  SOD  psig  DRIVER  PRESSURE. 
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photographs  of  the  dynastic  loading  te<ts  at  a driving  praaaura  of 
ISO  paig  ara  shown  in  Figuraa  17,  18,  and  19  corresponding  to  0%,  0.9%, 
and  3.1%  volumetric  air  loading.  Figure  17,  tha  0%  ullage  case,  shows 
a liquid  coluan  relatively  free  of  bubbles,  indicating  that  tha  cavita- 
tion phenosMnor.  is  greatly  reduced  at  this  low  driving  pressure.  The 
largest  disaster  bubbles  can  be  iteen  in  the  center  of  the  coluan  length 
bubbles  with  disaster  as  large  a 0.023  inch  can  be  observed.  The  effect 
of  the  low  driving  pressure  cn  ullage  break-up  is  sore  pronounced  in 
Figure  18,  the  0.9%  ullage  caee . Large  disaster  bubbles  and  even  undivided 
ullage  can  be  seen  superiaposed  or.  a field  of  well  dispersed  bubbles  of 
aaan  diasater  0.009  inch.  Bubbles  with  dissMter  as  large  as  0.045  inch 
can  be  s«en  in  the  coluan.  Note  that  the  "fog-like"  quality  of  the 
liquid  coluan  observed  in  higher  driving  pressure  tests,  i.e.,  300  psig  ond 
500  psig,  with  NOS -365 , and  0.9%  ullage  is  not  observed  at  this  lower 
driving  pressure.  The  subdivision  of  the  original  ullage  loading  in  the 
flow  guide  and  injector  port  is  not  as  intense  at  low  driving  pressure. 

This  saas  observation  is  borne  out  for  the  3.1%  ullage  case,  shown  in 
Figure  19.  The  aean  bubble  diaaeter  for  this  case  is  0.006  inch,  with 
bubbles  well-distributed  through  the  coluan  length.  The  "fog-like" 
opaque  liquid  coluan  mj  observed  for  the  3.1%  ullage  case  at  driving 
pressures  of  300  psig  and  500  psig  is  not  observed  at  150  psig  driving 
pressure . 


FIGURE  17.  NOS-365  TEST:  0%  ULLAGE,  150  ;.siq  DRIVER  PRESSURE 
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.4  Summary  of  Flow  Visualization  Test  Results  on  Dynamic  Loading  Characteristic* 

of  a Liquid  Charge. 

1.  The  combination  of  Pneumatic  Load  L.P.  Cylinder  and  Liquid  Propellant 
Compression  Chamber  produces  the  desired  physical  state  of  the  dynam- 
ically-loaded liquid  propellant  charge.  A pre -compressed,  multiple 
bubble,  NOS-365  medium  has  been  produced  representing  similar  conditions 
found  in  other  gun  simulation  studies. 

2.  The  versatility  of  the  system  so  designed  has  been  established.  Para- 
meters under  precise  control  are  pre-compression  liquid  pressure  level, 
fraction  of  initial  ullage  voIism,  mean  bubble  size  in  the  dynamically- 
loaded  column,  and  nature  of  the  occluded  gas. 

3.  The  physical  state  of  NOS-365  liquid  monopropellant  with  varying 
fraction  of  air  ullage,  upon  dynamic  loading  into  the  Transparent 
Visualization  Chamber  bore,  is  markedly  different  from  that  observed 
with  inert  (water)  charges.  In  general,  an  opaque  "fog- like"  appear- 
ance of  the  NOS-365  column  results  with  dynamic  load,  for  driving- 
pressures  in  the  Pneumatic  Load  Cylindar  greater  than  300  psig. 

Subsequent  compression-ignition  analyses  must  take  this  physical 
observation  into  account. 

4.  The  formation  of  bubbles  in  a dynamically- loaded  "neat'  propellant 
charge  appears  to  result  from  cavitation  effects  in  the  flow  guide 
and  injection  port.  It  is  estimated  that  the  resulting  void  volusw 
may  be  as  much  as  0.1%  of  the  total  available  chamber  bora  volume . At 
the  "fire  condition"  of  the  caecression- ignition  sensitivity  tests  to 
be  performed  in  the  next  phase  of  the  project,  it  is  anticipated  that 
the  cavitation  vapor  bubbles  will  not  have  returned  to  solution,  since 
the  time  of  "fire"  is  precisely  the  time  that  the  photographic  film 
was  exposed. 

5.  Observed  mean  bubble  diameter  and  distribution  of  bubbles  in  the 
dynamically-loaded  liquid  propellant  column  are  quite  repeatable  in 
the  300  psig  and  500  psig  driver  pressure  tests.  At  low  driving 
pressure,  i.e.,  150  psig,  tha  nature  of  the  ullage  break-up  is 
sensitive  to  the  condition  of  the  Pneumatic  Load  Cylinder  bore  surface 
and  resulting  "0"-ring  friction.  The  subdivision  of  the  initial  ullage 
upon  dynamic  load  is  lass  reproducible  at  this  extremely  low  driving 
pressure.  It  has  bean  observed  that  pockets  of  gas  are  superimposed 

on  a distribution  of  fine  bubbles  for  low  driving  pressure  tests.  The 
characteristic  "fog-like"  appearance  of  the  resulting  liquid  propellant 
column  is  not  present  at  this  low  driving  pressure. 

6.  it  would  definitely  tfe  worthwhile  for  experiments  to  be  performed  with 
water  and  added  surfactants  and  with  other  liquid  monopropellants , i.e., 
■’LGP  1845",  to  understand  the  factors  that  determine  bubble  size  and 
distribution  in  a dynamically-loaded  co.l\mn.  It  would  also  be 
instructive  to  alter  the  shape  of  the  cavitation  insert  flow  guide  to 
purposefully  alter  the  nature  of  the  cavitation  phenomenon.  This,  of 
course,  all  has  a bearing  on  compression- ignition  sensitivity  of  the 
resulting  dynamically-loaded  liquid  propellant  column. 
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3.0  Starter  Charge  Pressurization  Rata  Tests 

In  order  to  establish  the  threshold  condition  for  runaway  reaction 
in  a liquid  monopropellant  leading  to  severe  over-pressures  and  possible 
explosion,  attention  must  first  be  directed  toward  control  of  the  pres- 
surization rate  to  which  the  dynamically-loaded  liquid  monopropellant  column 
is  subjected.  Previous  studies  performed  by  Princeton  Combustion  Research 
laboratories,  Inc.,  Ref.  51,5%  indicate  that,  for  a given  mean  plateau 
pressure,  the  pressurization  rate  is  an  important  parameter  in  determining 
sensitivity  to  compression-ignition.  In  those  studies  we  found  that  the 
pressure-time  start-up  histories  could  be  brought  under  precise  control  by 
varying  the  type  of  smokeless  powder  charge  in  the  pressurizing  chamber, 
in  those  studies  the  range  of  achievable  pressurization  rates  was 
1.25  <_  dp/dt  (kpsi/asec)  £ 150.  Thus,  the  objective  of  this  phase  of  the 

study  was  to  establish  "recipes",  i.e.,  starter  charge  loading  density  and 
types  of  smoke loss  powders,  to  produce  various  pressure-time  start-up 
histories  to  which  the  dynamically-loaded  liquid  monopropellant  charge 
would  be  subjected  for  cotnpra ssion-igniton  sensitivity  analysis. 

3.1  Description  of  Test  Equipment 

Figure  20  shows  an  assembly  drawing  of  the  starter  Charge  Tester 
utilized  in  the  Starter  Charge  Pressurization  Rate  Tests.  The  schematic 
drawing  presented  in  Figure  21  identifies  electrical  connections  for  joth 
instrumentation  and  firing  circuit,  vent  gas,  smokeless  powder,  and  wadding 
locations.  The  Starter  Charge  Chamber  is  machined  from  a 3.5  inch  diameter 
x 8 inch  long  bar  of  type  17-4 PH  stainless  steel  heat  treated  to  condition 
H 1150  prior  to  machining.  The  chamber  bore  diameter  is  1.00  inch.  One 
end  of  the  ch*toer  bore  ha*  a 1 3/8-8M-2B  internal  threaded  portion  to 
accept  the  Electrical  Primer  Contact  Assembly  (PA  BUBBLE-1300) . The 
other  end  of  the  chamber  has  a "neck"portion  1.3725  inch  diameter  and 
1.250  inches  long.  This  "neck"  is  machined  with  an  0-ring  groove  for 
sealing  purposes  when  the  Starter  Charge  Chamber  is  mated  with  the  Blind 
Chamber.  A 0.88  inch  length  of  the  chamber  body  is  machined  with  3 1/2-8N-2A 
threads.  Theme  threads  mate  with  the  internal  threads  of  the  Chamber 
Coupling  Collar  to  join  the  Starter  Charge  Chamber  and  Blind  Chamber  together. 

The  Blind  Chamber  (PA  BQBB'LE-1005-A)  was  designed  to  allow  calif  it  ion 
testing  of  the  various  starter  charge  mixes  with  an  internal  geometry 
that  duplicates  that  of  che  more  expensive  and  coaplex  Liquid  Propellant 
Oaapression  Chamber.  It,  too,  is  machined  from  type  17-4 PH  stainless  steel 
in  the  beet  treated  condition  H 1150.  The  Blind  Chamber  is  3.5  inch 
in  dinar  by  2.5  inch  long  with  a 0.95  inch  length  of  3 1/2-8M-2A 
external  threads  for  mounting  to  the  Coupling  Collar.  The  central  bore 
of  the  Blind  Ovs^ser  , into  which  the  "neck"  of  the  starter  Charge  chamber 
fits,  is  1.375  inch  in  diameter  by  1.240  inch  deep.  A small  bore  0.25 
inch  diameter  by  0.50  inch  long  simulates  the  pressure  communication  passage 
between  the  Starter  Charge  Chamber  and  the  isolated  liquid  column  in  the 
Liquid  Propellant  Compression  "Chamber . A PCB  Type  119a  piezoelectric  gage 
is  fitted  into  the  Blind  Chamber  to  monitor  the  pressure-time  history  within 
this  passage  due  to  the  starter  charge  combustion.  This  P-t  history  is 
that  which  will  be  sensed  by  the  Separator  Piston  compressing  the  liquid 
ooiimn  during  compression-ignition  sensitivity  testing. 
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The  llwt  Shield  {P/M  BUBBLE-1008)  fits  ovtr  tha  1/4  inch  dia— ter 
pressure  n— isiicstior  puuqt  in  the  Blind  chamfear  and  ie  designed  to 
prevent  tbs  direct  blast  of  burning  powter  and  wadding  from  an  taring  tha 
pass. age.  It  ia  0.87  inch  in  dia— tar  by  0.50  inch  high  and  aachinad  from 
17-4  PM  stainlass  ataal  hardanad  to  condition  H 1025  aftar  —chining.  Six 
0.12  inch  diaawtar  holaa  are  evenly  spacad  around  the  shiald  pariphary  and 
oonnact  to  v.  osntar  bora  0.375  inch  in  dia— tar  by  C.25  inch  deep.  The— 
holaa  allow  gas  presaura  to  comnunicata  to  tba  flow  passage  in  tha  Blind 
Ch— bar . Two 0.152  inch  dia— tar  holaa  located  (1625  inch  apart  paaa  longi- 
tudinally through  tha  Blast  Shiald  to  acooandata  Mounting  screws. 

Tha  Chanbar  coupling  Collar  (P/M  BOBBLE -100 3)  ia  aachinad  from  a 2-inch 
thick  plata  of  AX8X  1020  carbon  ataal.  It  ia  5 inches  0.0.  with  3 1/2-8M-28 
thread*  spanning  tha  2"  length  of  tha  oollar.  A milled  flat  with  two 
3/8-16  UMC-2B  tapped  holaa  located  1 7/8  inch  apart  provide  a —ting 
surfaoa  for  tha  collar.  Whan  tho  Starter  Charge  chanbar  and  tha  Blind 
Chanbar  are  joined  together  and  threaded  into  tha  Coupling  Collar*  tho 
entire  taster  is  supported  by  this  mounting  surface.  Tha  Chanbar  Coupling 
Oollar  provides  a si— la*  high  strength,  fastening  ache—  for  tha  chanbar 
sections. 

tha  combination  of  tha  Primer  Gas  Intagrating  CPU—)  Chamber  (P/N 
BOBBLE -1006 ) , tha  Orifice  Plug  (P/N  BUBBLE-1007) , and  the  Electric 
Primer  (01 in  8M-52A3B1)  constitute  the  igniter  assembly  for  tha  starter 
charge.  When  inserted  into  tha  Starter  Chergo  Chamber,  this  assembly 
also  provides  a presaura  seal  for  tha  chamber.  Tha  Integrating  (Pun)  Chaatoer 
ia— chi— d from  AXSX  1020  carbon  steal  and  is  1.63  inch  long.  0m  and  has 
a —chined  bore  to  accept  tha  electric  pri— r.  The  M-52A3B1  pri— r ia  a 
canter  button  contact,  conductive  mix  type  with  a cup  a— 1 brass  casing. 

The  main  body  of  tha  Integrating  (Fuse)  Chamber  is  1.000  inch  diameter 
with  an  "O'* -ring  groove  -chined  0.63  inch  from  the  front  edge.  A shoulder 
1.225  inch  and  0.375  inch  thick  acts  as  a atop  whan  tha  Integrating  (Fuse) 
Chamber  assembly  is  inserted  into  tha  Starter  Charge  Chanbar.  Tha  —in 
bora  of  the  Integrating  (Fuse)  Chanbar  ia  0.55  inch  i.d.  by  0.75  deep.  Tha 
entrance  bora  ia  —chined  with  3/4-10  UK-2B  threads  for  —unting  tha 
Orifice  Plug. 

Tha  Orifice  Plug  is  —chined  from  a Grade  8 alloy  ataal  3/4-10  DM  bolt. 
Tha  plug  is  0.55  inch  long  and  has  four  (4)  0.0935  inch  diameter  thru 
holes  located  evenly  spaced  on  a 0.3  inch  dia— ta^r  circle.  Tha—  orifice 
holes  —ter  the  hot  pri— r ga— a emerging  from  tha  Integrating  (Fuse)  Chamber 
into  the  main  Starter  charge  Chamber. 

When  threaded  into  tha  Starter  Charge  Chamber,  the  Primer  Electrical 
Contact  Assembly  (P/N  BUBBLE- 1300)  retains  tha  Intagrating  (Fuse)  Chamber 
assembly  and  provides  tha  electrical  curren'  to  tha  prixmr . The  Contact 
Assembly  is  2 1/4  inches  long  with  a 1 1/4  inch  wrench  hex  and  1 3/8-0M-2A 
threads.  The  Pri— r Electrical  Contact  Housing  (P/N  BOBBLE- 1301)  is  made 
of  AISI  1020  carbon  stael  rod  and  ia  — chined  to  accept  a BMC  electrical 
connector  at  tha  hex  end.  A 0.076  inch  dia— ter  hole  in  tha  housing  face 
alio—  the  spring  loaded  Pri— r Electrical  Contact  Pin  to  make  contact  with 
the  center  button  of  the  pri— r. 
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Captation  products  of  tha  saoke lass  powder  r tartar  charge  ara 
continuously  van  tad  froa  tha  chanbar  thru  tha  Powdsr  Chaster  Vant  Plug 

(p/m  Mmaut-1009) . 

The  vant  consists  of  an  AXSX  1020  carbon  staal  cylinder  0.410  inch 
dimeter  by  0.635  inch  long.  Tha  sealing  cona  at  ona  and  is  88  dagraas 
includad  angla.  A 0.191  inch  dl  ana  tar  by  0.50  inch  long  hard  coppar 
insert  is  press  fit  into  tha  bora  in  tha  staal  body.  Tha  coppar  carrias 
a C.0S3  inch  disaw  tar  thru  hols  which  serves  as  tha  vant  hole.  Coppar 
is  used  because  of  resistance  to  high  tasperature  erosion.  A 1/2-13  hex 
bolt  with  a 3/16  inch  diaswtar  thru  bora  retains  tha  vant  in  tha  Blind 


A photograph  of  tha  sssaeblad  starter  Charge  chanbar  and  Blind  Chanbar 
is  shown  noun  tad  in  tha  test  stand  in  Figure  22  . 

3*2  Systsn  Operation 

Tha  Integrating  (Fuss)  Chanbar  is  first  loaded  with  aaokeleaa  powder. 
In  all  cases  DuPont  I MR  4198  was  used.  Tha  iwss  loading  can  be  varied, 
depending  on  tha  nature  of  the  start-up  pressurization  history  desired 
in  the  main  starter  Charge  Chanbar.  It  has  been  our  experience  that  pres- 
surisation rata  dalaya  (a  pre-presaurisaticn  "stew")  can  be  achieved  with 
such  an  igniter  configuration  by  dacreasing  tha  charge  loading  in  tha 
Integrating  (Fuse)  Chanbar . This  will  be  discussed  further  in  the  next 
section. 

Once  the  Integrating  (Fuse)  chsnber  is  loaded,  tha  Orifice  Plug  Is 
threaded  in  place.  The  gas  was  flow  into  tha  aain  starter  Charge 
Chanter  can  be  precisely  astsred  by  appropriate  choice  of  Orifice  Plug. 

A thin  layer  of  cellophane  taps  la  used  to  cover  the  orifice  holes 
to  retain  the  powder.  The  Integrating  (Fuss)  Chanter  and  Electric 
Priawr  Contact  Assanbly  are  than  fitted  into  the  Starter  Charge  chsnber. 
Finally  tha  raamlnlng  loading  of  aankeleas  powder  la  introduced  into  tha 
bore  of  tha  chanbar  and  retained  with  wadding.  Tha  aospluta  Starter  Charge 
Tbster  Assanbly  Procedure  is  presented  in  Appendix  D. 
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3.3  Bsaulta  of  itartar  Cteg*  T»«t» 

A cautious  approach  was  taken  in  the  a tar tar  char?*  tast  aarlas  in 
acdar  to  aaaur*  proper  function  in?  of  the  starter  char?*  ch— bar  with 
aaaociatad  alactric  primer  aaaaahly,  fua*  section , and  orifioa  plug. 

Supporting  tharnochasd  cal  equilibrium  caJ  culatlona  for  DuPont  MR  4198 
Smokeless  Powder  wore  perforate  with  the  ,>BLM2"  aosputar  oodo  to  aatabliah 
peak  chanter  prassur*  as  a function  of  loading  density  in  tha  chanter. 

The  results  of  this  computation  are  shown  in  Piguro  23  , a plot  of 
predicted  equilibria  cheater  prassur*  vs.  loading  density.  Star tar 
charge  taste  ware  parforaad  far  p"  • 0.075,  0.125,  end  0.165  g/cc.  in  all 
caaae  2.60  g of  INS  4198  powdar  ware  loaded  into  tha  Integrating  (Pus*) 
Chaster.  As suiting  peak  prsssur*s,as  sanaad  hy  tha  piasoalactric  gag*  in 
ths  Blind  Chaster,  are  plottad  in  Pigur*  23  , raprassntad  by  opan  circlss. 
tha  firing  hast  resulting  in  a peak  chaster  pressure  of  25  kpsi  correspond* 
to  a total  of  10.15  g of  DuPont  INK  4198  Sate* less  Powder  loading.  By 
reducing  tha  loading  of  tha  powder  charge  in  tha  Integrating  (Puna)  chanter 
the  onset  of  rapid  chaster  preasurisation  can  toa  delayed.  This  correspond* 
to  a low  prassur*  pr# -compressive  interval  in  the  p-t  history,  the  duration 
of  which  is  a function  of  tha  powdar  loading  in  tha  Pus*  section.  Like- 
wise, a faster  burning  saekeless  powder,  i eh  u Heroo  shotgun  Powder  or 
Baroo  Pistol  Powdar  will  result  in  a wore  . apld  chanter  pressurization. 

A sunnary  plot  of  p-t  start-up  histories  for  various  starter  charges 
and  loading  densities  is  shown  in  Figure  24  The  range  of  achievable 
-tressurizatian  rates  is  seen  to  be  25  <_  dp/dt  (kpsi/nsac)  100.  Tha  fastest 
w^ssurisation  rata,  i.e.,  100  kpsi/asec,  is  obtained  with  narco  Pistol 
~ 'tear  V ths  wain  chanter  and  an  MR  4198  fuse  (Cure*  E) . Ths  slowest 
j*sv»  ..tion  rata,  i.a.,  25  kpsi/saac,  is  obtained  with  0.65  g of  MR 
4198  and  10.57  q of  MR  4198  in  tha  aain  chaster  (Curve  A).  Curve  A 
also  displays  a low  pressure  pre-coapression  interval,  lasting  sppraxiM&eiy 
2.5  boss  prior  to  rapid  pressurization . 

Starter  charge  pressurization  rats  tests  ware  limited  to  swan  peak 
chanter  pressure*  of  25  kpsi,  although  the  apparatus  has  been  designed 
to  withstand  tasting  with  a naan  peak  chamber  pressure  of  50  kpsi.  It  wee 
felt  tost  all  the  phenomena  of  interest  would  be  present  at  25  kpai  naan 
peak  chanter  prassur*.  Also,  since  large  excursions  from  the  mean  pressure 
can  be  ar^ '^ip*.  In  the  liquid  propellant  compression  chamber  if  explosion 
wart  . -.cur-  * -alt  it  necessary  to  confine  sensitivity  testing  to  25 
kpsi  mean  peak  chaster  pressure  in  order  to  assess  the  magnitude  of  any 
pressure  spikes. 
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SMOKELESS  POWDER  LOADING  DENSITY,  g/cc 

pigur*  23. COMPARISON  OF  CHAMBER  FIRINGS  WITH  DuPONT  IR1  4198  SMOKELESS 
POWDER  AND  THERMOCHEMICAL  EOUIL'BRIOM  CALCULATIONS. 
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ai«a»Ti»G  aiMiu  km* 


(FP8K  tactjom)  

C.S3g  MWt  UOt  4198  10.57g  CMPont  HI*  4198 

1.30g  DuPont  M 4198  9.92g  Went  II*  4198 

2.80g  tMPont  XMft  4198  7.S5g  IXtPont  m*  4198 

2.80g  Dtfwt  XNR  4198  7.00g  — roo  Shotgun  Vendor 

2.60g  DuPont  TMll  4198  8.50g  Horco  Pistol  Pcwdor 
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4.0  Pffwrraaalon-Iqnition  Sensitivity  studies  of  4 Dynamically  Loaded 
Liquid  Propellant  Charge  " 

4,1  Description  of  Tut  Bquipaaut 

Figure  25  shows  «n  M’/  drawing  of  tha  L.P.  Cowpraaa ion- 

ignition  Sanaitivity  IMater  uaad  in  tha  experimental  atudias.  functional 
components  such  as  valves,  gaa  linaa,  liquid  fill  and  diacharga  linaa , 
gaa  vent,  and  electrical  connection*  nra  prevented  in  the  acheauitic  drawing 
figure  26  . A couplet*  description  of  the  Starter  Charge  Cheaper  and 
aseociatad  components  ia  given  in  section  3*1  The  L.P.  Pomps, aaelon 
Chamber  (P/S  SC— LB- 1002  - C ) as  tea  with  the  Starter  Charge  Oiaahnr  via  the 
Chamber  Coupling  Collar . The  Compression  dhaabar  ia  aechlned  froa  type  17 -4m 
stainless  steel  in  the  8 1150  heat  treated  condition  and  Measures  3.5  inch  in 
iiaw*ttr  by  7 3/16  inch  long.  vein  bore  of  the  rhnabar  ia  0.43S2  inch 

fllwartar  and  carries  tha  Separator  Piston  (P/N  scbslS-1100-a)  and  the  Projectile 
Piston  (P/M  SO—  LI- 1013) . The  Separator  Piston  prevent*  hot  combustion 
geees  froa  the  Starter  Charge  Chamber  froa  igniting  the  liquid  propellant 
charge  in  the  Compression  Chaaber  bore.  The  Projectile  Platon  is  utilised 
to  facilitate  chaaber  filling  under  dynamic  loading  and  prvvidra  the  trigger 
to  act' vats  tha  ignition  of  the  Starter  Charge.  The  Sant  Plug  (P/M  BURBLE- 1009- a) 
aounts  in  the  Compression  Chaaber  and  provides  continuous  venting  of 
starter  charge  combustion  gases  through  its  0.053  inch  disaster  copper 
orifice.  Sard  drawn  copper  ia  used  for  the  vent  orifice  insert  because  of 
its  resistance  to  hot  gaa  arosion. 

Tha  art  ream  and  of  the  chaaber  is  machined  to  accept  the  aluainua  Sheer 
Disc  (P/N  SO— LB-1015 ) and  the  Disc  Retainer  (P/N  BUMS LB-1016-A) . Tha  Sheer 
Disc  acta  ar  a step  for  tha  Projectile  Platon  during  dysssic  loading  and, 
in  tha  event  of  s;  Gvarpr»jw»ur#  (exceeding  65  kpsi)  due  to  L.P.  combustion , 
the  Projectile  Piston  will  shear  thru  the  disc  and  allow  the  chaaber  to 
aapty  and  depressurise  rmpldly.  Mounted  an  the  Disc  Retainer  is  tint  Micrc- 
svitch  Aaasably  which  fires  tha  electric  prime-  when  dynamic  loading  is 
complete . The  L.P.  Coepresaion  Chaabar  also  contains  a Bleed  Valve  for 
flushing  residual  ullage  froa  the  aystaa  prior  to  dynastic  loading  of  the 
liquid  propellant  charge.  Tha  injection  port  ia  designed  so  as  to  accept 
various  Plow  Guides  to  altar  the  cavitation  dynaaics  of  the  injected  liquid. 

The  Pneumatic  Load  L.P.  Cylinder  (P/N  SO— LB-1004)  ia  machined  from 
type  416  stainless  steel  heat  treated  to  Rockwell  C30/38.  The  0.50  inch 
diameter  bore  contains  a Pneumatic  Piston  (P/N  BUB6L2-1021)  which  transfers 
the  liquid  propellant  charge  through  a Poppet  Valve  into  the  injection  port 
of  the  L.P.  Compression  Chaaber.  The  Pneuaatic  Load  L.P.  Cylinder  acts  as 
a liquid  syringe  to  control  the  dynamic  loading  process.  A regulated  supply 
of  Nj  gea  fed  through  a Solenoid  Valve  provides  tha  driving  pressure  for 
acceleration  of  the  Pneuaatic  Piston.  Tha  injection  time  ia  controlled  by 
this  driving  pressure,  for  a fixed  volumetric  loading  of  liquid  propellant. 

The  Poppet  Valve  controls  tha  steady  state  (equilibrium)  pressure  in  the 
L.P.  COepression  Chaaber  by  adjusting  tha  Poppet  Valve  Spring  coopression. 

The  Pneuaatic  Load  L.P.  Cylinder  contains  a Silicon  Rubber  Septum  through 
which  a Gas  Ullage  Syringe  passes  to  introduce  a precise  loading  of  air 
ullage  (voIubm  percent,  STP' . Liquid  propellant  is  introduced  into  the 
Pneuaatic  Load  L.P.  Cylinder  through  a fill  port  equipped  with  a tube  fitting 
and  valving. 
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hydraulic  pressure.  to  200  kpai?  ar.  mount^d^  ^’'  ™gg!,dizad  withstand 
to  record  preesure-tiae  histories  durino  . h*  h'P'  ComPr«**ion  Chamber 

transducer  is  located  in  th.  pr.a.u^  °"  t9*tm  0na  pressure 

Starter  Charge  Chamber  and  the\  P ^ pa8sa<?*way  between  the 

driving  p-t  history  e.n.ed  bv \h.  h ,i0"  Chamber  to  measure  the 

transducers  and  three  light  Senior  a^bM01*11"*  Charge  1 Thr”  pressure 
1.0  inch  on  center,  along  th.  Ccxrpr.a.ion  cJLSir  ^ Wt  *"d 

ch^Ur,  «?^Ui1.tic^nT*w^t"if^“rS^r.L'j7.C0nPr*"iOn 

‘-1  gjototranalator  light  S.n.or  >.aa„bii.. 

*>— > M0D  1 h*'™ 

Of  such  a light-.en.inq  system^  J.STJ  Compression  Chamber  design,  use 

the  liquid  oropellant  charge  i.  a usef^l^ia  **  ?f  datactia9  ignition  of 
ignition  ut^Iies.  The  design  provides  for  (i)°!!iC  f001  f°r  tha8a  co®Ptession- 
in  the  chaster,  (ii)  . seal  aoain.t  aaay  ins®rtion  *nd  removal 

aansitiv.  Phowtr-i^to"  °f  th' 

l-Ui.,1  aapphir.  trindov  whan  •^Wc^^thTh?*."*0*1’''’'  °f  the 
hot  gaa  environment.  uoiacted  to  the  high  pressure, 

Tha  light  aanaor  assembly  has  two 

translator  hou.ing  which  contains  thTactive  f°™f?nentS:  (1>  the  photo- 
(2)  the  high  pressure  window  hou«in„  ■>,  k ^"‘ght  sensin9  element  and 
from  the  high  pressure  (25-50  kp.i)9  L*C  Protects  th#  Phototransistor 

•«*  ss.*55s*s2ssa,,s:Jirlir  i*»qth  °f  *-m  ^ «*>■ 

connector.  The  opposite  and  of  th.  f ting  in  a Mlcrodot  fumale 

tf*  phototr“,i*^r  l.n.  ccMd  dor  proJ“^“'r  hou,in'!  i5  °P*"  »ith 

thr.^*toi,is„”i*n*“rtn;0i”s: ha;;;i;.!,iu°  “ch  - ^ w w. 

V2  inch  wench  cl  act  for  tiqhtanino  u#ln9'  Th*  window  housing  his 

«■— ? to  .ccpt  th.  „d  TToTIT,  i"'""*1  bo~  with 

x«.  Sight  hoi.  to  ad  it  light  Cton  th.  ..‘JZ  P^Stit^istot. 

"“•fs  *iSrfy~3~^  “-400 

hemiapharical-ahaped  lens.  Tha  oD.ti  harmatically  -sealed  case  with  a 
exposed  is  -65  c to  125  r*  Th.  °P**ating  temperature  range  with  elements 
to  fully  saturate.  The  outer  dimensiL**9*1^*3  9mW/cm'  of  incident  radiation 
in  length  and  0.08  inch  diameter.  The  collect*  pho’Jotrar‘sistor  are  0.6  inch 
«•  accomplished  by  mean.  of  two 'fine  coppi*^^  Connections 

brass  femala^cr^d^t^^ictor^fter^iT3^^ ®Ctor  Wir*  iS  soldered  fc°  a 
TVo  Ml  Teflon  alaeve  collars  i/sul.^^L  *ith  heat  shrink  tubing, 

transducer  housing.  The  brass  oroundin  k em*le  Microdot  for  the  photo- 
to  the  housing  when  th*  phototrmadUf^ ^ Passes  the  emitter  lead 
transistor  element  is  potted  into  the  L . f lnally  assembled.  The  photo- 
silicone rubber  RTV-lOfh  lnt°  th*  housin9  using  General  Electric 
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A lock  nut  is  tightmad  to  hold  the  photo  transducer  in  place.  Tha 
•a^phire  window,  which  atuurti  0.125  inch  diameter  by  0.200  inch  long  and 
is  manufactured  by  Adolf  Mailer  Go.,  haa  optically  poliahed  faoaa  for  optimum 
light  transmission.  Tha  window  ia  placed  into  tha  high  pressure  housing 
from  tha  pressure  side  and  ia  sealed  by  naans  of  a Parker  Mo.  S-17M-R5S2-90 
•0"-ring  in  a dove-toil  groove  in  the  housing  face. 

The  relative  spectral  response  of  the  IS  400  photo translator  ia  shown 
in  rigure  28  , with  tha  response  of  the  human  aye,  output  of  e 2870  X 
Tungsten  source,  and  output  of  GaAs  sources  displayed  tor  coeperison. 

4.3  System  Operation 

Following  tha  complete  assembly  of  the  L.P.  Compression  Chamber  and  the 
Pneumatic  Load  L.P.  Cylinder,  the  Starter  Charge  Chamber  is  fitted  with  an 
electric  primer  mad  loaded  with  the  prescribed  smokeless  powder  air.  The 
Starter  Charge  Chamber  and  the  L.P.  Compression  Chamber  ara  than  joined 
together  with  the  Chamber  Coupling  Collar  on  the  test  stand.  The  liquid 
propellant  filling  procedure  can  now  begin. 

Tha  Phemutic  Load  L.P.  Cylinder  and  tha  L.P.  Compression  chamber 
are  first  bled  to  remove  residual  unwanted  air.  This  ia  accomplished 
by  first  positioning  tha  Separator  Piston  and  the  Projectile  Piston  with 
ends  on  either  side  of  the  Compression  Chamber  fill  port,  thereby  forming 
a 5/32  inch  wide  passageway  for  flow  of  liquid  from  fill  port  to  Bleed 
Valve.  The  Poppet  Valve  is  than  manually  held  in  tha  open  position  as  a 
gravity- fed  supply  of  liquid  propellant  flows  into  the  Pneumatic  Load  L.P. 
Cylinder,  through  the  Poppet  Valve,  Flow  Guide,  fill  port,  the  5/32  inch 
slice  of  Compression  Chamber  bore,  and  out  the  Bleed  Valve.  In  this  nanner 
the  entire  system  including  the  initial  bore  volume  between  the  two 
Compression  Chamber  pistons  can  be  flushed,  with  the  Poppet  Valve  released, 
free  to  return  to  its  normally  closed  position,  a prescribed  volumetric 
loading  of  gas  ullage  is  introduced  into  the  Pneumatic  Load  L.P.  Cylinder 
with  the  gas  syringe.  The  displaced  liquid  propellant  volume  is  permitted 
to  flow  out  through  the  Isolation  Valve.  All  valves  are  then  closed.  The 
liquid  propellant  reservoir  is  removed  from  the  test  area,  all  electrical 
connections  completed,  and  the  safety  enclosure  closed.  Upon  enabling  the 
instrumentation  to  record,  the  system  is  ready  to  fire. 

The  ignition  circuit  which  ultimately  initiates  the  electric  primer  Is 
energized  by  first  closing  the  Ignition  Circuit  Safety  Interlock  key  switch 
and  then  the  Ignition  Circuit  Arm  switch  on  the  control  panel  (Refer  to 
Figure  29  , Firing  Circuit  Schematic.)  This  applies  +500  VDC  to  the  open 
leg  of  the  microswitch.  The  firing  voltage  is  not  sent  to  the  electric 
primer  until  the  microswitcn  is  closed  at  the  completion  of  dynamic  loading. 

A liquid  propellant  load  circuit  is  now  energized  by  closing  the  Load 
System  Safety  Interlock  key  switch.  Refer  to  the  schematic  drawing,  Figure 
29  , for  the  following  description.  The  dynamic  loading  sequence  is 

initiated  by  closing  the  Load  System  Firing  Switch.  This  action  opens 
the  Solenoid  Valvs  which  releases  the  driving  N2  pressure  into  the  Pneumatic 
Load  L.P.  Cylinder  to  accelerate  the  Pneumatic  Piston.  The  motion  of  the 
Pneumatic  Piston  forces  the  liquid  propellant  and  its  gas  ullage  past  the 
Poppet  Valve  once  the  force  balance  establishes  Poppet  Valve  opening.  The 
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Figura  28.  8alativa  Spectral  Response  Char  sctsn  sties  ot  l*  400 
Phototrantiitor. 


liquid  propellant  and  it*  gas  ullage  flow  through  the  Plow  Guide , possibly 
resulting  in  cavitation  depending  on  the  nature  of  the  orifioe , into  the 
bore  of  the  L.P . Compression  Chamber.  The  Projectile  Piston  is  then 
driven  to  the  right  allowing  the  charge  to  fill  the  bore  until  the  maximum 
stroke  is  achieved.  The  contact  of  the  Projectile  Piston  with  the  pushrod 
of  the  Micro switch  Assembly  closes  the  micros* itch  and  applies  the  + S00  VDC 
to  the  Electric  Primer  in  the  Starter  Charge  Chamber.  The  resulting 
pressure  generated  by  the  uambustion  of  the  starter  charge  mix  is  sensed 
by  the  Separator  Piston  which  is  free  to  move  in  the  Compression  Chamber 
bore.  As  the  pressure  in  the  Starter  Charge  Chamber  rises,  ths  Separator 
Piston  accelerate*  compressing  the  liquid  propellant  charge  along  with 
its  bubble  distribution.  The  nylon  portion  of  the  composite  Separator 
Piston  expands  radially  outward  to  the  bore  to  prevent  hot  gases  from 
reaching  the  liquid  propellant  charge.  This  is  a high  pressure  seal . Ms 
rely  on  "0" -rings  fitted  to  the  Separator  Piston  for  low  pressure  sealing. 
Ooaibustion  gases  from  the  starter  charge  are  continuously  being  vented  through 
the  0.053  inch  diaamter  orifice  in  the  Vent  Plug.  If  the  compression 
cycle  does  not  produce  an  explosion  of  the  liquid  propellant  charge,  the 
Separator  Piston  will  return  to  its  original  position  as  Starter  Charge 
Chamber  pressure  declines.  If,  however,  the  liquid  propellant  charge 
undergoes  runaway  mention  at  sons  point  in  ths  compression  cycle,  the 
Projectile  Piston  will  shear  through  the  alixiinum  Shear  Disc  whan  the 
bore  pressure  exceeds  65-70  kpsi.  As  ths  Projectile  Piston  emerge a from 
the  bora,  the  combustion  gases  from  the  liquid  propellant  charge  are 
vented  out  though  the  bore  hole  in  the  Disc  Detainer  bolt.  The  venting 
is  so  rapid  that  a positive  pressure  gradient  between  the  Starter  Charge 
Chamber  and  the  L.P.  Compression  Chamber  drives  the  Separator  Piston  from 
the  bore  also.  Both  the  Separator  Piston  and  the  Projectile  Piston  must 
be  replaced  in  the  event  of  a liquid  propellant  explosion. 

4.4  The  Data  Acquisition  System 

Princeton  Combustion  Research  Laboratories',  Inc.  data  acquisition  system 
includes  a Biomation  Waveform  Recorder,  Model  1015,  Tektronix  Model  7704a 
Oscilloscope,  DATAcap  Coupler  Model  B403,  and  Texas  Instruments  Modal  733 
ASR  Electronic  Data  Terminal. 

The  Biomation  Waveform  Recorder  records  four  simultaneous  analog  input 
signals,  with  resolution  of  each  input  signal  being  one  part  in  1024  (batter 
than  1/10  percent).  The  four  input  signals  recorded  are  : Channel  1: 

Starter  Charge  Chamber  pressure , p£ ; channel  2 : Light  sensor  output  , LI, 
Channel  3 1 Liquid  Compression  Chamber  pressure,  p ; Channel  4:  Liquid 
Compression  Chamber  pressure,  p . The  recorder  is  operated  in  the  Pre- 
trigger  Mode.  This  mode  of  data  recording  permits  the  recording  of  the 
event  prior  to  the  trigger  signal,  if  desired,  and  tha  subsequent  signal 
history.  This  is  accomplished  by  having  the  unit  record  continuously  until 
trigger  is  received,  then  using  the  trigger  to  stop  the  recording  process, 
freezing  the  data  that  was  recorded  prior  to  triggering  in  memory.  Each  of 
tha  four  channels  store*  1024-10  bit  words.  The  sampling  rate  utilized 
is  100  TOjz  (0.01  msec/word),  resulting  in  a data  recording  window  of  10.24 
msec.  'Lie  window  is  quite  adequate  to  capture  the  events  of  interest.  The 


wawsfoi m noortar  ia  opantad  1a  eha  pratrifqac  aada  with  a tricar  daisy 
of  900  uapla  intarrala.  The  triflar  signal  irtiliaad  ia  tha  aatanal  toltaga 
aaaoclatad  with  tha  Starter  Charye  Chamber  pressure  transducer  response . 

Mwn  tha  trigger  circuit  of  tha  waveform  reaordar  detects  tha  asternal 
voltana  signal  of  0.6  V (or  6000  pai) , it  raeoris  900  words  of  ensuing 
data,  tola  provides  for  sufficient  resolution  of  tha  rapid  praaauriaation 
' of  tha  p*t  start-up  history. 


The  stored  data  on  tha  wevofora  rwoordar  can  be  output  to  an  oscillo- 
scope tor  visual  display  and  photographic  recording,  and/or  to  tha 
printer  of  a teletypewriter,  and/or  directly  to  a data  file  on  tha  tian- 
ahaglag  onaputar  (XflN  370-156)  tor  subsequent  output  in  tabular  ton  and/or 
Chinan  plate,  xn  order  to  acoenpliah  tha  two  latter  none  of  digital 
data  display,  a PhThoap  Oouplnr  and  tanas  merman ta  Electronic  Data 


tha  OWBhoap  Couplar , Nodal  640 3,  la  an  Integrated  circuit  unit  specifically 
designed  to  interface  tha  It  naaMon  Mesa  fora  ha  carder  ta  serious  recording 
and  oaaputlng  units,  by  providing  tha  circuitry  and  oahla  naoaaaary  for 
oaanoeting  tha  digital  output  of  tha  wawaton  recorder,  tha  coupler  is  sat 
up  to  accept  tha  digital  data  fra  tha  wawaton  recorder  an  a 
(conand-type)  basis.  Digital  data  ia  straand  free  tha  wave  fora  recorder 
into  tha  Dhthoap  Coupler  and  than  onto  a oeaoatta  tape  an  the  tone 
Instruments  713  w Electronic  Onto  terminal,  whan  tha  entire  data  file 
has  been  stored  an  ngnetic  tape  it  ia  than  sf  is  snail  to  a file  on  «n  IBM 
370-156  tine-sharing  oonputar  tor  euhaaguant  calibration,  tabulation,  and 
plotting. 


Functional  to 


mart  aimlant  (Water)  tests 


functional  tests  of  tha  ays  ton  wars  performed  with  water  to  determine 
that  ail  aspects  of  system  operation  rhachad  out.  theae  functional  tests 
with  water  also  provided  a baseline  tor  comparison  of  oompreeaion-ignlton 
sensitivity  studies  performed  with  N06-365  liquid  nonopropellssit.  in  ell 
functional  testa  the  line  pressure  to  tha  Pnararatic  Load  L.F.  Cylinder  was 
300  prig  and  "neat'  liquid  wea  dynamically  loaded  in  tha  Ooapra salon  Taafii 
Two  different  preasuriaatian  p-t  start-up  carves  war r seployad.  Type  C and 
Type  E. 


A verification  of  tha  systan  reaponaa  to  tha  pressure  loadings  ganarata 
by  Type  C and  Type  B starter  charges  was  obtained  by  smercising  the  "E0WR- 
SCEPTEE"  Oonputar  Program  54  for  a liquid  spriag/aess  system  with  physical 
properties,  a.g. , density  and  aowd  spaed,  representative  of  H06-365  liquid 
nonopropellant. 


In  this  simplified  approach,  no  account  is  taken  o*  systan  damping  and  tha 
liquid  is  assumed  "beat".  The  ness  of  tha  liquid  charge  is  lingted  with  tha 
nass  of  the  Separator  Piston,  tha  sum  being  Ml.  The  applied  pressure  produced 


by  tbs  atvtir  sharps  an  tbs  cruse cHwil  nm  of  the  toparstor 

non  results  la  tbs  driving  tores  r,  input  to  tha  Bfiirir  pcsgrar  ft* 
spring  noitat  Ml  ter  tha  Mat  liquid  asaapropallant  chart*  is  obtained 
as  follows i 

ft*  lasa tropic  balk  aadulus  of  tha  liquid  is  fim  byi 


dr  - k di 


mmi.  - (c*0,Bo«-Jd5 

• (1.193  a 10*  oa/aac) 2 (1.391  g/cc) 

- 0.302  a 10u  dyaas/oaz 
or 

k - 7. 314  a 10*  dynaa/os  S K1 

sharo  v • 1.402  eaj 
1 - 1.104  on 

rigors  30  aboos  tbs  predicted  Ports*  la  tbs  liquid  etcnopropallant 
"spring*  ia  rs spaas*  to  tbs  driving  fores  loading  associated  with  startar 
ebargo  proasurlastloa  cars*  Type  S.  Largo  saplitodo  ooclllatioas  ars 
prodletad  ia  tbs  liquid  "apriag”.  ta  fact,  tbs  liquid  ia  ia  taasioa  at 
t - 0.11  none.  of  ooarss,  tha  response  of  tbs  raal  systsa  iacltadas  tbs 
offsets  of  daapiag.  tharsfsrs  tbs  teres,  or  pr scaurs,  seoursioos  ia  tbs 
liquid  ars  not  mtlolpatsd  to  bs  as  violent,  liailarly,  rigors  31 
sbong  tbs  fnHetsd  teres  ia  tbs  liquid  opriag  ia  xeapouss  to  tbs  driving 
tores  leading  associated  with  startar  charge  prossorisstioa  cures  Typo  C. 
Motion  that  tha  predicted  liquid  pressure  response  is  lass  violsat , tbs 
avlituds  of  tbs  faros  eaoureieu*  bsiag  sigaifiesatiy  rsduosd.  Also,  tbs 
liquid  is  slasyo  ia  osapraaaiee,  i.o. , a asgativs  pressure (tension)  doss 
not  result. 


ta  of  tba  teuctieuoa 


. J,  p , c*  23  kpoi  ia 

firing  records  of  tie 


a are  identified  ia  Figures  32  aad  33 

d.  those  figures  display  digit  load 
history  ia  tbs  Starter  Charge  Cbaaber , 
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p (t)  , and  in  the  bora  of  the  L.P.  Ooapression  Chan bar . P.(t)  and  9.(t)  . 
Pressure  tranaduoar  p la  locatad  0.60  inch  fit*  tha  rear  faca  of  tna 
Saparator  Platen  and  la  locatad  2.00  inchaa  downs tra an  of  p,  • tha 
vortical  scala  in  thaaa  figures  ia  20  kpsi/inch,  with  a hi f tad  Baselines 
for  aach  pressure  racord.  tha  oacillation  frequency  of  5000  Ha  in  tha 
p - t racord  ia  asaociatad  with  axcitation  of  tha  fundamental  acouatic 
■Sda  of  tha  a tartar  charge  chaabar  dua  to  tha  conbustion  procaaa.  tha 
liquid  reaponee,  p and  p, , (essentially  a damped  spring-mess  ayatan)  to 
tha  driving  presaufe  p,  associated  with  a tar  tar  Typa  K diapl&ya  rather 
large  amplitude  praaaura  excursions  of  the  seam  fraquancy.  Although  tha 
■a an  paak  liquid  praaaura  ia  aaan  to  be  25  ftpei,  a praaaura  wave  of 
50  kpei  amplitude  ia  generated  in  tha  dynamical  ly-loadad  liquid  oolumn, 
at  tha  location  of  p_.  Motioa  that  tha  nagnituda  of  tha  praaaura  spikes 
in  tha  liquid  oolumn  are  not  as  aavara  for  "ZMDtT  ROM  3*  sinoo  the 
pressure  axcuraiona  developed  by  the  macks  lass  powder  starter  charge 
Type  C are  leas  savers.  This  was  predicted  by  the  systaaa  response  cal- 
culations of  tha  "SOPO-SC*PTI«"  Coda. 

In  both  of  thaaa  .Inert  loading  teats  ayatan  components  and  inatru- 
nantation  functioned  perfectly.  Tha  integrity  of  all  seals  was  estab- 
1 is  had . As  anticipated,  all  light  sensors  renained  an  baselina  for  tha 
duration  of  tha  inert  loading  teat,  indicating  that  coapcaaalon-lgnltlon 
of  tha  dynamically  loaded  neat  water  charge  did  not  result.  Conpreasian- 
ignitian  sensitivity  studies  with  MOM- 365  liquid  nonopropsllar-C  are 
discussed  next,  in  Section  4.4. 


4 * k»-h>  Liquid  d— copellant  tuti 


OHvnMloi'ipitlea  Muitiflty  itullta  of  dynwicallr-lotM 
■04-3*5  liquid  urmopcoffallant  «u«  pogfOK— d with  thb  ippantu*  pco- 
Tiosaly  diieriM.  Iba  psasaurlsstian  rate  to  «hieb  the  liquid  pro- 
pa  llont  charge  with  associated  ullaqa  la  ad)MtQd  is  as  Uportat 
par— tar  la  dataculalag  the  thra ahold  for  runaway  raactioa  ad,  there- 
fora , tbs  Mb  operating  Units  of  start -up  pressor i sat ion.  Tha  netrln 
of  taata  parfwd  is  p— ad  ia  labia  3 . Figaros  34  thru  41 
display  fir inq  records  of  tbs  aanartiwity  taata  parfocaad.  Za  aaeh 
oa as  starter  charge  obaabar  proa  aura  and  two  liquid  props  11  sat  ooapraaaioa 
ebaabar  praaaoraa  ara  diaplayod  as  a function  of  tiaa  on  a atafparad 
hasaliaa  seals.  IN  tfcaawlaa  raspaoaa  of  tbs  pbototraaaiator  llqbt 
aanaors  ia  tha  ongpco salon  ebaabor  bora  ia  aboao  aa  aa  inaart  ia  aach 
figure. 


Firing  taata  ■WHI‘  Kb  MO'S  1,  2 and  4 wars  parfocaad  with  oaat 
propellant  and  a dynarlc  loading  dr i war  praaaurs  of  300  paig . m all 
on—  tha  a—  paak  smarter  obarga  praasora  ana  25  kpsi.  tha  preeour- 
iaatian  rata  was  wnriad  bataaan  taata,  oorraapending  to  those  aaaociatad 
with  start-up  earn*  Type  C , X and  0 raapactiwaly.  rigors  34  displays 
taat  results  of  *MMI  MM  1*.  Iba  raaponaa  of  tha  dynaMlcally- loaded 
oaat  paopallant  cbarga  to  a proaaurlsation  rata  glwan  by  TYPX  C start-up 
charge  is  banign  This  can  ba  aaaa  by  inspection  of  rbo  p-t  history  ia 
tha  liquid  cbarga.  a oonpariaon  of  thoao  racords  with  thosa  aaaociatad 
with  taat  *nw  BOB  3".  and  by  tha  laek  of  rospanao  of  tha  photo  tran- 
sistor light  *h«h; 


figure  35  displays  toot  ra suits  of  KMU  KM  2*.  Iba  raaponaa  of 
tho  dya  Ml  nally-loodod  oaat  props  11  snt  cbarga  to  a pressurization  rata 
giwm  by  TUB  t start-up  eurwe,  a faster  pressurization  rata  than  TOP*  c 
curse,  ia  aa  eaplaaien.  Light  aaaaoc  U da  tacts  ooabustion  ia  tha 
li^id  cbarga  at  approalaataly  0.40  anac  after  ignition  of  tba  starter 
charge  (asset  of  rbaadiar  praaaurisatiao) . fallowed  by  a raaponaa  of  LI 
at  t - 0.70  anac  — than  L3  at  t - 0.45  jwac.  iba  prasaura-tian  histories 
ia  tba  li^tid  cbarga  show  large  seals  oscillations  with  paak  proaauras 
ia  aaoaas  of  40  kpsi . Cbarga  anpllflar  gains  wars  increased  in  snbaequant 
*••**  o*dar  to  capture  tho  onset  praaaurs  peak.  seweral  points  ara 
aurtb  ootiagi  (i)  tha  departure  of  tha  Uquid  pressure  raaponaa  fron  that 
of  an  inart  loading  taat.  i.a. , Figure  32  *XMtr  MM  2”,  occurs  ia  tba 
aaigbhorbood  of  t • 0.50  nose  after  ignition  of  tba  starter  cbarga.  and 
(ii)  at  appruniautaly  t • 1.0  mac  after  peasant laatlon  coast  the  shear 
disc  ruptures  and  rapid  depressor  I ration  of  tba  liquid  propallMt  oom~ 
P«»44lan  obaabar  bora  natlu.  ibis  was  our  first  indication  that  tba 
daai^  of  tba  pnaaua  relief  blow-out  ayatan  would  ausurs  tha  integrity 
°f  tba  onap  r a an  i un  - ip»  i t ion  taat  apparatus  as  a result  of  an  onploaion. 
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Figure  39. 
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rigun  37  diapUqri  tMt  multi  of  1QI1I  MM  4",  an  intimdtita 
pro— urlmatlon  rata  oaaa  corresponding  to  Wl  D starter  charge.  The 
dynamically-lomded  noat  propellant  chary*  omderyoes  delayed  naamy 
reaction,  approximately  4.25  mac  after  starter  chary*  iyaitloa.  fb* 
response  of  tbs  liquid  prays 11 an t chary*  la  essentially  bo* if*  up  to 
this  point.  Ob  fortunately,  no  light  a*asor  output  data  was  obtain* d 
duu  to  tba  delayed  liquid  response . Thia  telaynd  explosive  raapoaa* 
of  tba  liquid  r/ropall— t auyyosta  that  a critical  thro*  be  Id  oandition 
for  aaapr*  salon- ignition  has  boon  achieved.  ter  pro— arlaatiaa  rates 
1 na  than  that  associated  with  TTN  D starter  cbarya*  tba  liquid  reape*— 
ia  benign;  for  faster  prasauriaatiaa  rates,  —pl—ion  result*. 

Figure  42  ia  a s— ry  plot  displaying  tha  corridor  for  aaf*  start-up 
operation  of  tha  dynamically-loaded  —at  propellant  chary*. 

9m  naxt  sequence  of  cmynMim-lgiltion  sasaitlvity  touts  parfornad 
was  directed  at  date  min  log  the  affect  of  deliberately  introduced  air 
ullage  (3.1%  by  vol— a)  on  ' re  tendency  towsrd  explosion.  These  tests, 
testified  as  "80— UK  MM  3"  and  "MMSLS  MM  5"  wars  performed  with  a 
dynamic  lovl  driver  pressure  of  300  pelg  and  pcsssuriMtion  rates  given 
by  starts^  clary*  rmt  C and  TTK  D respectively . Figure  34  display* 
digitised  firing  records  of  *IIMU  MM  3*.  9m  response  of  tha  dynamically- 
loaded  *08-345  propellant  charge  with  its  associated  bubble  distribution  to 
a pressuriaat ion  rata  given  by  TT7I  C start-up  curve  is  benign.  Again  we 
not*  the  lack  of  response  of  the  phototranaistor  light  sen  tore  indicating 
tha  absence  of  oasprem ion-ignition. 

Figure  38  displays  digitised  firing  records  of  *IDMU  MM  5*.  9m 
response  of  tha  propellant  cbarya  with  associated  ullage  tn  * faster 
pressurization  rata  given  by  TYPE  D start -up  curve  is  an  explosion.  Light 
sensor  L3  detects  ooabuutiun  in  tbs  liquid  chary*  at  t ■ 0.90  msec  follow- 
ing ignition  of  the  starter  charge.  9m  delay  time  to  explosion  is  consid- 
erably reduced  by  the  introduction  of  3.1%  ullage  a*  compared  to  the  —at 
liquid  reap—  to  TTP*  0 starter  chary*  pressurization . ,9m  pressure 
generated  in  tba  liquid  coluem  ia  seen  to  ba  in  excess  of  100  kpsi  (charge 
amplifier  signal  saturation) . Rapid  depressurization  of  tha  propellant 
compression  chamber  bore  is  seen  to  begin  1.6  nsec  after  the  onset  of 
pressurization.  Also  notice  that  pressure  transducer  p shorted  during  the 
firing  due  to  tha  excessive  pressure.  This  transducer  we*  scrapped  after 
the  firing  due  to  diaphregn  rupture. 

The  next  group  of  taste  was  performed  with  a dynamic  loading  driver 
pressure  of  500  peig  and  a neat  props  1 lent  cbarya.  At  this  high  driver 
pressure  tba  cavitation  phenomenon  is  anticipated  to  be  more  severe  than 
the  300  peig  driver  case , but  the  resulting  bubbles  (from  cavitation  — d 
residual  ullage  in  liquid)  tend  to  be  mailer  in  dim star.  Figure  40  , 
"BtM*lX  mm  7"  - -hc_-  the  bssiy.  , ■■iionas  or  cue  dynamically- loaded  neat 
propellant  charge  subjected  to  TTP*  C starter  chary*  pressurization  history. 
Again  the  light  sensors  show  no  sign  of  compression-ignition  in  tbs  charge. 
However,  as  show  in  Figure  39  for  -BOBBIE  MM  6-,  explosion  results  when 
the  neat  propellant  charge  is  subjected  to  the  start-up  pressurization 
history  given  by  a TYPE  0 curve. 
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Fiqura  42.  Corridor  for  Safa  Start-Up  operation  of  a Dynamically 
Loadad  Naat  MOS-365  Liquid  Monopropmllant  Charm. 
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Xn  «il  of  th*  ibovt-nftmcad  compress ion -ignition  toots,  tho 
dynamically- looted  liquid  propellant  charge  dioplayod  a benign  response 
to  tha  prossuriaation  rata  givan  by  startar  charge  TTR  C.  Zt  was  fait 
that  explosion  could  bo  produced  by  compression- ignition  by  reducing  tho 
dynamic  load  driver  pres aura  to  150  peig  thereby  increasing  tha  aim 
bubble  diameter  distributed  through  tha  liquid  propellant  charge . fho 
results  of  tho  Plow  Visualisation  studies,  discus sod  in  taction  2.3, 
indicate  that  large  bubbles  are  created  with  this  low  driver  pressure . 

Zt  is  our  experience  that  larger  diameter  bubbles  are  more  sensitise  to 
compression-ignition  than  Mailer  diameter  ones.  Therefore  *1011121  MM  8” 
eo\*  conducted  with  a liquid  chary*  containing  3.1%  by  volume  air  ullage. 

As  anticipated,  explosion  anausd.  Light  sensor  L3  detected  combustion 
in  tho  boro  of  tho  liquid  propellant  compression  chamber  at  t ■ 0.70  ease 
after  tho  onset  of  starter  charge  pre sour last ion.  Urn  pressure-ties 
histc  ry  in  tha  liquid  charge  indicates  severs  pres 'Jure  peaks  (p  > 100  kpsi) 
at  approximately  1.00  esse  after  starter  charge  ignition.  The  severe 
overprsssv re  at  axial  location  P3  resulted  in  the  rupture  of  another 
transducer  diaphragm. 


Sensitivity 


Lon  end  Discussion  of 
Suudis s 


Tn  considering  whether  rapid  compression  of  t liquid  monopropsllsnt 
charg*  with  associated  ullage  can  lead  to  secondary  ignition  sites  (hot- 
ignition)  ard  resulting  explosion  of  the  charge,  w*  have  to  consider 
the  detailed  processes  that  may  initiate  a runaway  exothermic  chemical 
decomposition  reaction.  Heat  produced  by  the  coaqpression  process  in  the 
neighborhood  of  a collapsing  bubble  is  believed  to  th#  triggering  agent, 
not  the  pressure  itself. 


When  th#  liquid  monopropellrnt  is  rapidly  compressed  to  c a.  50  kpsi  an 
estimate  of  tho  -siaperature  rise  shows  that,  if  th#  liquid  is  debubbled 
(de-aorsted)  and  if  the  compression  process  is  taken  to  be  isantxopic,  it 
can  be  as  much  as  10  to  20  8c,  with  some  uncertainty  in  th*  calculation 
due  to  la  -Jc  of  orecisa  data  on  compressibility  and  thermal  expension 
coefficient,  among  other  properties.  A completely  bubble  free  monopro- 
pellant is  therefore  unlikely  to  explode  under  sudden  compression. 

However,  with  one  or  h»re  saall  bubbles  in  the  liquid,  a "hot  spot" 
can  be  generated  by  the  adiabatic  compression  of  the  bubble,  either  air 
or  vapor,  and  if  the  bubble  in  initially  large  enough  and  if  the  compres- 
sion i ocess  it,  rapid  enough,  the  sharp  local  rise  in  temperature  can 
cause  the  bubble  to  retain  its  heat  long  enough  to  initiate  exothermic 
chemical  reaction.  The  resulting  rapid  gas  generation,  which  can  lead  tc 
rupture  of  a gun  chamber  or  to  a gun  breech  blow  if  it  is  rapid  enough 
or  if  the  volume  expansion  does  not  acconmodate  the  gas  production,  leads 
to  rupture  and  expulsion  of  the  specially  designed  shear  disc  assembly 
from  the  liquid  propellant  compression  chamber  bore  in  the  PCRL  experiraerts. 
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lbs  following  MMriatt  th*  rtaalta  of  tho  .onpri aaloo-1  jnitiun 
studios i 

(1)  For  pn— tintioi  ratM  loo*  Ohs  20  kp«l/baoc  • dywicaUr 
1 — -«  aaat  ni-MS  liquid  aoMfctfailMt  dm<a  1*  baalga.  ter  pm 
■urfaatioa  ratal  is  tuMaa  of  this  tfemhald  ooailtloa,  aaploataa  of  tbs 
liquid  mapm^laM  ebarqa  go  ml  to.  idithir  iv  bo  mot  or  oaataiainp 
paqoaafaUy  iatr^md  diatrlbatad  alloy . ihit  establishes  tbo  baas  In 
of  safe  atarttq  for  practical  dyuaic  load  LK  ayataaa. 

(29  tbo  severity  of  tbo  oaplooioo  oaad  ooooeiotod  doiay  tiao  to 
taploaiio  ia  o om  flood  liqaid  aaoepmallaat  oharpo  appoan  to  bo 
pcopcrtlaaal  to  tbo  opawo  doaaoity  of  tbo  bubbloo. 

(3)  for  o given  proaauriiatiea  roto , co.  20  hpei/nmc,  babble  oioo 
ploys  o dominant  roll  ia  dotormlnlng  tho  tsdancy  toward  explosion.  Mna 
tbo  omb  bubble  vlta—tor  ia  iaorooood  by  aa  ardor  of  aoplhidi , firm  10"* 
inch  to  10~*  inch,  tbo  thro  oho  Id  for  ruMsay  reaction  shifts,  tbo 
dynmloslly- loaded  liquid  aonopcops lloat  charge  is  no  lonqnr  booiqo  dan 
subjected  to  jrsoowtiootloo  rofeoo  of  20  kpoiytooec  tear  distributed  bubbloo 
of  asn  dimotor  10  1 inch.  This  io  of  practical  significance ■ Lino  (fill) 
praumn  ia  a dynastic- load  rapid-firs  UK  cob  bo  controllod  oo  oo  to 
product  tbo  assignor  bubble  brook  up.  Also  volumetric  ooatant  of  residual 
ullage  cao  probably  bs  minimised  to  tho  raago  omoidorod  ia  tho  afttasan- 

t lowed  toots. 

(4)  Mo  detonation*  turn  occurred  in  tho  caapreaalau-iqnitian 
oonaitirity  stadias. 

(5)  xt  io  apparent  that  explosion  doos  not  occur  daring  tbo  firvit 
proosnro  polos  an  start-op.  Stornaway  exothermic  nactioc  onauoa  os  sau'ly 
•o  0.50  aooc  after  pcooaur ii\at ion  onset  to  as  lots  oo  4.00  nose.  This 
mggoots  that  bubble  collapse  tim  compered  to  hoot  rotoation  (nr  escape) 
tim  during  tbo  preoeuriaatioK  cycle  ia  key  to  understanding  ooepcssauon- 
igaition. 

(6)  Nanapsaant  of  tbs  start-up  pressurization  history  is  of  parcaount 
importance  from  the  standpoint  of  compression -ignition  sensitivity  analysts. 
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5.0  Ifcwiiml  Appromcfc  to  tmitimtlo*  of  Bmplomicm  by  Dynamic  Compression 

SaSEaBte 

5.1  Pam  ■wothmtloml  Situations  to  bo  Oonsiasrod  Theoretically 

A oa^wfhulw  theoretical  aMlrti*  of  tho  Mniitlvlty  probloa  duo 
to  bwtablm*  should  toko  into  account  at  looot  four  dlfforont  physical 
(itwtioH.  0m  is  tho  caoo  of  tho  bufablo-froo  or  noat  liquid.  This 
sight  bo  aonsHsrsil  tho  bsooliao  ossa.  As  a result  of  ita  nan-soro 
oatproosihility , tho  liquid  taut  up  onorqy  in  a rapid  compression,  and 
this  onorqy  nani foots  i too If  as  boat  or  tsqi orators  rioo.  This  rloo  ia 
an  tho  outer  of  a fmw  ton  a of  dogmas  C,  not  onouqh  for  soot  propal  lasts 
to  ossos  a rssoHsy  reaction  ia  the  milliseconds  of  gun  cycle  tins,  however , 
a oontalnsr  of  aaat  liquid  can,  if  tho  applied  rate  of  pressure  rise  is 
tailored  appropriately , (o.g. , a sharp  r empress  Inn  pulse  to  plateau) 
wriUto  during  tho  tramsioe  of  s reflected  expansion  wove,  such  cavitation 
bubbles  can  contain  air  drawn  out  of  solution  as  well  as  tho  propellant 
vapor.  It  requires  detailed  fluid  mechanical  analysis  plus  consideration 
of  finite-rate  cavitation  theory  in  ordoi  to  estimate  theoretically  the 
time  delay  required  for  the  appearance  of  cavitation  bubbles  of  a sirs 
sufficient  to  initiate  the  runaway  reaction . This  question  of  a safe  no  i- 
osvl bating  pressure  rise  profile  has  not  been  investigated  for  a neat 
liquid  to  the  staoe  of  deciding  just  what  the  safe  limits  are,  to  our 
knowledge. 

A second  case  of  importance  is  that  of  undistributed  sir  ullage,  a 
large  cavity  or  "bubble" . it  can  be  shown  by  analysis  that  a cavity  of 
as  little  as  1%  of  a large  caliber  LP  gun  chamber  volume  constitutes  a 
largo  bubble  in  the  context  of  this  theoretical  discussion,  tha  - is,  ora 
whose  collapse  time  is  comparable  to  the  applied  pressure  rise  tine 
and  whose  hast  retention  time  ia  also  long  compared  with  the  pressure 
onset  time  . in  other  words,  a large  bubble  is  one  extreme  cf  the  hot-spot 
initiatior  theory,  and  undistributed  ullage  usually  fits  that  case. 

The  third  case  for  theoretical  consideration  is  that  of  the  isolated 
mall  bubble,  initially  perhaps  9.1  to  0.01  cm  in  size,  or  even  smaller, 
m the  gun  chamber  compression  process,  the  bubble  can  "inplods"  to  a diameter 
as  sms 11  as  a few  microns,  creating  pressure  levels  within  the  bubble 
and  in  the  immediately  surrounding  liquid  enormously  higher  than  the  field 
pressure  in  the  LP  chamber,  and  -e aching  temperatures  in  the  range  of 
many  thousands  of  degrees  Xelvin.  this  is  the  case  treated  here.  The 
bubble  is  taken  as  filled  with  both  sir  and  vapor;  we  have  ignored  for 
simplicity  tha  dissolution  of  the  sir  in  the  surrounding  liquid  but  we 
have  explicitly  taken  into  the  analysis  the  finite  rate  of  evaporation 
and  condensation  of  vapor  in  the  bubble.  Itie  genera,  conclusion  can  be 
anticipated,  that  a runaway  reaction  can  be  provoked  by  a combination 
of  conditions:  a large  enough  bubble,  a rapid  enough  compression,  and 
• sensitive  enough  propellant.  The  formulation  of  the  equations  is  such 
as  to  admit  the  possibility  of  runaway  exothermic  reaction  in  either  the 
vapor  phase  or  the  liquid  phase.  It  should  be  noted  that  monoprope Hants 
composed  cf  water  solutions  of  inorganic  explosive  compounds  that  do  not 
vaporize  cannot  produce  explosive  vapors  in  the  bubble , so  it  must  be 
inferred  that  the  runaway  reaction  in  such  cases  taker  place  o>  1 v in  the 
liquid  phase,  i'f  so,  the  theory  would  suggest  that  the  react  n goes  most 
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rapidly  right  at  tha  fr  iable -gas  interface,  where  tha  temperature  reaches 
the  highest  levels . I tha  liquid  charga  in  tha  gun  chamber  containc  many 
bubbles,  a raaliatic  at  t inept ion , cauaad  by  break-up  and  acattaring  of  the 
unavoidable  air  ullage  introduced  in  tha  filling  process  or  brought  out 
of  solution  by  cavitation,  than  tha  question  can  be  asked,  whether  the 
isolated  bubble  theory  is  applicable.  It  can  ha  shown  theoretically  that 
the  nost  important  affects  of  the  bubble  in  causing  a runaway  reaction 
are  contained  within  a few  radii  of  tha  bubble,  so  that  if  tha  ullage  is 
lass  than  1%,  there  is  a good  chance  that  most  of  tha  bubbles  can  be 
treated  by  single -bubble  thecxy. 

This  brings  up  tha  fourth  cm,  that  of  bubble  clusters.  No  theory  < 

is  offered  hare  for  this  case.  Obviously,  it  can  becona  tha  dominant  case  | 

whan  tha  anount  of  ullage  is  sore  than  1%.  A bubble  cluster  can  lead  to  j 

runaway  reaction  wore  readily  than  tha  isolated  bubble,  since  it  can 
c rents  a situation  in  which  a wall  field  of  propellant  is  surrounded  by 
heat  sources , thus  being  heated  no re  rapidly  to  decomposition  temperature 
levels  than  by  a single  central  hast  source.  But,  unless  the  ullage  is 
large  (which  ought  to  be  corrected  in  the  design)  , the  number  of  such 
clusters  would  be  much  smaller  than  the  number  of  single-bubble  reaction 
sites,  and  sc  clusters  nay  not  be  the  dominant  problem.  However,  the 
cluster  situation  does  deserve  analysis.  It  has  never  been  treated,  to 
our  knowledge . 

5 . 2 Scope  of  the  Present  Work:  Fundamental  Processes  and  Theoretical 
Considerations 

With  the  preceding  discussion  as  background,  the  present  work  was 
initiated  dimply  as  an  antry  Into  the  problem.  Limitations  of  contract 
time  and  very  severs  limitations  of  contract  funding  precluded  more  than 
a definition  of  tha  problem  and  an  indication  of  the  manner  of  solution, 
rar  the  reasons  given  above,  attention  was  focused  on  the  3ingle-bubble 
problem. 

Many  of  the  ideas  that  underlie  the  following  theoretical  formulation 
were  drawn  from  the  extensive  work  that  has  been  done  on  the  theory  of 
bubble  dynamics  aver  since  the  pioneer  work  of  Rayleigh  1 , most  of  it  ! 

with  an  eye  on  application  to  cavitation  damage.  Our  problem  differs  j 

from  the  cavitation  work  in  one  important  respect:  the  liquid  is  an 
explosive  medium  that  can  decompose  rapidly  and  exothermically  when  initiated 
by  a collapsing  bubble  serving  as  a hot  spot.  j 

Before  entering  into  the  analysis,  it  is  worthwhile  to  discuss  the  j 

physical  picture.  We  visualize  a small  air  bubble  at  some  moderate  initial 
pressure,  containing  some  vapor,  perhaps  one  millimeter  in  diameter  or 
even  less,  located  in  the  middle  of  a field  of  liquid  monopropellant.  J 

r At  the  starting  instant,  the  applied  pressure  comes  on,  not  simply  as  a 

i step,  bu  as  an  arbitrary  function  of  time,  rising  to  its  operating  level 

| of  ca.  5u  kpsi  in  a time  of  the  order  of  1 msec.  The  "wavelength"  of  the 

> applied  pressure  is  much  larger  than  the  diameter  of  the  bubble,  so  that  ; 

it  is  legitimate  to  treat  the  liquid  field  around  the  bubble  as  havinq  j 

a uniform  value  of  "pressure  at  infinity",  which  is  taken  to  vary  with  i j 

time  in  the  manner  of  the  applied  pressure.  I ] 
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An  "implosive"  motion  of  the  liquid  ensues  and  the  bubble  starts 
to  contract.  Due  to  the  momentum  acquired  by  the  liquid  in  its  inward 
motion,  the  cressure  in  the  bubble  can  reach  levels  on  the  order  of  a million 
atmospheres.  The  adjacent  liquid  reaches  similar  pressures.  The  temperature 
of  the  gas  in  the.  bubble  can  reach  peaks  on  the  order  of  tens  of  thousands 
of  degrees  Kelvin.  At  the  severe  pressures  in  the  micro-domain  around 
the  compressed  bubble,  the  liquid  must  be  taken  as  compressible.  This 
moderates  enormously  the  bubble  pressure  that  would  have  been  predicted 
on  the  basis  of  incompressible  liquid  theory;  and  correspondingly,  the  gas 
temperature  in  the  bubble  is  moderated  enormously  by  the  introduction 
of  compressibility.  One  of  the  examples  given  below  shews  that  the 
gas  temperature  decreases  by  a factor  of  five  when  the  compressibility 
is  introduced.  A numerical  example  shows  that  as  much  as  95%  of  the 
energy  of  the  implosion  resides  in  the  compression  work  absorbed  by  the 
compressible  liquid,  at  the  moment  of  peak  bubble  pressure, and  only  the 
other  5%  in  the  compressed  gas.  Conduction  of  heat  from  the  hot  bubble 
gas  tc  the  surrounding  liquid  also  moderate*,  the  pressure  rise  and  the 
temperature  rise  during  the  collapse.  Tiiis  non-adiabaticity  is  difficult 
to  predict  in  the  absence  or  information  about  the  motion  of  the  gas  in 
the  t ny  collapsing  bubble.  In  the  present  work,  it  is  assumed  that 
the  gas  is  simply  in  radial  contracting  motion  and  that  the  pressure  is 
approximately  uniform,  i.e. , that  the  bubble  wall  Mach  Number  is  much 
less  than  one.  The  formulation  for  heat  loss  ests  therefore  on  molecular 
heat  conduction  in  a contracting  (moving)  gas.  A similar  formulation 
is  used  to  describe  the  heat  conduction  in  the  su-  ^ -Minding  liquid. 

The  heat  conduction  processes  within  the  bubble  and  outside  the  bubble 
are  incorporated  in  the  gm  eral  energv  equation  of  the  fluid  in  each  case. 

It  will  be  noted  that  each  energy  equation  includes  a ehem: car  reaction 
term  to  describe  the  dependence  of  the  exothermic  reaction  heat  generation 
on  the  local  temperature, that  is,  the  formulation  allows  the  possibility 
of  reaction  in  both  the  gas  (vapor)  and  the  liquid..  Associv'  . d with  the 
question  of  the  flow  of  heat  from  the  bubble  into  the  surrounding  liquid 
is  the  rate  of  vaporization  of  the  liquid  to  -o*-m  vacor  in  the  bubble. 

In  the  present  work,  finite-rate  vaporization  kinetics  ir  included.  In 
some  o tner  published  work,  a spec; ally  in  the  cavitation  theory  literature, 
it  has  been  commonplace  to  take  the  vapor  density  in  the  bubble  as  being 
in  thermodynamic  equilibrium  with  the  surrounding  liquid  surface  temperature. 
We  have  reviewed  tbe  supporting  order  of-magnitude  rguinents  and  we  axe 
not  convinced.  Accordingly,  tlw  fin  Lit— rate  evaporation-condensation 
equations  are  introduced. 

One  of  the  most  interesting  features  of  the  problem  arise”  frorn 
consideration  of  the  tine  scales  .inherent  in  the  problem.  Five  r vn  be 
identified:  (1)  The  induction  time  for  cheidcal  reaction  in  the  vapor 

phasv , determined  at  the  high  temperature  at  the  moment  of.  peak  bumf  le 
pressure;  (2)  The  corresponding  chemical  reaction  time  in  the  liquiu  phase; 
(?)  The  duration  of  heat  retention  in  the  bustle  (heat  conduction  time;  ; 

(4)  Tne  rise  time  of  the  imposed  Dressure  'brain.  ";  (5)  The  colla  ;e  time 
of  the  bubble  in  response  to  the  imposed  pressure  field. 

It  *s  possible  to  anticipate  some  of  the  theoretical  result If 
the  chemical  reaction  induction  tire  is  ici  7 cocrured  with  " e :v>at  reter  tioc 
time  (i.e.,  slow  kinetics  and  tiny  hubbies)  there  is  a uo od  .nance  that 
the  -unaway  reaction  can  be  avc lied . 'This  effect  ■-*"  ieaorstratec  1; 
the  PCRL  ernpressicr.  ir.itiat:  _r.  tests  . tcrtec  i.r  L~~?  . f ;; 
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Consideration  of  the  relative  times  of  bubble  collapse  and  pressure 
rise  leads  to  an  interesting  theoretical  problem  that,  as  far  as  we  :an 
find,  has  not  been  treated  in  the  literature  but  has  a very  important 
bearing  on  our  problem  of  explosion  avoidance.  It  can  be  shown  theoretically 
that  the  collapse  time  of  a typical  bubble , at  driving  pressures  on  the 
order  of  1 or  more  kpsi,  is  only  on  the  order  of  a few  microseconds. 

On  the  other  hand,  realistic  pressure  rise  times  (an  LPG  system  or  a 
warhead  explosive  subjected  to  setback  forces)  are  generally  cn  die  order 
of  one  or  more  milliseconds.  Thus,  there  are  hundreds  mud  paihaps 
thousands  of  cycle3  of  bubble  collapse,  rehound,  and  collapse  again  within 
the  rise  tins-  The  heat  retention  time  is  generally  also  long  compared 
with  the  coll -tpse  time.  This  means  that  each  successive  collapse  intro- 
duces more  arc*  more  compression  heat  into  the  bubble,  so  that  the  danger 
point  for  nr  away  reaction  would  be  well  along  in  the  rise  of  pressure, 
perhaps  even  after  the  applied  pressure  has  reached  die  highest  level, 
and  thus  the  point  of  interest  would  be  after  some  hundreds  of  cycles, 
respite  tliie  observation,  pteviour  theoretical  work  has  been  aimed  at 
computing  tf<e  pressure , temperature , anu  reaction  rate  trajectories  of  a 
bubble  during  the  first  cycle  only,  usually  be  a vary  painstaking  numerical 
procedure  cn  a computer.  This  seemed  to  give  plausible  results,  but  only 
because  i:  was  commonly  assumed  that  the  question  of  time  scales  was  an 
unimportar  c one  and  that  the  pressure  rise  could  therefore  be  taken  to 
be  instantaneous.  It  was  thought  that  a judden  junp  is  just  as  good  a 
model  as  any  other  pressure  rise  function.  This  is  not  so.  The  control 
of  pressure  tie  rate  is  probably  the  most  practical  way  to  avoid  explosions. 
The  many-cycle  model  must  not  be  ignored,  but  this  means  that  the  theoretical 
methodology  has  to  be  altered.  It  is  simply  not  feasible  to  compute  at 
length  the  entire  history  of  the  bubble  through  all  the  hundreds  and  perhaps 
thousands  of  cycles  required  to  reach  the  theoretical  explosion  point, 
home  approximate  method  for  summing  all  the  effects  of  successive  cycles 
is  needed.  More  on  this  problem  and  possible  approaches  is  given  in 
Section  9,0  . 

As  pointed  out  at  the  beginning  of  this  section,  the  original  limited 
task  specified  in  the  contract  was  to  develop  a methodology  for  solving 
the  system  of  equation  developed  below.  Although  not  simple,  it  seemed 
straightforward  at  the  time  the  statement  of  Work  was  written.  We  conclude 
from  the  above  arguments  that,  for  a ramp-type  pressure  rise  or  for  any 
other  time  function  of  practical  nature,  a new  method  has  to  be  sought. 
Straightforward  computation  of  the  problem  is  not  applicable  to  real  LPG 
systems  or  to  the  important  problem  of  explosion  avoidance. 
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6 • 0 literature  Analysis 

'6 . l Collapse  of  Bubbles  in  Inert  Liquid* 

The  collapse  of  a cavitation  bubble  was  first  analyzed  by  Rayleigh1 , 
who  addressed  the  problem  of  a spherical,  evacuated  bubble  in  an  infinite 
liquid  domain.  Neglecting  viscosity,  surface  tension  and  beat  transfer, 
as  well  as  liquid  compressibility,  Rayleigh  used  the  momentum  equation 
to  derive  the  equation  for  bubble  surface  motion: 

R d2R/dt2  + j-(dR/dt) 2 - [p(R)  - pQ]/p  (6>1) 

whtre  R(t)  denotes  the  bubble  surface,  P.  ■ const  is  the  driving  pressure  at 
infinity,  and  P(R;  ' » const  is  the  Bubble  pressure.  The  mass  continuity 

for  the  incocpressibleBliquid  leads  to 

u/ U * (R/r)2  (6.2) 

where  u » dR/dt,  which  was  used  in  the  derivation  of  the  bubble  surface 
equation.  The  time  for  collapse  was  shown  to  be 

‘col  ■ c-915  V^Po1'’  t6-3> 

expressed  in  terms  of  the  initial  bubble  radius. 


Considering  the  final  bubble  surface  velocity,  Ray laigh  also  introduced 
the  effect  of  a permanent  gas  in  the  bubble,  undergoing  isothermal  compression 
during  collapse  and  calculated  the  conditions  for  rebound,  as  a certain 
maximal  pressure  is  reached.  In  a concluding  rroaxi,  Rayleigh  naeutior.ed  the 
effect  of  liquid  comprussibility , noting  that  a satisfactory  theory  of 
collapse  sho  Id  include  it  at  the  early  stages  of  development. 


Addressing  the  original  problem  considered  by  Rayleigh,  namely  constant 
bubble  pressure  and  incompressible  liquid,  PoritsJcy2  investigated  the 
effects  ol  liquid  viscosity  and  surface  tension  upon  bubble  collapse.  He 
has  shown  both  effects  to  enter  only  through  the  bubble  surface  balance, 
(viscosity  in  the  liquid  varnishes  under  the  assumptions  of  incompressibility 
and  irrotational  flow' . Incorporation  of  viscosity  at  the  bubble  surface 
results  in  the  balance  of  pressures. 


PB°  = p(R)  + 4y (dR/dt) /R 


(6.4) 


and  the  bubble  surface  equation  is 

R d2R/dt2  -f-  |-(dR/dt)2  + (dR/dt) /R  = 


(6.5) 


where  u is  the  liquid  viscosity,  and  the  added  term  is  clearly  one  of  viscous 
damping. 


The  effect  of  surface  tension  further  modifies  the  balance  of  pressures 
at  the  bubble  surface. 
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p # « p(R)  ♦ 4u(d*/dt)  /*  ♦ 23/*  «6.6i 

B 

so  that  the  forcing  Cara  on  tha  right  land  *ida  of  Eg.  6.5  bncoaai; 

(pBJ  - 20/*  -P0>/o  . (*.?) 

Clearly  surface  tana  ion  ia  important  whan  tha  bubble  radius  is  —all . and 
tanda  to  accelerate  collapee  and  dacalarata  axponaion.  Vlsooaity.  on  tha 
othar  hand  tanda  to  slow  do—  tha  procaaa  of  collages  as  wall  as  rebound . 
Within  tha  ranga  of  para— tars  oonaidarad,  Poritaky2  shows  a reduction  in 
both  collapsa  and  rebotatd  ti—  a by  a factor  -0(10) . Tha  rasults  wars 
obtained  by  m— arical  intagration  of  Bq.6.5  aftar  tr  ana  for— t ion. 

Shu’,  in  a not*  following  tha  analysis  of  Poritaky  obtained  asymptotic 
relations  describing  tha  behavior  of  bubble  undergoing  collapae  or  rebound, 
through  tha  phase  plana  (t)  - R plane)  analysis.  The  results  point  out  tha 
funds— mtal  importance  of  viscosity  and  surfaca  tension  in  tha  probl— 
oonaldorad.  It  should  be  stressed,  however,  that  both  analyses  are  for 
constant  bubble  pressure  and  hence  apply  most  closely  to  the  cavitation 
bubble  probl—  and  not  to  tha  case  where  tha  cavity  contains  permanent 
gas,  in  which  case  tha  driving  pressure  difference  is  ti—  dependant,  and 
oscillatory  collapse-rebound  motion  is  anticipated,  further,  the  consider- 
ation of  compressible  liquid  would  have  necessitated  the  incorporation  of 
viscosity  in  the  liquid  mo— ntun  equation,  resulting  in  a such  higher 
degree  of  analytical  complexity . 

The  formation,  expansion  and  collapse  of  vapor  bubbles  in  a cavitating 
liquid,  upon  passage  over  an  ogive  shaped  body  in  a water  tunnel  wet*  the 
subject  of  an  analysis  by  Plesset.1*  This  work  combines  experimental  observa- 
tions, involving  pressure  measure— nts  along  the  body  and  high  speed  photo- 
graphy of  the  bubble  motion,  with  analysis  of  the  bubble  surface  equation. 

The  following  assumptions  were  —ployed  in  the  analytical  model:  (1)  incom- 
pressible liqu-'d,  (2)  pure-vapor  bubble  (no  air  or  permanent  gas)  , with 
PB  ■ Pq  * Pv  * const,  where  p is  the  vapor  pressure,  (3)  the  bubble  radius 
is  small  compared  to  the  length  scales  ar  • -iated  with  longitudinal  and 
radial  pressure  gradients,  (4)  bubbles  are  c^-i#d  with  the  local  (approxi- 
mately constant  axial)  liquid  velocity,  i.e.,  no  velocity  slip  between  the 
cavity  and  the  liquid,  and  (5)  bubble  is  outside  of  the  boundary  layer.  These 
assumptions  allow  consideration  of  a vapor  bubble  within  an  "infinite"  but 
time -dependent  pressure  field,  with  pQ  * pQ(t)  determined  by  the  observed 
(stationary)  pressure  field  and  the  longitudinal  liquid  velocity  component. 
Thus  the  relationship  of  Eq. 6.1  has  been  used,  with  the  prescribed  forcing 
term 


P(R)-  P0  Pv  ~ 2cr/R  - pQ(t) 

p = p 


(6.3) 


which  includes  the  surface  tension  effect.  In  addition  to  the  aforementioned 
Pq(t)  variation,  R_,  the  maximal  bubble  radius  observed  experimentally, 

— — zr.  ir.puw ; solutions  were  octainea , and  comparison  with 

experimental  data  regarding  R versus  t was  good. 
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rt should  be  mentioned  that  at  tha  pressures  cans  ids  rad  (p  up  to 
8 pais)  and  tha  typical  bubble  radius,  R - 0(0. l)in,  tha  characteristic 
coliapaa  time  is  -0(1) msec,  and  co^>araBli  to  tha  tins  ovar  which  PQ(t) 
varies,  interestingly,  tha  plotted  results  show  tliat  tha  instant  of 
— «is«i  radius  (full  expansion)  lags  appreciably  behind  that  of  minimal 
p (t) , and  occurs  somewhat  prior  to  tha  attainment  of  navi as 1 pressure. 

However,  no  attempt  was  aada  to  calculate  the  pressure  field  about  tha 
bubble,  *ad  tha  extant  to  which  it  inflt^sncas  tha  pressure  profile  about 
tha  iassrsad  body. 

Several  important  aspects  of  cavitation  bubble  dynamics  ware  discussed 
in  a further  study  by  Pleseet.9  Actual  solutions  ware  not  included  in  this 
work;  rather,  the  physics  of  the  cavitation  problan  ware  illuminated  through 
comparison  of  characteristic  scales . Dynamics  of  vapor-bubbles  and  of 
bubo las  containing  permanent  gas  are  treated  separately,  and  tha  (secondary) 
importance  of  rectified  diffusion  in  the  gas,  (whare  aasa  exchange  between 
the  phases  is  considered)  was  pointed  out.  The  questions  of  uniformity 
of  tha  temperature  field  within  the  bubble,  and  limiting  asymptotics  of 
i so  thermal  (at  tha  low  and  the  high  frequency  i imits  of  bubble  oscillations ) 
and  adiabatic  ( intarsmdiats  frequencies)  bubble  processes  ware  discussed  in 
terms  of  the  thermal  diffusion  scales. 

Tha  question  of  stability  of  the  spherical  bubble  shape  was  addressed 
a a well;  using  linearised  small  perturbation  reasoning,  it  was  shown  that 
the  stability  results  for  a planar  phase  interface  are  reverted  for  the 
case  of  a spherical  cavity,  namely,  that  if  the  surface  is  accelerated 
outward  (into  the  liquid)  the  surface  shape  perturbation  term  approaches  a 
constant  or  aero,  while  it  tends  to  grow  out  of  bounds  in  an  oscillatory  manner 
if  acceleration  is  in  the  direction  of  the  gas  (as  in  collapse) , and  R -*•  0. 

The  question  of  stability  of  tim  spherical  bubble  shape  was  addressed  as 
well.  Using  linearized  smell  perturbation  reasoning,  it  was  shown  that  the 
stability  results  for  a planar  phase  interface  are  reversed  for  the  case  ot 
a spherical  cavity,  namely,  that  if  the  surface  is  accelerated  outward  (into 
the  liquid)  the  surface  shape  perturbation  term  approaches  a constant  or 
zero,  while  it  tends  to  grow  out  of  bounds  in  an  oscillatory  manner  if 
acceleration  is  in  the  direction  of  the  gas  (as  in  collapse) , and  R -*  0. 

This  re  a;  on on g is  repeated  with  greater  analytical  detail  in  the  review  by 
Plesset  and  Prnsperetti 6 leading  to  the  conclusion  that  collapsing  bubbles 
are  in  general  unstable  (converging  streamlines  in  incompressible  flow)  and 
would  result  in  shape  distortion  and  subsequent  breakup.  Bubbles  containing 
sufficient  amounts  of  permanent  gas  (relative  to  liquid  vapor)  would  tend 
to  be  more  stable  in  this  respect,  on  the  premise  that  they  might  not  reach 
very  small  radii  at  the  fully  compressed  state. 

These  observations  have  led  to  the  statement  that  the  onset  of  insta- 
bility during  violent  collapse  of  vapor  bubbles  occurs  before  the  effects  of 
liquid  compressibility,  viscosity,  surface  tension  or  heat  transfer  become 
important,  thereby  justifying  analysis  of  the  events  leading  to  bubble  breakup 
through  the  classical  Rayleigh  equation.5  However,  two  major  carats  should 
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be  rilud  r •yarding  tha  proceeding  stability  results > (1)  The  linearised 
an*  lysis  naturally  breaks  down  as  amplitudes  of  the  pertur national  harmonic 
components  erceed  saw  small  valita  compared  to  R(t).  Thus  the  final  reacts 
remain  unrsaolvad,  since  nonlinear  coupling,  »»h  ai  lare  variations  of  R 
itself  during  collapse,  viscous  effects  or  heat  transfer  could  become 
important  and  counteract  to  indefinite  amplitude  growth.  Yet  the  vary  same 
effects  were  intrinsically  neglected  to  facilitate  linearisation.  (2)  The 
effect  of  liquid  compressibility  we*  precluded  in  tin  analysis,  yet  it 
reams  to  beoome  important  at  least  at  the  asms  extents  of  bubble  compression 
as  those  associated  with  instability,  for  driving  pressures  p » 1 atm  are 
considered.  Admittedly,  the  convergent  efcreeml loss  during  ocepressian  would 
enforce  s net s-s table  process  at  bast,  but  the  (particularly  sensitive)  high 
pressure  region  near  the  bubble  surface  would  tend  to  be  stabilised  by 
0p/3t  ji  0,  as  sane  of  tha  liquid  anas  displaced  Inward  is  taken  up  in 
accumulation. 

Moat  treatments  of  liquid  compressibility  were  associated  with  the  forma- 
tion end  propagation  of  hydrodynamic  shocks  in  the  liquid,  and  hence  relate 
to  processes  following  complete  oollepee  of  vapor  bubbles  or  underwater 
explosions,  ^n  tha  latter  context  the  subject  wee  extensively  treated  in  e 
book  by  Cole7,  who  summarized  the  knowledge  accumulated  from  re as arch  prior 
to  and  during  World  War  ZI.  An  exceptionally  sue  essful  approximation  due 
to  Kirkwood  and  Beths*  has  been  advanced,  under  tuts  assumption  that  tha 
entropy  change  in  tha  liquid,  9S  • 0,  (unlike  the  enthalpy  change,  density 
change,  pressure  jump,  etc.).  This  leads  to  the  explicit  expression  of  the 
conservation  of  total  enthalpy  along  an  outgoing  (radial)  characteristic  line, 
viz.  , 


r(h  + uV?)j  ■ const. 

where  I denotes  the  characteristic.  This  approximation  forma  tha  basis  of  a 
theory  by  Gilmore  *,  and  was  latsr  coopered  favorably  with  numerical  solutions 
of  the  compressible  equations  of  motion  (concerning  continuity  and  momentum, 
neglecting  viscosity  and  heat  transfer)  by  Hickling  and  Plesset.10  These 
- athors  considered  the  problem  of  a bubble  containing  permanent  gas  under- 
going compression  and  their  rebound.  Upon  rebound,  a compression  wave  is  3 ant 
into  the  liquid,  steepening  into  a weak  shock  within  a distance  comparable 
to  the  original  bubble  radius  R , and  with  decreasing  strength  due  to 
spherical  dilatation.  Numerical  solutions  were  obtained  by  the  method  out- 
lined  by  VonKeumann  and  Richtmyer  , for  the  equations  in  Lagrangian  form. 

The  solution*,  were  not  pursued  beyond  wave  steepening, as  numerical  instability 
was  encountered. 

Heat  transfer  from  a collapsing  gas  bubble,  provided  the  initial  radius 
is  small  enough,  nas  a primary  effect  on  sonolumine sconce  considered  in  an 
analysis  by  Hickling.12  The  major  effect  of  heat  loss  from  the  compressed 
bubble  gas  to  the  liquid  is  that  of  decreasing  the  emitted  light  intensity. 

The  phenomenon  itself  results  from  subjecting  a liquid  with  dissolved  gas  to 
cavitation  by  a sound  field-  The  luminescence  emanates  from  cavitation  bubbles. 
The  different  eaittance  properties  associated  with  various  gases  were  ex- 
plained in  the  analysis  on  basi3  of  differences  between  their  thenxjphysicai 
properties.  Interestingly,  the  equations  of  sr-  t-—-.  'or  the  brittle  gas  coat  nr. 


(beside  continuity  aad  energy)  a aonantua  aquation  with  an  intarnal 
pressure  gradient,  despite  indications  narta  in  a scale -analysis.  At  the 
initial  pressure  ratios  considered,  p * /p.  — 1/40,  with  p_  • 3 ata, 
it  would  seem  important  to  compere  the  effect  of  liquid  eoaprassibility 
on  tbs  — i temperatures  attained  (the  effect  could  be  cosparable 
to  that  due  to  heat  loss  fro®  the  gas  by  conduction) . 

Beat  transfer  was  also  considered  oy  plasaat  and  Belch1*,  in  their 
treatment  of  Taper  bubble  growth  in  an  initially  uniformly  superheated 
liquid,  end  later  by  Sklmer  end  Benkoff 1 * in  a similar  treatment  with 
general,  spherically  symmetric  tseperature  fields.  Tbs  letter  authors  also 
published  s second  report  car  earning  a binary  liquid,  and  added  & species 
conservation  equation  to  the  liquid  model.  In  both  instances  the  pressure 
within  the  bubble  was  considered  uniform  and  constant , equal  to  the 
vapor  (or  vapor  mixture)  pressure  This  and  aany  other  aspects  of  bubble 
dynamics  ware  treated  ia  e unified  analytical  discussion  by  Hsieh1  *,  where 
e rethsr  conpletw  formulation  of  the  problem  was  derived. 

the  excellent  revise  by  Pleaset  and  Proaperetti*  illuminates  nany  of 
the  physical  proc eases  involved  end  serves  as  s good  source  of  references 
on  the  various  subjects  associated  with  bubble  dynamics . An  additional 
general  explanation  to  the  process  of  cavitation  ia  provided  in  the  book 
by  Batchelor1 7 ♦ 


6.2  Hot  Spot  Initiation  Theory 

Thermal  initiation  of  explosives  can  be  c tainvd  in  several  Banners: 
uniform  heating  deposing  a constant  tseperature  on  the  saxple  surface)  , 
heating  by  shocking  the  explosive  (heating  behind  a propagating  shock  wave 
in  the  substance)  and  by  formation  of  hot  spots. 

In  the  first  instance,  a constant  temperature  it  Applied  tc  an  explosive 
sample  surface.  Depending  on  tha  shape  (slab,  cylinder  or  spherical)  and 
a characteristic  length  of  the  sanple,  as  well  a?  on  the  thennocheaical 
and  thermophysical  properties  of  the  explosive,  runaway  reaction  may  occur 
for  a surface  temperature  beyond  a typical  critical  value.  The  criterion 
is  simply  whether  a particular  temperature  (and  history)  would  allow  a self- 
sustaining  exothermic  reaction  at  some  inner  location.  A study  of  this  form 
of  initiation  was  made  by  Zinn  and  Mader1®  covering  all  three  sampla  qeo- 
me tries  mentioned  above,  and  four  typical  explosives.  Theoretical  analysis 
involves  solution  of  the  nonsteady  Fourier  heat  transfer  equation  with  an 
exothermic  reaction  term  of  first  order  overall,  for  which  solutions  were 
generated  numerically.  Theoretical  predictions  were  compared  with  experi- 
mental measurements  on  a plot  of  time  to  explosion  versus  inverse  surface 
temperature  and  show  good  agreement  in  general  (except  for  TNT,  where  the 
decomposition  kinetics  were  uncertain).  The  authors  identify  three  parameter 
groups  which  determine  the  explosive  behavior  upon  initiation,  along  with 
the  geometry: 


where  a end  C an  tha  thermal  diffusivity  and  specific  haat  of  the 
propellent  respectively,  Q,  A and  t ara  tha  haat  of  ra  act  ion,  tha  pra- 
f actor  and  activation  energy  respectively  and  a la  tha  characteristic 
aasapla  dimension. 

shock  initiation  of  explosives  was  discunaed  racantly  by  Malker  and 
Maalay1’  of  Lewrenoe-Livermore  Laboratory,  presenting  a unifiad  nodal  for 
solid,  liquid  and  gaseous  explosive  subatancaa.  Tha  nodal  dram  froa  a large 
nunbar  of  publiahad  works  on  tha  subjact,  both  thaoratical  and  sapor  imsntal, 
aa  mil  aa  froa  experiments  par  form  d by  tha  authors.  Three  stages  in  shock 
initiation  mra  pr opoaad  as  follows! 

(1)  Production  of  ions,  radicals  and  fraa  ataas  aithar  by  tha  mechan- 
ical ahaar  affect  of  tha  shock  or  by  thaiaal  hasting  behind 

tha  shock  front  through  molecular  collisions,  or  by  coabination 
of  ahaar  and  collisions.  Tha  first  applies  aostly  to  solids, 
where  molecular  notion  is  limited,  tha  second  to  gases  where  tha 
translational  energy  gain  behind  a shock  is  substantial,  and 
tha  coabination  to  liquids  which  have  sore  aoleculsr  fraadoa 
than  solids  but  also  haws  nora  intanaolacular  structure  than 
gaaas.  These  reactions  involve  usually  vary  snail  activation 
energies. 

(2)  Formation  of  statistically  random  reaction  sites  where  exothermic 
reactions  are  sustained.  The  number  of  these  areas  of  concentrated 
reaction  was  found  directly  proportional  to  the  shock  strength 

or  pressure  jump. 

(3)  Input  of  critical  energy  flux  is  required,  for  the  distributed 
exothermic  reaction  sites  to  gr~w  into  an  explosion  or  detonation. 

In  a previous  article20,  the  authors  derived  a correlation  for 

the  critical  energy  flux,  based  on  energy  conservation  relationships 
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where  t is  the  shock  pulse  duration,  Ap  the  mean  shock  pressure 
rise,  p the  density  in  the  unshocked  explosive  and  u the  shock 
front  speed.  Obviously,  this  critical  energy  input  Is  from  the 
shoex  itself  and, except  for  the  relation  to  explosive  density 
and  shock  velocity, contains  no  detail  of  the  actual  reaction 
considered,  unlike  the  surface  heating  model  discussed  previously.18 
In  this  respect,  or!y  a posteriori  estimates  can  be  made  for  given 
shock  conditions  leading  to  explosion  in  experiments,  and  the 
correlation  of  Eq.  6.9  does  not  offer  predictive  capability  outside 
of  a tried  explosive  formulation. 

A simplified  explanation  of  hot  spot  initiat  xi  cf  explosives  was  given 
in  a short  analysis  by  Zinn21,  considering  a spherical  region  of  radius  a 


tt  • unif«—  tlmtwl  fifintutt  T . with  t single,  first  or ter  owrtll 
exothermic.  eh— i rul  reaction.  The  Boost— dy  hast  transfer  aquation  la, 

JT/St  ♦ aV2T  - (Q/Cp)  YA— pt-B/RJ)  (6.10) 

with  Y ™ 1 to  approsi— to  ignition.  A cooling  ti— , during  which  T ^ 
within  tha  aphara  of  radiua  a can  ba  aatahlishad. 

tg  • 0.04  a1.'®  (6.11a) 

dtta  to  diffusion  in  aba— ca  of  —action,  by  lotting  A ■ 0 in  Bq. 
tin  tha  othar  hand,  a haat-up  or  induction  ti—  can  ba  calculated,  during 
v^iich  temperature  r— alna  unchanged  until  tha  an  not  of  rapid  banting  by 
— action, 

asp  (E/R^)  (6.11b) 
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obtainad  by  naglacting  tha— nl  diffusion,  lotting  72t  - 0 in  Eq.  6.10. 

If  t*  - T*,  ignition  la  poaaibla.  This  is  tha  basis  for  tba  calculation 
of  a critical  (minimal)  radius  of  tba  hot  xpot. 
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The  author's  calculatad  result  (numerical  solutions  of  Eq.  6.10  and  a 
species  aquation  where  molecular  diffusion  was  neglected)  agree  well  with 
the  result  of  Eq.  6.12  , na— ly  ignition  was  found  possible  whan  a > a 
However,  conditions  for  ignition  could  not  ensure  propagation.  Indeed1 
Rosen22  has  shown  that  unless  E <<  30  kcal /mole, tha  —action  in  tha  problem 
formulated  above  would  extinguish  due  to  heat  losses  exceeding  tha  gain 
due  to  cxothersicity . 


Although  the  analysis  gives  a clear  physical  picture,  tha  results 
seem  strongly  dependent  upon  the  hot  spot  model,  involving  a sharp  tempera- 
ture fall  at  r - a (outside  of  which  T « T < T.  is  kept)  initially,  leading 
to  abrupt  cooling.  A >sore  realistic  model  may  involve  a Gaussian  initial 
temperature  distribution  with  being  the  peak  temperature  and  the  radius 
a taken  at  sooe  mean  value  between  T,  and  T. . Of  course,  enhanced  analytical 
complexity  can  be  expected. 


A treatment  oi  shock-bubble  interaction  in  gases  (neon  bubble  in  helium 
and  vice-versa)  was  presented  by  Evans,  Harlow  and  Meixner23.  The  details 
of  the  two-dimensional  finite  difference  algorithm  used  were  given  in  a 
previous  report  by  Harlow.21*  The  results  depict  the  process  of  bubble 
collapse  upon  shock  passage,  and  the  deformation  caused  to  the  initially 
planar  shock:  for  the  neon  bubble  in  heli’.nn,  shock  velocity  in  the  bubble 
is  smaller  than  in  ci»e  neiium  outside,  and  hence  the  shock  curves  locally 
backward;  the  bubble  shape  (initially  spherica1'  is  compressed  from  front 
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to  back.  Tho  rtnrM  is  shown  for  hsliun  bubbls  in  nson.  Isotherm  plots 
show  high  tesgmrature  rsgions  within  ths  deformed  bubbles  and  outsids , sftsr 
shock  passage . A high  degree  of  internal  non -uniformity  in  temperature  is 
indicated.  Calculations  with  nitron* than*  (liquid)  containing  an  air  bubble 
apparently  resulted  in  nuswrical  instability. 


Hot  spot  initiation  of  explosives  by  interaction  between  shock  and 
inhaax>geneities  in  the  explosive  (such  as  bubbles)  was  recognised  and  inves- 
tigated in  a preliminary  work  by  Nader. 2 s The  calculations  were  aimed  at 
correlating  available  experimental  detonation  times  in  nitrone thane,  liquid 
TNT  and  single  crystal  PBTN  which  ware  -Odfcjeec  after  the  shock  passage. 

The  one -dimensional  analysis  (formulation  was  not  reported)  was  concerned  with 
two  type*  of  prescribed  hot  spherical  spots,  in  ths  domain  behind  the  planar 
shock,  as  fellows j (1)  thermal  hot  spot,  in  which  tsaperature  was  initially 
elevated  uniformly  over  the  hot  spot  region,  while  density  was  taken  equal  to 
that  of  ths  surrounding  (shocked)  material,  and  (2)  pressure  hot  spot,  which 
in  addition  to  the  elevated  temperature  included  a density  increase.  The 
elevated  temperature  in  both  caeca  was  prescribed  such  that  decomposition 
would  occur  et  a much  faster  rats  in  ths  hot  spot  region  than  in  the  surrounding 
explosive  ( - 1400K  for  nltroam thane) . Thus,  unlike  the  treatment  by 
Evans  et  al,21  no  attempt  was  made  to  assass  the  formation  of  the  hot  spot 
itself.  Several  distinct  equations  of  stats  were  considered,  for  ths  differ- 
ent hydrodynamic- regions , • .g. , for  the  unde tone ted  explosive,  for  the 
detonation  products  (gas) . for  the  mixture  of  detonation  products  and  unreacted 
explosive.  The  work  also  addressed  the  simpler  problem  of  uniform  heating 
behind  e shock  in  homogeneous  explosive.  In  fc*h  instances  good  agreement 
with  the  experimental  data  of  Campbell,  Davis  and  Travis2  6 was  demonstrated, 
regarding  both  induction  times  and  critical  hot  spot  diameters.  Mader  concluded 
that  (1)  the  one -dimensional  (hoaaogeneous)  shock  Initiation  model  was 
successful  and  (2)  that  transport  of  heat  from  hot  spots  under  shock  inter- 
action occurs  primarily  by  hydrodynamic  flow  (i.e.,  convection)  with  thermal 
diffusion  negligible,  as  indeed  indicated  by  the  short  induction  times. 


This  work  was  criticized  by  Enig  and  Petrone27,  who  questioned  the 
validity  of  the  model  predictions  on  grounds  of  (1)  use  of  two  incompatible 
(e«P»P)  equations  of  state  to  describe  the  same  unreacted  explosive  region, 
and  (2)  use  of  unreasonable  values  for  therraophysical  parameters  in  the 
equations  of  state,  which  seemed  necessary  to  obtain  agreement  with  the 
experimental  detonation  delay  times.  The  parameters  in  question  were  Cp  and 
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for  the  unreacted  explosive,  being  the  isobaric  specific  heat  and  the  thermal 
expansion  coefficient  respectively.  Mader25  assumed  these  were  constant 
in  the  process,  i.e.,  c_(T)  ■ C„  , etc.,  and  also  that  the  specific  internal 

? p0  27 
energy  cf  the  explosive  is  e 2 C^T  with  Cy  = Cv  » const.  Enig  and  Petrone 

argue  that  if  -ealistic  experimental  values  were  substituted,  agreement  between 
the  eodel  res  .s  and  experimental  predictions  would  break  down.  The  reply 
Mader  these  comments.  ~ 
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nw  bubble-shock  interaction  problem  was  addressed  again  by  Hader28, 
extending  tha  previous  one -dimensional  modal  to  two  dimensions.  The 
part icla-in -call  method  of  Harlow2’  was  used  for  finite  difference 
approximation,  modified  to  incorporate  exothermic  reaction  in  the  liquid 
(represented  by  Arrhenius  kinetics).  Ir  contrast  to  the  previous  model2’, 
this  work  is  concerned  with  tha  process  of  hot  spot  formation.  To  include 
also  the  possibility  of  shock  interaction  with  (tensity  inhosnogeneities  in 
the  explosive,  the  study  incorporates  an  aluminum  cylinder  in  place  of  the 
bubble  as  well,  the  results  are  generally  in  agreement  with  the  predictions 
of  the  one -dimensional  modal.  For  the  particular  caae  of  tha  bubble,  with 
a shock  strength  of  85  kbar  and  1200k,  temperatures  of  4000K  were  found  in  the 
compressed  bubble.  Mechanisms  of  detonation-failure  mm  discussed,  with  the 
possibility  of  rapid  cooling  by  the  expending  flow  (convection)  from  the 
hot  spot,  me  well  as  interaction  with  rarefaction  waves)  the  second  of  these 
wee  outside  of  the  capabilities  of  the  two-dimensional  model. 

Calculated  numerical  results  for  a variety  of  cases  were  published  by 
Mader  in  a sarias  of  reports 10  1 k , end  e description  of  a ruler ran  reactive 
flow  model  end  finite  difference  algorithm  was  reported  sure  recently1’. 

This  followed  from  tha  conclusion  by  Fader  that  tha  particle -in -cell  model 
had  insufficient  accuracy  for  modeling  of  reactive  flows.  This  improved 
algorithm  was  used  to  determine  cases  of  failure  to  detonate1 11  by  hot-spot 
interaction , and  involves  the  analysis  of  rarefaction  waves  passing  through 
the  reaction  sane  in  nitrasM thane. 

Recently,  experiments  were  performed  at  BRL  with  hydroxy laemonium 
nitrate  liquid  monopropellant1 s,  using  a compression  chamber  in  a blow-out 
gun  configuration.  A single  0.1  cm3  bubble  was  introduced  into  the  mono- 
propellant  under  test.  Pressure  transducers  and  photo  diodes  were  used 
to  detect  the  magnitude  and  duration  of  the  pressure  pulses  emanating  from 
the  compression  charge,  and  to  track  the  possible  secondary  ignition  in  the 
test  propellant.  Such  evidence  was  indeed  found,  but  tha  results  tend  to  be 
admittedly  inconclusive  as  to  the  exact  causo  and  mode  of  secondary  ignition 
by  the  occluded  bubble  'or  other  causes).  Nevertheless  the  results  indicate 
that  further  testing  and  analytical  evaluation  are  required.  It  should  be 
pointed  out  that  pressure  stressing  of  the  test  propellant  does  not  occur 
impulsively  (unlike  the  shocks  considered  by  Mader 2 ’ ' 2 8 ' 3 °~35  but  rather 
gradually,  as  in  the  process  of  initiation  of  tha  actual  propellant  charge 
in  the  gun,  i.e.,  pressure  ramp  durations  ~0(103)ysec. 
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6 . 3 Detonation  h Liquid  Honopropellant i The  Work  by  Stanford 
Research  iaztltute  (SKI) 

Tha  r*Mtrc:'  effort  by  SRI  wu  iuaurix«d  in  two  report* 5 7 ' ! * on  the 
subject  of  detonation  in  liquid  nonoprops Hants.  In  essence,  two  separate 
problsns  were  edcreaewd  in  each  report)  first,  investigation  of  tbs  condi- 
tions for  initiation  and  propagation  of  low  velocity  detonation  (LVD)  > 
second,  on  a micro-icale , the  formation  of  a hot  spot  in  the  liquid  propel - 
ant  by  a collapsing  bubble,  in  the  pressurized  region  behind  the  shock  front. 

the  two  subjects  mentioned  were  not  combined  within  a detailed  analysis 
oi  shock -bubble  interaction  in  the  manner  considered  by  Evans,  Harlow  and 
Maimer* 1 or  b*,  Mader**'1®'*®  1 * (neither  was  mention  made  of  the  latter 
work) , although  essentially  the  same  problem  was  Addressed.  Therefore  the 
review  herein  constitutes  of  two  parts. 

6.3.1  Low  velocity  Detonation  (LVD) 

This  form  of  detonation,  compared  to  the  cooeeon  mode  for  explosives 
termed  hJ gh  velocity  detonation  (HVD) , involves  lower  peak  pressures 
-0(10) khur,  lower  shock  front  propagation  velocities,  -0(2000)  a/s  and 
can  be  sustained  by  partial  reaction  only,  in  the  propellant  region 
behind  the  shock  front. 

Experiments  were  car  ;ied  out  with  open  steel  tube  holding 
liquid  pi  opellant, placed  upright  on  a pad  under  which  an  explosive 
charge  was  set  off.  Strain  gauges  in  the  liquid  recorded  the  progress 
of  the  shock  f row  bottom  to  top.  In  one  node37,  the  gauges  were  placed 
on  the  inner  tube  wall  and  measured  axial  and  hoop-strains,  presumably 
recording  peak  strains  upon  passage  of  the  shock  front.  In  the  second 
mode  , the  gauges  were  placed  on  fiberglas  partitions  (perpendicular 
to  the  axis)  blocking  partially  the  tube  inner  bore.  It  would  be  inter- 
esting to  assess  the  degree  to  which  normal  shock  propagation  and  axial 
flow  symnetry  were  affected  by  the  presence  of  these  thin  bulkheas; 
unfortunately  no  such  analysis  was  attempted. 

It  has  been  presumed  in  these  experiments  that  the  initiation  shock 
strains  the  tuba  walls  in  iuch  a manner  as  to  put  sufficient  regions  of 
the  liquid  in  tension  (or  sub-pressure) , causing  vapor  bubbles  to  form 
by  cavitation,  since  the  speed  of  the  shock  in  steel  is  higher  than  that 
in  the  liquid,  the  norma  shock  in  the  latter  would  then  pass  through 
cavitated  regions,  formin  hot  spots.  Successful  attempts  to  create 
LVD  were  reported39,  but  seems  that  the  experiment  was  not  designed 
to  verify  the  interaction  between  the  shock  and  the  cavitation  bubbles 
(or  the  number  and  size  of  these  bubbles)  even  qualitatively. 

A two-dimensional  computer  simulation  was  carried  out  to  determine 
conditions  for  LVD  in  a confined  liquid  propellant  charge,  simulating 
the  effect  of  hot  spots  by  varying  amounts  of  uniforr  heat  release  behind 
the  shock  front.37 
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K rn«  11—ntl — ml  analysis,  was  reported3*,  to  determine  charec- 
taristic  LVC  curvas  on  tha  p - i/p  plana,  again  traating  tha  reaction 
as  an  overall  hast  ralaasa  par ana tar.  This  formulation  involves  a 
homogeneous  energy  aquation  for  tha  region  behind  tha  shock  (where 
heat  ralaasa  due  to  decomposition  should  occur) , while  tha  Rankine- 
Hugoniot  relations  (jump  conditions)  contain  a term  denoting  heat  of 
reaction.  Thus,  whereas  the  modal  may  describe  certain  final  states 
of  low  velocity  detonation  in  mi.  approximate  manner,  it  certainly 
may  not  describe  the  process  by  which  this  form  of  detonation  evolves, 
and  these  is  no  assurance  that  the  final  statea  described  may  indeed 
by  achieved  physically. 

In  addition,  the  drgree  of  reaction  in  the  model  is  represented 
by  a generalised  "reaction  coordinate"  (given  as  s correlation  with 
pressure)  which  has  no  direct  relationship  with  the  chemical  process 
at  hand.  Further,  the  question  of  different  equations  of  state  for 
the  unreacted  liquid  propellant  and  for  the  mixture  of  decomposition 
products  and  unreacted  liquid  was  not  raised.  A single  (caloric) 
equation  of  state  was  considered,  relating  (e,p,o,f),  where  e is  thn 
specific  internal  energy  of  the  liquid  end  V the  reaction  coordinate. 
Following  tha  treatment  by  Nader 2 s and  the  comments  by  Enig  and 
Patrons2  , this  problem  seems  to  be  of  primary  importance  in  the 
codding  of  initiation  by  shock-hot  spot  interaction. 


G.3.  2 Collapse  of  a Vapor  Bubble 

The  analysis  of  collapse  of  a vapor  bubble  was  treated  in  an 
Appendix  in  both  reports37'38  and  was  carried  out  by  Levine  and 
Wooten  of  Ultrasystems . The  following  comments  can  be  made,  regarding 
the  more  recent  report : 3 8 

(1)  For  liquid  field  pressures  pQ  > 1000  atm,  (as  behind 
the  shock  front) , the  inertial  time  scale  associated  with  bubble 
collapse,  derived  by  Ra  yleigh1  is  Ro^pL//pn^>5  ^ 1 ysec  for 
R <f  1 mm  and  * 1 g/cm3.  Therefore  treatment  of  the  process 
of  collapse  is  relevant  within  the  time  scales  afforded  by 
shock  passage,  and  when  interaction  with  a shock  is  sought. 

However,  as  shown  in  the  works  of  Evans  et  al.23  and  later 
by  Mader 25 ' 2 8 ' 3 0 3 11  (both  analytically  and  in  experiments)  two 
important  effects  prevail  under  these  conditions.  First,  as 
an  initially  planar  shock  interacts  with  a spherical  bubble, 
the  shape  distorts  appreciably  during  collapse  and  becomes 
mushroomlike.  Second,  heat  transfer  due  to  conduction  would 
be  negligible,  as  the  thermal  diffusive  tire  scale 
R 2/a  >>  R_(p  /p  ) **  for  10  " ' R < 0.1  cm  and  hydrodynamic 
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convaction  is  expected  to  dominate.  Thar* for* , on  the  micro- 
ictlai  associated  with  bubble  collapse  (Rq  and  the  character- 
istic collapse  time)  fcwo-dimansional  effects  assume  relative 
importance.  Unfortunately,  these  effects  ware  excluded.  In- 
stead, spherical  syassetry  was  assumed,  the  bubble  center 
taken  stationary  relative  to  the  liquid  at  infinity,  and  tha 
(non -reacting)  heat  diffusion  equation  in  the  liquid  was 
appl ied . 

(2)  The  question  of  liquid  compressibility  within  this 
analysis  was  not  Mentioned , and  the  liquid  was  assumed  in- 
compressible as  indicated  by  the  use  of  the  bubble-equation, 
in  the  for*  of  5q.  6.1  . As  is  wall  known,  extremely  high 
pressures  can  evolve  in  the  liquid  during  bubble  collapse, 
under  which  all  known  liquids  exhibit  compressibility.  Thus, 
there  seems  to  be  no  physical  justification  to  assume  incom- 
pressible liquid  presently,  while  considering  compressibility 
for  shock  propagation  in  tha  LVD  analysis.  As  will  be  shown 
in  later  sections,  this  effect  has  a profound  influence  upon 
the  maximal  bubble  pressures  (and  temperatures)  achieved  in 
collapse . 

(3)  Reaction  is  assumed  to  occur  ir.  the  vapor  (not  in 
th*  liquid) , yet  the  assumption  of  equilibrium  vapor  pressure 
is  maintained,  and  a Claus iua-Clapeyron  expression  (with 
constant  heat  of  vaporization,  Q ) was  used  to  determine  the 
instantaneous  vapor  pressure  from  the  bubble  surf -ce  temper- 
ature. Clearly  this  '.isumption  is  not  expected  to  hold 

when  reaction  is  appreciable,  or  over  a wide  range  of  (p  ,Tq) 
over  which  Q is  changing  appreciably.  A more  appropriate 
form  for  the  mass  balance  between  liquid  and  vapor  (non-equi- 
librium) was  given  by  Plesset  and  Prosperetti. 4 *  6 * * * * * 

(4)  The  bubble  properties  were  assumed  uniform, although 

(regarding  the  temperature  distribution)  there  might  be  con- 
ditions under  which  this  assumption  breaks  down,  and  a thermal 
profile  must  be  considered,  as  in  the  work  by  Hickling.12 

The  enproy  equation  for  the  uniform  bubble  was  taken  from 
Hsieh  , and  appears  to  contain  a derivation  error.  The  first 

law  of  theiflvodynamics , for  the  bubble  volume, 

<5e  = 5q  - p6v 

The  variation  in  total  internal  energy  is 

oE  * 5(&e)  = oV5e  e5(cv) 

for  tne  case  of  ncnvanishing  mass  exchange  between  the  bubble 

and  the  liquid  cor.sidersd  c*-'  Hsieh 
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The  second  term  was  apparently  represented  by  0^6  (p  ) , but 


« - e°  + jfcv( T')dT'  * 

for  the  temperatures  considered,  with  e°  being  the  specific 
internal  energy  at  the  reference  state  (TO) , and  C the 
isochoric  specific  heat  of  the  gas. 


In  conclusion,  the  simulation  was  intended  to  represent  the  process 
of  bubble  collapse  in  a shocked  liquid.  The  various  deficiencies  pointed 
out  iu  the  proceeding  discussion  seen  to  indicate  that  the  actual  process 
watt  Inferred  only  in  an  approximate  way,  with  may  important  physical 
aspects  excluded.  In  this  respect,  the  analysis  of  the  same  process  by 
Mader28'30  35  seems  superior,  although  it  involves  a much  higher  level 
of  computational  effort. 


The  analysis  was  clearly  not  intended  for  modeling  of  secondary 
ignition  (due  to  hot  spot  formation  through  bubble  collapse)  during 
the  start-up  transient  of  actual  gun  charge  initiation,  in  the  manner 
considered  by  Knapton  et  al.3*  The  time  scale  for  this  process  is 
0(0.1  - 1)  msec  much  longer  than  the  shock-bubble  interaction  times, 
and  may  not  be  investigated  within  the  realm  of  a single  bubble-collapse. 

As  mentioned  earlier,  no  attempt  has  been  made37'32  to  combine 
the  analyses  of  LVD  propagation  and  bubble  collapse,  consequently, 
information  concerning  LVD  evolution,  failure  of  hot  spot  initiation 
or  LVD  stability  cannot  be  obtained  directly  from  the  model. 
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7.0  Physiol  Considerations  in  Theory  of  Liquid  Propellant  Compression  with 
isolated  Bubbles 

7. 1 chemical  Reactions  and  Vap  ization 

In  the  loading  of  liquid  propellant  (LP)  in  a gun  system,  bubbles 
might  be  preset t either  as  a result  of  unavoidable  ullage  in  the  supply 
system  or  due  to  cavitation  somewhere  in  the  flow  passages.  Such  bubbles 
may  contain,  therefore,  varying  amounts  of  air  and  LP  vapor.  The  presence 
of  air  under  most  practical  rapid  loading  conditions  must  be  expected. 
Vfcpor  bubbles,  even  without  air,  could  emanate  from  rapid  passage  of  the 
LP  over  sharp  edges  in  the  loading  orifice,  as  pointed  out  by  Batchelor17 
in  a discussion  of  various  cavitation  mechanisms  in  the  motion  of  liquids 
at  high  speeds.  Indeed  the  Flow  Visualisation  studies  of  Section  2.0  for 
dynamic  loading  of  the  liquid  monoprcpellant  display  this  behavior.  For 
these  reasons , realistic  modeling  of  hot  spot  initiation  due  to  collapse 
of  bubbles  should  incorporate  both  air  as  a permanent  gas  and  the  condens- 
ible vapor  of  the  LP  substance  within  the  bubble. 


Chemical  reactions  may  occur  both  within  the  gaseous  bubble  and  in 
the  high  temperature  liquid  layer  surrounding  the  bubble.  Depending  on 
the  particular  type  of  LP,  one  or  the  other  may  be  prominent.  If  the 
liquid  consists  mostly  of  volatile  propellant  material,  the  reaction  may 
proceed  within  the  bubble  as  well  as  in  the  liquid.  On  the  other  hand, 
if  the  liquid  consists  mostly  of  Inert  carrier  material  and  its  vapor  phase 
is  nonreactive,  the  likelihood  of  vapor  phase  reaction  diminishes,  and 
decomposition  can  be  expected  to  occur  mostly  in  the  high  temperature 
liquid  layer  surrounding  the  bubble.  To  present  a formulation  adequate 
for  both  LP  types,  reaction  terms  must  be  placed  in  both  the  liquid  and 
gas  phase  energy  conservation  equations.  Of  course,  this  greatly  com- 
plicates the  analytical  model  of  the  process,  since  in  addition  to  the 
"stiff*'  reaction  terms  incorporated,  it  becomes  also  necessary  to  follow 
chemical  species  concentrations,  unavoidable  in  any  realistic  treatment 
involving  chemical  change. 


The  process  of  vaporization  or  condensation  at  the  bubble  surface 
can  be  expressed  in  terms  of  the  Knudsen  relationship. 


ni  “ a[m  . - ra.  ] 
A cut  in 


(7.1) 


where  m is  the  net  mass  flux  of  the  vaporizing  species  A,  and  0 < a < 1 

A 

is  the  adsorption  coefficient.  The  outwa^ * mass  flux  (into  the  gas)  is 
given  oy  the  equilibrium  expression. 


B - f(p*q  ) 

out  vrA 


RTS* 

' 2iTW  ' ma 


(7.2) 


with  ’(T#)  being  the  equilibrium  vapor  density,  which  can  be  obtained 

from  a Clausius -Clapeyron  equation,  and  the  square  root  term  is  the  molecular 
velocity  of  species  k perpendicular  to  the  surface.  The  incoming  mass 
flux,  however,  depends  on  the  actual  concentration  cf  vapor  off  the  surface, 
end  hence  on  the  physical  and  chemical  processes  in  Che  bubble. 


f(p J ■ (t^t) 


where  (-)  denotes  values  at  the  gas  side  of  the  interface.  Actual 


(7.3) 


values  of  m.  and  m _ are 
in  out 


large  compared  with  due  to  the  large 


molecular  velocity.  Therefore,  at  quasi-steady  state,  the  assumption 
of  liquid-vapor  equilibrium  is  reasonable  for  calculation  of 

p ~ “ p*^(Ta) , although  the  liquid  may  undergo  vaporisation , and  m^  0. 

Under  dynamic  conditions,  however,  r relatively  small  change  in  p~  can 

A 

result  in  an  appreciable  change  in  m(t)  , which  car.  not  be  anticipated  when 
the  liquid- vapor  equilibrium  assumption  is  employed.  Thus,  :he 
combination  of  the  last  three  equations, 

R T h 

■>  _ _ V~  1 / t A \ 


A 


(7.4) 


where  p 


= p Y and  Y is  the  mass  fraction,  serves  as  an  adequate  gas- 

A A 


liquid  interface  condition  in  the  dynamic  process. 

In  view  of  the  large  variation  expected  in  bubble  temperatures  and 
pressures  daring  the  process  of  bubble  compression,  the  Claus ius-Clapeyron 
aq 

expression  for  should  be  used  in  its  integral  form  to  accommodate 

Qv  » Qy (Tg).  Ftor  instance,  note  that  Qv  i 0 in  the  neighborhood  oi  the 

critical  point  (P  . T ) 

B B 


CRIT. 


Characteristic 


Scales 


Heat  transfer  between  the  bubble  gas  and  the  surrounding  liquid  has 
an  extremely  important  role  in  proper  assessment  of  hot  spot  initiation. 
Therefore,  detailed  modeling  of  the  bubble  processes  can  not  be  avoided 
in  the  analysis,  since  the  thermal  war's  and  pressure  field  structure 
within  the  bubble  determine  the  gas  phase  side  of  the  thermal  balance  (the 
pressure  field  is  important  in  the  presence  of  gas  phase  reaction) . The 
question  is  whether  the  bubble  environment  can  be  ideally  considered  as 
uniform  at  all  times  (yet  time  dependant)  during  the  process,  without 
serious  loss  of  physical  significance  in  the  results.  Further,  within 
the  time  frame  of  start-up  pressurization  of  the  LP  charge,  the  question 
is  what  relationships  exist  between  the  rates  of  bubble  collapse,  thermal 
equilibration,  chemical  relaxation,  and  the  rate  of  pressurization  itself. 
To  obtain  some  preliminary  insight,  the  following  length  and  time  scale 
analysis  was  carried  cut. 


The  length  scale  associated  with  the  thermal  wave  within  the  bubble 
is  the  diffusive  characteristic  distance, 

*TH1  "*  ^aBTCOI^  ~ ‘(p„C  ^ Ro  (pi/^o^  ^ {7*S) 


whez » T ia  the  characteristic  collapse  time  taken  from  Rayleigh's 
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incostpressible  liquid  theory.  A second  associated  characteristic  thermal 
length  is  the  convective-diffusive  distance, 

y „ - a /q  - (-L-)(p  /P  )**  (7.6) 

XTH?.  TT  PjjCp  L ° 

where  U S Ro/TC0L  “ (p^Pj)*1  i*  the  characteristic  collapse  velocity. 

The  relationship  between  and  is 

XTH2/^™i“  XTHl/Ro  <7’7> 

These  scalar  can  be  calculated,  depending  on  the  driving  pressure  Pq, 

the  bubble  pressure  r \tio  p /p  and  the  reference  bubble  radius  R . 

so  o 

Utilising  the  properties  of  water  and  nitrogen,  the  following  values  evolve: 


(ata) 

► 

■V*o 

R (cm) 

O 

XTRl<Cn) 

XTH2/XTHl 

1 

10 

o.l 

0. 056 

0.56 

1000 

100 

10-2 

-4 

1.2X10"4 

-5 

0.012 

1000 

100 

10 

1.2x10 

0.12 

The  results  can  be  summarised  as  follows: 


(a)  At  low  pressures,  as  in  the  tarly  stages  of  coopression, 

(1)  relatively  large  bubbles,  with  R - 0(jL  net,  possess  highly 
nonuni  for*  thsrmal  fields,  since  ~ »0* 

(2}  under  these  conditions,  small  bubbles  with  R - 0 (1-10) um 
would  have  uniform  temperature  distributions. 

(b)  At  high  driving  pressures  and  pressure  ratios, 

(1)  a thin,  high-gradisnt , thermal  layer  near  the  bubble  surface 
is  inferred  tor  relatively  large  bubbles,  down  to  Rq  - 0(100))j, 

(2)  a nonunifora  thermal  fiald  is  indicated  again  for  small 
bubblas,  Rq  ~ 0(l-10)ya,  with  y^^  ~ R . 

The  uniformity  of  the  treasure  fiild  within  the  bubble  can  be  inferred 

fro*  a comparison  of  characti  ristic  time  seals s associated  with  fluid 

dynamic  equilibration,  t__  , and  bubble  collapse,  T : 

Ca  COL 


Vc 


R {p./Tp 

o L o 


JW 


(7.8) 


where  c denotes  the  adiabatic  speed  of  sound  in  the  gas.  For  the  properties 
of  water  and  nitrogen,  the  time  scale  ratio  can  ba  calculated: 


1000 


100 


0.02 


Evidently,  at  low  driving  pressures  and  low  compression  ratios  ths 

uniform  pressure  field  assumption  is  very  good,  while  at  higher  pressures 

it  tends  to  be  invalidated;  nevertheless,  the  interior  pressure  field 

can  be  assumed  uniform  to  a good  approximation  even  ac  the  high  values 

of  (P  j P_/P  ) shown. 

O BO 


Start-up  of  the  LPG  system  involves  an  onset  of  pressure  which 
applies  roughly  1000  atm  within  a period  of  0.1  to  1 msec.  During  this 
process,  gas  bubbles  imbedded  in  the  LP  charge  can  collapse  at  a rate 
which  may  differ  significantly  frasi  the  rate  of  pressure  start-up. 
Utilising  the  Rayleigh  time  scale, 


T 


COL 


rk(p/p) 

O L O 


*1 


the  calculation  yields: 


P (atm) 
o 

R (cm) 
o 

TC0L(mMc) 

1 

0.1 

0.1 

-2 

-4 

1000 

10 

3x10 

1000 

** 

1 

O 

H 

3xl0~6 

The  calculations  show  that  at  the  early  stages  of  tne  process,  for 
relatively  large  bubbles,  the  bubble  collapse  and  the  applied  external 
pressure  may  interact  dynamically,  so  that  the  time  variation  of  the 
imposed  PQ(t)  ®*y  not  be  neglected  in  the  oollapse  analysis.  However, 
for  sisall  bubbles,  with  Rq  ~ 0(1-10) jam  at  practically  all  stages  of 

external  pressurization,  the  rate  of  bubble  collapse  is  much  faster  than 
the  rate  of  PQ(t)  variation  and  PQ(t)  can  be  considered  as  quasi-steady 

during  any  collapse-rebound  cycle  period  In  this  instance,  the  bubble 
collapse  process  end  the  external  prasaure  variation  are  effectively 
decoupled  dynamically.  More  importantly,  if  the  viscous  dissipation  in 
the  liquid-gas  system  is  sufficiently  low,  the  proceeding  calculations 
indicate  that  (for  the  relatively  small-size  bubbles  considered) , practically 
hundreds  or  thousands  of  collapse-rebound  cycles  may  evolve  during  a single 
P0'Pressurization  period. 

The  possibility  of  rebound  of  the  bubble  following  full  compression , 
as  wall  as  the  periodicity  of  collapse-rebound,  were  mentioned  already  by 

Rayleigh 1 . To  assess  the  implications  of  this  process  upon  hot  spot 
initiation,  the  characteristic  collapse  time  should  be  compared  with  the 
bulX  thermal  relaxation  time  for  the  bubble  and  the  liquid,  defined 

TTOB  * Ro/aB  " Ro/(X/0bV 


(7.10) 
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and 


T 


THL 


RO/0,L 


(7.11) 


respectively. 

Por  the  properties 

of  nitrogen 

and  water, 

tthl 

tthb 

Isa  * 10-3 

tthl 

'o(atm) 

VPo 

R (an) 
o 

(msec) 

(msec) 

T X 10 
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7xl03 
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2 
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n 
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67 
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56 

70 

1000. 

100 

1 

O 

7xl0"3 

0.15 

47 

7 

The  results  show  chat 


(1)  Thermal  relaxation  within  the  bubble  has  a rate  much  slower  than 
the  collapse  rate. 

(2)  Similarly,  the  rate  of  liquid  thermal  relaxation  is  also  much 
slower  than  the  collapse  rate. 

The  last  two  observations  regarding  the  relationship  between  rates 
of  collapse,  start-up  pressurization  anu  thernml  relaxation  may  ba  r unsnarl zed 
as  follows: 


(1)  Bubble  collapse  and  rebound  occur  much  faster  than  characteristic 
start -up  pressurization, 

(2)  Thermal  relaxation  within  the  bubble  and  in  the  surrounding 
liquid  are  much  slower  than  typical  bubble  collapse. 

7.3  The  affects  of  Liquid  Compressibility 

The  relevance  and  importance  of  liquid  compressibility  to  the  present 
bubble  compression  problem  cannot  be  over  emphasized.  It  strongly 
influences  two  major  aspects  of  the  physical  process : the  dynamics  of 
pressure  waves  in  the  liquid,  and  the  balance  of  energy  between  the  liquid 
and  the  bubble. 


Neglecting  for  the  moment  the  (highly  localized)  conductive  term  in 
the  energy  balance  for  the  liquid,  and  considering  a non-reacting  liquid 
phase,  the  dynamic  process  can  be  expressed  in  terms  of  three  nonlinear, 
coupled  differential  equations  which  'orm  a hyperbolic  set.  As  pointed  out 

by  Lax  **  , such  nonlinear  hyperbolic  systems  can  obtain  weak  solutions, 
which  allow  the  formation  and  propagation  of  hydrodynamic  shocks,  observed 
also  in  practice.  Previous  treatments  of  bubble-liquid  dynamics  which 
account  for  liquid  compressibility  were  concerned  with  the  problem  of 

shock  propagation  in  liquids,  as  a scuased  in  the  book  by  Cole7. 

Prior  to  the  advent  of  high  speed  computers, many  of  the  analyses  sought 
closed  form  solutions  and  employed  suitable  approximations.  The  tost 
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well  known  Approximation  is  dua  to  Kirkwood  and  Batha  ' , baaed  on  tha 

assumption  that  5s  * 0 across  tha  liquid  shock  (nagligible  entropy 
change) , and  leads  to  an  expression  of  tha  total  (static  and  kinetic) 
enthalpy  distribution  following  an  outgoing  characteristics  line, 
denoted  by  I s 

r(h  + u? /2)  • const. 

This  result  was  utilized  in  a first  order  perturbation  theory 

by  Gilmore v , tnd  later  compare-*  gainst  numerical  solutions  of  the 
exact  equations  of  notion  (cont  * uitv  momentum,  assuming  isentropicity 

in  tha  liquid)  by  Hickling  and  Pl<  >set 1 0 , where  the  approximation  was 

demonstrated  to  be  in  good  agreeir  nt  with  the  numerical  results. 

Rsgarding  the  influence  on  hot  spot  initiation,  two  effects  may  be 
considered  in  this  instance.  First,  the  passage  of  hydrodynamic  shock 
(resulting  from  rebound  of  one  bubble)  over  a neighboring  bubble  undergoing 
compression  - or  at  early  stages  of  rebound  - can  significantly  contribute 
to  the  thermal  energy  accumulation  in  the  latter,  and  may  lead  to  hot  spot 
initiation  after  rebound.  Second,  wave  steepening  was  shown  to  occur10 
within  a distance  of  0(R  ) of  a rebounding  bubble  crater.  Since  com- 
pressed bubbles  are  generally  much  smaller  than  their  original  size, 

R(t)/R  « 1,  the  effect  of  shock  front  encountering  a small  disturbance 
(e.g. ,°another  bubble)  may  lead  to  enhancement  of  the  liquid  phase 
reaction,  a mechanism  for  detonation  investigated  by  Nader  on 

the  other  hand,  shocks  emanating  from  a rebounding  bubble  are  expected 
to  be  weak,  at  r/R  » 1 as  indicated  by  the  results  of  Hickling  and 
Plesaet10,  since  tfiey  are  driven  by  a small,  finite  amount  of  energy 
and  diminish  in  strength  due  to  spherical  dilation. 


The  second  physical  aspect  of  liquid  compressibility  can  be  invoked 
in  terms  of  the  energy  balance  in  a system  containing  a single  garious 
bubble.  For  simplicity,  only  the  fully  collapsed  state  (minimal  bubble 
radius)  is  considered,  but  the  results  apply  to  all  states  of  bubble 
surface  motion.  If  the  liquid  field  pressure,  P0  - Const.,  the  driving 

ttroke  work"  for  the  collapsing  bubble  is  given  in  the  incompressible 

liquid  case  , by 


- P 5V_  - -P  -7  (R3  - R3) 

STROKE  o B o J min  o 


(7.12) 


However, 


when  the  liquid  is  compressible,  the  above  expression  is  modified,  as 

00  p (r) 

"sitoo:  ■ po4vb  • Po  J04w,dr  j d<1/0>  |7-131 
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whare  the  second  term  cot  responds  to  P Av  k the  pert  of  the 

O Xs 

stroke  work  absorbed  by  liquid  compression.  At  the  fully 
compressed  state,  the  stroke  work  is  converted  entirely  into 
compression  work  (considering  viscous  and  other  losses  negligible) 
so  that  for  the  incompressible  liquid t 


STROKS 


Xln  fc> 

■ X * V. 


4irPQ(R)  R^dR 


(7.14) 


where  the  functional  dependence  of  p (R)  is  determined  by  the  nature 
of  the  bubble  compression  process  B (isothermal,  adiabatic,  or  inter- 
ned i ate)  . On  the  other  hand,  for  the  compressible  liquid  case, 


STROKE 


P d 

B 


7 °/r) 

-Ip  4irr2dr  I pd(l/p) 


(7. IS) 


where  the  additional  term  arises  from  liquid  compression  work.  Now 

Fines  p(r)/P  > 1 for  R < r < <*, 
o — 

7 (Ur)  &it) 

pr2dzj  pd(l/p)>  - P pradr J d(L4>)  (7.16) 


y Jir) 

prJdr j d?l/p) 


- r uvv 
o L 


where  Pa  is  a constant  mean  pressure  in  the  liquid,  PQ  < P^  < p(R). 

Therefore  the  liquid  compression  work  always  exceeds  the  associated 

driving  stroke  work  P Av  . This  means  that  when  the  liquid  is 

O L 

compressible  (as  all  liquids  are) , less  of  the  stroke  work  is  available 
for  bubble  compression.  In  particular,  this  conclusion  leads  to 
the  following  comparisons t 


COMPRESSIBLE 


INCOMPRESSIBLE 

LIQUID 


Indeed,  the  calculated  results  of  Appendix  A for  the  incompressible 
case  show  enormous  maximal  bubble  pressures  »und  temperatures)  compared 
with  the  approximate  solutions  for  the  compressible  liquid  case,  given 
in  Appendix  B ; the  above  results  tend  to  be  strongly  verified. 

o 

Comparison  at  several  initial  bubble  pressure  ratios  PB/p0  end 
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d.;ivlng  pressures  P is  given  in  ths  following  table  along  with  the 
fraction  of  coo^rosSion  work  absorbed  in  ths  liquid  to  its  compressibility » 
;.ll  the  tabulatsd  results  are  for  ths  adiabatic  bubble  approximation. 


TAR  LX  4 


COMPARISON  Or  COMPRESSIBLE  AMD  INCOMPRESSIBLE 
LIQUID  SOLUTIONS  AT  MAXIMAL  COMPRESSION  STATE 


o 

PB 

P (atm) 

o 

*S«in 

PBaax 
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. *CL 

P 

o 

R 

o 

P “ “1 
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2.4 

(a) 

:>oo 

.277 

21 

4.67 

13 

( 

1000 

.298 

15 

4.3 

29 

NOTEi (a) 

compressible  liquid 

approximate 

solution  is  in  Appendix 

7.4  SumsMtiry 


Tbs  following  preliminary  conclusions  may  be  draw.,  based  on 
ths  physical  considerations  of  ths  precssdir.g  sections. 

(1)  Liquid  compressibility  has  a profound  effect  upon  thw  actual 
maximal  pressures  and  temperatures  obtain*!  in  a collapsing  gaseous 
bubble.  Corns ».red  with  results  for  incompressible  liquid,  final  bubble 
presaurss  are  roducsd  by  two  orders  of  magnituda  or  nx>r«  and  tsmp% raturss 
by  at.  least  one  order  of  magnitude,  when  the  bubble  is  considered 

ideally  adiabatic  (fast  compression) . This  imp~i.es  that  liquid  cooipressibility 
must  be  taken  into  account  in  the  analysis  of  hot  spot  initiation. 

(2)  The  applied  start-up  pressurization  tine  ia  typically  much 
longer  than  the  inertial  bubble  collapse  tiw*  c,_.  ■ R (P.P  )**  for 

the  range  of  bubble  radii  and  pressures  considared,  and  almost  always 
riuch  ihcrter  than  ths  bulk  thermal  relaxation  time  of  the  bubble-liquid 
•yetis*.  Therefore,  for  sufficiently  low  viscous  damping  in  tha  liquid, 
many  collapse-rebound  cycles  can  be  expected  during  start-up  pressurization, 
over  which  heat  will  accumulate  in  the  bubble  and  in  the  immediately 
surrounding  liquid.  Thus,  whereas  a single  bubble  collapse  may  be 
too  brief  to  lead  directly  to  hot  spit  initiation,  s aerie*  of  hundreds 
or  thousands  of  collapae-rebound  cycles,  during  which  ths  external 
driving  pressure  is  continuously  rising,  may  under  favorable  interaction 
with  tha  chemical  kinatica  lead  finally  to  runaway  reaction. 


? 
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8 * 0 Analytical  Modal 

The  analytical  modal  concarni  a single  gaa  bubbla  in  an  infinita 
liquid  medium.  Conaaquantly , the  treatment  ia  dividad  into  thraa  parts: 
tna  liquid  phase,  tha  gaa-bubble  intarior  and  tha  interface  conservation 
conditiona  at  tha  bubbla  surface.  Oaapita  *-ha  indicat iona  that  apharical 

bubbla  ahapaa  are  sometime#  unstable,  aa  discussed  by  Pleaset*  and 

Plaaeat  and  Prospereiti  , * tna  assumption  of  apharical  bubbla  ahapa 
ia  mada,  to  obtain  a mathamatically  tractabla  nodal  at  this  stage. 

8.1  Tha  Liquid  Phaaa 

Assumptions  ragarding  tha  liquid  phaaa  ara  aa  follows: 

a)  inviacid,  homogeneous  and  iaotropic  liquid 

b)  body  forces  ara  naqllqibla 

c)  irrotational  liquid  motiona 

d)  no  aourcaa/ainka  of  rnaaa,  momentum  and  anargy 
a)  apharical  symmetry 

f)  conpraaaibla  liquid  phaaa 
a)  surface  tanaion  ia  naqliqibla 

In  spherical  polar  coordinates,  tha  aquation  of  motion  ara  written 
in  Euler ian  form: 


3p/3t  + div  (cu)  - 0 


(8.1) 


3(pu)/3t  ♦ div(puu)  « - div  (P) 


(8.  V) 


3 (pa  )/3t  + div(pue,„  ♦ uP)  - div(Agrad  T)  ▼ q 

T A R#  L 


(8.3) 


u « ue^  ia  tha  only  component  of  velocity  conaidared  (radial) ; tha  total 
apacific  internal  energy  ia  - e(T)  ♦ uz/2  and  tha  pressure  tanaor  ia 
q ia  tha  la at  release  by  chemical  reaction  in  tha  liquid.  Tha  above 
formulation  follows  those  derived  by  Williams  k0  or  Batchelor  17 

After  taking  the  assumptions  into  consideration , tha  aquations  of 
motion  for  continuity,  momentum  and  energy  can  be  written,  in  simplier  form: 


ft 


1 9,2 

+ 7*  Ji (r  Du) 


(8.4) 


9u 


3u 


9t  * U3r 


u 1 3p 

~ ♦ ~ T*  - 0 

r o 3r 


! 


! 


| 


i 
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3e 

5t 


dr 


* ut — + 


p 3u 


o r 


3/-  1 ill)  + 

?7<r 


Vi/° 


(8.6) 


for  t > 0 and  V < r < ®. 

The  consideration  of  chemical  reaction  requirar  the  addition  of  species 
conservation  equations;  for  the  liquid  phase  mass  fractions  Y , assuming 

all  binary  dif fusivitiea  to  be  equal  to  the  mean  pD: 


av  3Y  , 3 9Y. 

5t  + U5T  ’ SF  57<r,‘LD5T)  + aLj  Vp 

j ■ 1 , 2 * * • N 


(0.7) 


where  a.  5 t(.(yj  - v')  , v'  and  v"  being  the  stoichiometric  coefficients 
Lj  3 3 L 3 3 

of  species  j in  the  liquid  reaction,  and  N is  the  total  number  of  such 
species  participating,  for  the  case  of  a single  overall  reaction  step, 
Shvab-Zeldovich  coupling  terms  stay  be  formed;  defining 

S.  = Y./d,  j-1,2,3  • • -N  (0.8a) 

j 1 S 

such  that 

♦j  3 8l  - 8j  , J-2,3,*  • * N (8-8b, 

and  the  differential  operator  for  these  0 terms  is  free  of  the  reaction 
term,  vis. , 3 


JT  * V*5T‘  * pT*  3?(rJpD§TL)  ' 3"2'3*  ‘ 'N  <8-9) 

Consider  an  overall  one-step  decomposition  reaction  of  a liquid 
phase  species  A,  forming  products  which  remain  dissolved  in  the  liquid, 

A ( iq)  — ►B1(g)  > B^g)  ♦ ■ * * ♦ 

with  Q being  the  heet  of  reection.  If  we  represent  approximately  all  the 

u 

products  B^(g)  as  a single  mean  species  (inert),  the  reaction  scheme 
simplifies  to 


A(liq)  — ► Big)  ♦ Ql 


(8.10) 


If  the  liquid  consists  of  reactant  A and  an  inart  diluent  C,  its  masr 
fraction  reuins  invariant  at  all  time,  Yc<r,t)  - const.,  and  only  two 

species  equations  have  to  be  solved,  in  the  form  aiven  bv  Eqs.  (C.7  ) 
end  ( 0.9  ) . 


f 


-101- 


The  liquid  co ipressibi lity  9tat-*!T'*nt  (i.e.,  equation  of  state)  is, 

CO0  - P ♦ B (0.1 

which  has  bean  usad  by  Hsieh1  * as  wail  as  by  Hilling  and  Plasset10. 

Tor  watar,  tha  values  of  n ■ 7 and  B » 3000  atm  * 're  proposed.  This 
corresponds  approximately  to  a bulk  modulus  of  30 3,000  psi  at  one 


atmosphere  pressure.  More  general  expresoions  in  /olve  polynomials 
including  powers  lower  than  n • 7.  The  above  compressibility  equation 
is  for  an  isothermal  process,  and  hence  not  strictly  applicable  to  the 
present  general  formulation;  a more  adequate  expr>  ision  can  be  derived 


from  basic  thermodynamic  principles,  cf.  Moalwyn-*.  ighes1' 1 . 
energy  statement  is 

e (T)  - Cv T 

and  the  reaction  ten.'  is 


q 


R,  L 


Vlit> 


The  internal 

(8.12) 

(8.13) 


where  W (T)  is  the  reaction  »ta,  which  can  be  fitted  to  an  Arrhenius 
expression . 


It  should  be  noted  that  thermal  transport  ia  expected  to  bn 
negligible  over  most  of  the  liquid  region,  except  within  a thin,  high 
gradient  layer  adjacent  ot  the  bubble  surface,  where  temperature  is  also 
high  due  to  transfer  from  the  bubble.  It  would  be  also  safe  to  assume 
that  the  reaction  term,  q , is  relatively  important  within  the  same 

layer  (provided  tha  activation  energy  is  suf f icien>.ly  high)  anr  negligible 
outside.  Ttms,  for  tha  liquid  domain  outside  of  this  thin  thermal  layer, 
the  assumption  of  isentropic  liquid  may  be  employee,  in  line  *ith  the 

analysis  of  Gilmore  * and  Hickling  and  Plesset10  . A suitable  length 
scale  associated  with  the  thermal  layer  thickness  is 

XTHL  “ (8.14) 

which  arises  naturally  from  the  energy  equation,  and  where  the  characteristic 
inertial  collapse  time  was  defined  in  Eg.  (7.9).  Explicitly, 


utilizing  the  properties  of  water, 


P (atm) 

o 

R (cm) 
o 

XT. 

i 

0.1 

4x10  4 

1000 

10~2 

2xlO~5 

1000 

1 

O 

H 

2x10" 6 

it  is  demonstrated  that  y 

THL/Ro  « l- 
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The  fo ) lowing  boundedness-conditions  are  relevant : 

u (°°,t)  m 0 (8.16*.) 

o(°°,t)  • P (t)  (8.16b) 

o 

where  P (t)  is  the  imposed  far  field  pressure,  which  may  vary  with  time 
o 

in  a specified  manner (*  j- , ramp  function);  however,  the  initial  pressure 
distribution  in  the  field  is  "v-.ifortn: 

p(r,o)  - P (0)  (8.17) 

o 

Cle surly,  tna  liquid  field  problem,  as  posed  by  the  proceeding 

formulation  and  boundary  conditions  is  quite  formidable,  and  considerable 

simplificavior  may  be  achieved  if  P (t)  • oonst.  is  specified.  This  type 

o 

of  simplification  has  been  used  in  cavitation  thaory. 

8.2  The  Gas  Phase 


The  following  assumptions  are  made,  regarding  the  bubble  interior: 

a)  the  bubble  environment  consists  of  a mixture  of  permanent  gases, 
as  well  as  vapor  from  the  liquid 

b)  the  gas  mixture  may  undergo  chemical  reaction,  represented  by 
a sxngle  overall  reaction  step 

c)  assumptions  (a)  through  (f) , made  for  the  liquid  phase,  are 
applied  to  the  gas  phase  as  well. 


Consequently,  the  gas  phase  system  of  equations  has  exactly  the  same 
form  as  the  liquid  phase  system,  employing  Eulerian  form. 


An  additional  simplifying  assumption  can  be  introduced,  regarding 
the  pressure  field  within  the  bubble.  During  the  important  final  stage 
of  the  process  of  collapse,  the  pressure  within  the  bubble  is  expected  to 
be  much  larger  than  the  liquid  pressure  at  infinity,  while  the  liquid 
density  would  always  be  higher  than  that  of  the  gas.  Therefore,  it 
can  be  expected  that  the  gas-dynamic  equilibration  txme  within  the  bubble 
would  be  always  nuch  smaller  than  the  characteristic  process  (bubble 
collapse)  time,  vis. , 


R/c  « R (p_/P  )**  (8.18) 

g o L o 

where  c =*  (yp  /p  J*5  is  the  adiabat:'c  speed  of  sound  in  the  ges  and  y the 
cj  go  ‘ 

specific  heat  ra.j.o.  This  affords  the  following  simplification. 


p(r . t)  = p (t)  , 0 < r < R 
“ 


(8.19) 


namely,  the  pressure  is  taken  uniform  within  the  buhllo  at  all  trmes. 


The  system  of  equations  describing  the  process  in  the  bubble 
interior  is  therefore, 


II-  * p - ° 

(8.20) 

3u  . 3u  n 

5T*"5F  -° 

(8.21) 

3a  , 3t  P 3a  1 9,  i.9T.  , qR 

5t  + + o 3?  " pF  3?<r  *5?  r 

(8.22) 

for  t > 0 and  0 < r < R. 

lb  accommodate  the  elevated  pressures  expected  during  the  period  of 
collapse,  the  Noble-Abel  equation  of  state  is  utilized, 

p<£-  nB>-  RUT/W  (8.23) 

ii  2 ■ 

where  nB  is  the  covolume  # typically  1 cm/g , and  ¥ the  Man 
molecular  weight  of  the  gas. 


The  reaction  in  the  gas  phase  may  be  first  order  or  second  order 
overall , 


A(g) 

— ► B(g)  + Qg 

(8.24a) 

a,  (g) 

Jk, 

+ A^(g)  — B(g)  +0 

* g 

(8.24b) 

where  B(g)  in  both  cases  denotes  the  Man,  inert  product  mixture.  The 
presence  of  a diluent  C(g)  is  also  anticipated,  again  with  invariant  mss 
fraction.  The  reaction  rates  for  the  chemical  Mchaniama  shown  are, 
respectively: 

oj(i) 

g 

- pA^xpt-Ij^/T)  (y^/wA) 

(8.25a) 

N 

w (p 

3 

- oiA2*xp(-i2/r)  (ya  /¥x)  (va  A2) 

(8.25b) 

The  species  conservation  equations  are  of  the  same  form  as  those  in  the 
liquid  (of  course,  using  gas  phase  properties)  and  shall  not  be  repeated. 
The  Man  molecular  Might  in  the  equation  of  state  is  defined 
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Zt  should  bs  enphaeiaed  that  dus  to  ths  presence  of  permanent  gas 
in  ths  bubble , ths  total  interior  pressure  is  not  equ.1  to  ths  vapor 
pressure,  and  is  not  .lets rained  by  the  phase  interface  temperature. 

The  associated  boundary  conditions  are  given « 

u(0,t)  - O 

3r/3r(0,t)  - 0 


8.3  The  phase  Interface  pan nervation  Conditions 

Following  Plesaet  and  Prospers tti  * , the  pressure  force  balance 

at  the  bubble  surface  is 


- P»  - Tf  - ^ u (8-27> 

where  is  the  interior  bubble  pressure  ( ar f umed  uniform  in  the  gas 

phase  treatment) , and  9 and  y are  tho  surface  tension  constant  and  the 
liquid  viscosity  respectively.  For  the  range  of  pressures  and  bubble 
si. see  considered,  p (R+)  - o( 1000) atm,  R - 0(2)  ym  utilising  the  properties 
of  water  ( a • 150  dyne/cm  and  y * 0.01  g/cm-s)  the  two  last  tens  can 
be  shown  to  be  negligible  coopered  with  the  pressure  p(R+) . Hence, 


- PB(t)  (8.28) 


- T(R~,t)  - Ts(t)  (8.29) 

i.e.,  teaperature  is  continuous  across  the  interface. 

Heat  and  mss  exchange  across  the  phase  interface  uay  occur  during 
the  cavitation  process  through  conduction  and  vaporisation.  Considering 
vaporisation,  the  bubble  surface  velocity 

U(t)  = dR/dt  - 0 + v (8.30) 

e B 

where  0 is  due  to  coaqaression  or  expansion  of  ths  bubble  gas,  and  v (t) 
i*  the  surface  regression  rate  due  to  vaporization,  defined 

vg  - a(t)/PL(R+,t)  (8.31) 

with  ■ the  vaporisation  mass  flux.  Consequently,  ths  velocity  on  the  gas 
side  of  the  interface  can  be  expressed  as 

u (R“,t)  - U - m/p  (R~,t)  (8.32) 

9 


P(R+,t) 
In  addition , 

T(R+,t) 


(8.26a) 

(8.26b) 
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taring  tha  rapid  oollapaa  pcooaaa  oonsldarad,  the  a priori 

asswption  of  equilibritsi  vaporisation  aa  nada  in  tha  SRI  mart"  ia 

not  warranted.  Sanaa,  tha  mogm  general  nonequilibrium  expression  ia  propoaad, 

% • '".-"a-1*’  <»r  - v 

as  given  by  Xqa.  (7.1)  through  (7.4). 

Furthor,  for  tha  case  of  water,  tha  latent  heat  of  v^artsatiou 
variee  wivh  taaparatura,  and  vanishes  at  tha  critical  point  Therefore, 
a general  interface  energy  constraint  would  involve  ■ Q^CT  ) . following 

Scale  and  Sutton  '*  and  WUliaaa  M the  conservation  conditions  for 
enthalpy  and  species  ness  fractions  are  given  respectively, 

ah(K~)  - (Vcp)9  3h/3r(R~)  - (3.33) 

ah(R+)  - U/C 3h/3r(R+)  - »*IQV(Y.) 


«d»re  Yl  denotes  the  evaporating  liquid  aass  fraction  (on  tha  liquid  side, 
a*-  R*)  and  Qv(Ta)  is  the  latenc  heat  of  vaporisation. 


mY  (r“)  - ((.D)  3f./3r  (R~)  - 
J 9 3 

■Yi(R+)  - (0D)l  3Y1/3r(R+) 

j-1,2-  • -N 


(8.34) 


In  conclusion,  it  should  be  enphasised  that  the  foraulatian  given 
in  the  lest  two  sections  is  strictly  for  the  spherically  sywnetric 
case  where  the  field  pressure  at  infinity,  P (t) , stay  vary  in  a contin- 
uous wanner,  and  is  intended  to  simulate  ths°rela lively  slow  start-up 
transient.  As  pointed  out  earlier,  this  wight  involve  many  collapse- 
rebound  cycles,  during  which  h-st  is  accvanilstsd  in  a region  roughly  the 
six*  of  the  initial  bubble.  so~* 5 
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9.0  Method  of  Solution 


9.1  flw  Two  Time  Seal*  Problem 

Within  the  realm  of  an  investigation  of  propellant  initiation  during 
the  start-up  transient  by  hot  spots,  through  gas-sapor  bubble  compression , 
the  foregoing  physical  considerations  have  pointed  out  the  poeaible  exis- 
tence of  two  distinct  tins  scales  (at  least)  for  the  problse.  Disregarding 
for  the  ac—nt  chemical  kinetic  ties  scales,  the  two  distinct  characteristic 
tisws  are  i 

(1)  lbs  typical  start-up  transient  time,  Tq - 0(0.1 -1)  msec 
characterises  the  variation  of  the  field  pressure , pQ(t). 

(2)  The  characteristic  bubble  collapse  time,  T ■ T defined  earlier 
0(1)  msec. 

Evidently,  rQ  » T.  for  the  problem  considered.  It  should  be  stressed 
that  an  oscillatory  pressure  field,  with  p (t)  having  both  isonotonically , 
varying  as  well  as  periodic  components  is  not  excluded,  but  the  condition 
an  the  wavelength  of  these  oscillations  is  that  they  be  much  logger  than  k_, 
the  initial  bubble  radius,  and  their  period  be  much  larger  than  T___ 
[otherwise,  if  the  opposite  extreme  is  considered,  the  short  wavelength , 
shock-bubble  interaction  problem  addressed  by  Nader2* ,ss  **  is  recovered  .] 

Thus,  with  viscous  dissipation  sufficiently  low  and  with  bubble- 
stability  permitting,  the  overall  problem  concerns  many  collapse -rebound 
bubble  cycles  ever  the  period  of  a single  start-up  transient.  Over  these 
periods,  pressure  and  temperature  within  the  bubble  would  oscillate  with  a 
relatively  high  amplitude , but  their  mean  values  are  expected  to  vary  as 
well,  on  a longer  time  scale.  This  is  based  on  the  observation  that 
thermal  dissipation  in  the  liquid  surrounding  the  bubble  occurs  on  a time 
scale  much  longer  than  x , and  therefore  heat  nay  accumulate  in  the 
bubble  over  several  periods  during  which  the  driving  field  pressure  pQ (t) 
somewhat  increase*. 

When  chemical  reactions  (in  the  bubble  and/or  within  the  high-temper- 
ature propellant  region  surrounding  it)  are  considered,  the  associated 
kinetic  time  scale*  have  crucial  importance  regarding  a process  that  may 
or  may  not  lead  to  secondary  ignition.  Of  course , these  chemical  time  scales 
cannot  be  considered  outside  of  the  fraamwork  of  the  process  at  hand,  as 
they  are  strongly  temperature  (end  pressure)  dependent. 

Depending  on  the  magnitude  of  the  kinetic  time  scales,  the  processes  of 
bubble  collapse  end  chemical  heat  release  may  interact  over  a short  time 
scale  (within  a period  of  collapse-rebound)  or  a long  time  scale  (within 
many  periods  of  collapse-rebound) . Even  with  a given  set  of  kinetic  para- 
meters and  chemical  mechanism,  it  would  be  difficult  to  assess  a priori  what 
type  of  interaction  should  be  expected,  owing  to  the  highly  nonlinear 
coupling  between  the  chemical  and  bubble -collapse  processes.  Thus,  analysis 
is  proposed  along  two  major  lines,  as  follows. 
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(1)  li»«»tlytt loo  of  the  short  tiam  seal*  prooess.  This 
analysis  is  si  nail  at  sol  Tin?  tit*-  bubble-propellant  inter- 
action  protolan  over  a nasi  1 number  of  aollap— -rubound 
parioda,  daring  which  the  driving  field  pressure  variation 
is  naall  or  vanishing. 


(2)  invastigation  of  tha  long  tins  seals  behavior,  or  the 
naan  prooeas.  this  part  of  tha  analysis  requires  re-easting 
of  the  problem  in  ta^s  of  two  ties  scales  (t  ,t  ) , and 
involves  transforming  the  aquations  given  in  ths^preoeeding 
section.  Them,  following  qaasi-1  inear  last  ion,  a mathematical 
averaging  method  would  be  appli  d,  ao  that  tha  anall-seale 
oscillations  are  considered  only  in  terms  of  their  mean 


ivalopa , which  nay  vary  with  the  long  tine  scale,  T 


Although  put  in  two  distinct  categories,  both  analyses  should  be 
carried  out,  in  tha  order  indicated.  The  solutions  to  tha  nicoro-scele 
process  are  required  to  obtain  a firm  guidance  for  the  linearisation  in- 
volved in  the  two-time  scale  analysis,  over  a wide  enough  range  of  initial 
conditions,  thermophysical  and  kinetic  parameters. 


*•2  Outline  of  tha  Solution  Procedure 

9.2.1  short  Tima  Scale  Analysis 

The  system  of  equations  for  tha  liquid,  gas  bubble  and  phase 
interface  are  utilised  as  given  in  Section  8.0.  A single  time  scale  and 
length  scale  are  proposed,  to  obtain  dimensionless  coordinate a: 

< 5 "V  " 3 t/'Toca.' 

The  sets  of  oonfarvation  equations  i-or  the  liquid  and  for  tha  gas  are  than 
solved  separate  y by  the  general -implicit  predictor-corrector,  finite 
difference  algo*ithm,  applied  to  the  • yatam  in  Culcriar.  form.  At  each  step 
during  tha  time  .tarching  procedure,  taa  liquid  and  gas  solutions  are 
matched,  to  satisfy  tha  bubble  interface  conservation  condition#.  A similar 
procedure  waa  successfully  applied  to  tha  problem  of  nitrsmine  deflagration 
by  Ban  Bauven  et  al.,"'  5 solving  simultaneously  for  tan  independent 
unknown  variables ; time-marching  was  need  for  obtaining  a final  stoady 
stats  solution  from  a sat  of  initial  data. 

9.2.2  Two  Tima  Scale  Analysis 

This  analysis  la  intended  for  elucidation  of  tha  long  time  scale 
behavior  of  tha  bubbla-propallant  systmi.  in  order  to  carry  out  tha  analysis, 
properties  which  oscillate  (a.g. , pres  ure  and  temperature  in  tha  liquid  and 
tha  gas)  are  designated  two  seta  of  independent  time  coordinates,  corresponding 
to  Tq  and  T , such  that  each  variable  has  an  amplitude  and  phase  which  are 
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slowly  varying  (~tq)  relative  to  tha  rapid  oscillatory  action  of  collapse- 
rabound . The  objective  ia  to  aaparata  in  this  aannar  tha  rapid  (short- 
tiaa  a cals)  behavior  froa  tho  alow  behavior  and  investigate  only  the  latter, 
naawly , tha  evolution  of  tha  af or  Mentioned  aaplitude  and  phase  in  r and 


Pro  blase  of  this  nature,  involving  both  ranid  and  aloe  tins 
variation,  arise  aoaetiass  in  wave  propagation  (dispersive)  through  slowly 
varying  aadia.  A xuique  Lap  angian-app  roach  for  waves  in  adiabatic  systaae 
was  developed  by  Nhlthan, * * "®  using  an  averaging  aethod.  Luke"*  investi- 
gated a similar  problsa  utilising  a singular  perturbation  technique,  with 
the  soall  perturbation  quantity  being  the  ratio  of  tiaa  scales  t./tq. 
Knakal1'  suggested  e nonlinear  MKB  approach  for  tha  classical  Kortiweg-de 
vriaa  (shallow  water  wave)  equation. 

Tha  scope  of  the  proposed  analytical  effort  is  outlined  as  follows: 

(1)  Oarivatlon  of  a hierarchy  of  partial  differential  equations  in 
tanas  cf  tha  alowly  varying  properties  ( amplitude  and  phase) 
of  the  impendent  variables,  utilising  a perturbation  aethod, 
the  snail  perturbation  quantity  being  the  ratio  T /xQ  « 1. 

(2)  Nuswrical  solution  by  finite  difference  to  t(ia  seroth  end 
first  order  problems . 

This  soda  of  analysis  i»i  necessary  for  understrnd.'.r>g  of  hot  spot 
fvtrsation  during  the  start-up  transient  of  LPG  systems,  which  uay  lead  to 
secondary  ignition  within  tha  propellant. 

although  it  is  possible  in  principle  to  extend  the  numerical  short 
time  seals  analysis,  proposed  in  tha  preoseding  section,  over  aany  periods 
of  bubble  collapse-rebound , this  is  highly  undesirable . it  would  involve 
a vary  large  cost  per  single  caaiputer  solution  and  the  errors,  accumulated 
over  aany  time  steps  of  integration  would  greatly  reduce  tha  reliability 
of  the  final  results. 
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APPVNDXX  A 


Pressure  Distributions  During  Rubble  Collapse  for 


The  purpoM  of  the  present  calculations  it  to  obtain  tha 
instantaneous  pressure  fiald  about  a gaseous  bubble , in  tha  prooaaa  of 
collapaa  in  an  infinite,  incompressible  liquid.  Tha  original  spherical 

bubble  ooilapaa  probier.,  as  poaod  by  Rayleigh 1 la  addxesaed.  *ralytical 
solutions  for  tha  bubbla  wall  valocity , U S dR/dt,  and  tha  praaauue  fiald 
in  tha  liquid,  p(r  ,t),  ware  obtainad  by  Rayleigh,  for  tha  oaaa  of 
isothermal  bubbla  procaaa.  Tha  bubbla  wall  valocity  ia  given  by 


(l/n1  - i>  - — (iai*) log  (i/n), 
PL 


(A-l) 


p V - Const. 

B B 

wham  P « Const,  is  tha  initial  fiald  pressure,  ia  tha  initial  bubbla 

O B 


pressure,  p is  tha  liquid  danaity  and 

n = R/R 


tha  dijnansionlass  instantaneous  bubbla  radius. 


for  tha  adiabatic  bubbla  procaaa,  a straightforward  calculation  yields 


2P 

o'  -TT* 


« si..  tiW**1’  i 


PL 


\A-2) 


pbvbY  " 0on,t> 

following  tha  analysis  by  Raylaigh,  whan  y -c^/c^  for  tha  gas.  In  tha 

derivation  of  Eg.  (A-2)  , tha  same  idealisations  used  by  Raylaigh  have  bean 
employed,  in  particular,  that  no  mass  exchange  occurs  between  tha  bubbla 
and  the  liquid,  and  that  tha  bubbla  porperties  arc  spatially  uniform . 

Tha  instantaneous  liquid  fiald  pressure  distribution  (for  incompressible 

liquid)  corresponding  to  tha  Raylaigh  problem  was  given  by  Plasaat  and  Proaparetti ' 
in  tha  form  of  tha  general  Bernoulli  equation : 


p ■ Po  . ,*»  - Po 


L,  ..a, l - 


(A-3) 


whara  Q 3 r/R(t) 


* * 
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and  p ( c , n ) is  t.—  liquid  pressure,  while  p (n ) is  the  bubble  pressure 

0 

given  by 


ISOTHERMAL ■ A = <p.  - P ) /P  - <p°/P  )/n3-l  <A_4> 

B BOO  Bo 

AOIABATICi  ♦_  1 (p_  - P )/P  - (p°/P  )/n3Y-l  (A“5) 

B BOO  BO 


The  state  et  full  coeipreesion  is  defined  by  U • 0.  Thus,  Kqe.  (A-l) 
end  (A-2)  cen  be  used  to  define  *or  isothermal  end 

ediebetic  eases  respectively  (in  preetice,  due  to  the  algebraic 
complexity  of  the  u -equations , the  inverse  calculation  is  Much  easier 
to  perform  nasMly,  finding  pg/PQ  as  a function  of  n ^).  Than  the 

staxiaal  bubble  pressure  ratio  s-g  can  be  calculated  by  employing  the 

values  of  , ;p°/P  ) in  Bqs.  (A-4)  and  (A-S).  The  results  of  these 
■in  b o 

calculations  are  given  in  Figures  A-l  and  A-2.  Indeed  it  is  shown  that 
the  isothemil  case  leads  to  more  violent  collapse,  nasMly  to  significantly 
higher  dg  and  lower  than  the  adiabatic  case,  as  expectsd. 


The  pressure  distribution  in  the  liquid  >(C;H)  during  bubble  collapse 
can  now  be  calculated,  by  substitution  of  the  velocity  U2  from  Eqs.  (A-l) 
and  (a-2)  into  Eq.  (A-3) , for  specified  <_  n < 1»  <*»*••  Thi* 

will  yield  distributions  due  to  isothermal  and  adiabatic  bubble  process 
respectively.  Note  that  the  actual  tiM  corresponding  to  the  particular 
values  of  n • r/*(t)  specified  is  not  of  intsrast  here.  Ton.  isothamel 
pressure  distributions  are  plotted  in  Figure  A-3  as  functions  of  C = r/R(t) , 
and  a similar  aet  of  adiabatic  pressure  distributions  is  plotted  in  Figure 
A-4.  In  these  diagrams,  the  pressure  maximum  die.ussed  by  Rayleigh  is 
clsarly  shown,  and  the  distribution  at  the  instait  of  full  collapse 
la  given  by 


♦(W?> 


max 


(A-6) 


It  should  be  mentioned  that  the  C-coordinete  is  stretched  with  time, 
since  R (t)  decrease;  in  the  process  of  collapse. 

The  pressure  profiles  resemble  in  shape  those  obtained  by  Hickling 

end  Plesaet10  (for  collapse)  for  compress^  >le  liquid.  They  demonstrate 
that  within  e distance  of  Rq  from  the  bubble  center  , pressure  in  the  liquid 
can  reach  values  several  orders  of  magnitude  abov*  PQ. 


vs.  the  initial  bubble  pressure  ratio 
.1  radii. 


0.  2.  4.  6.  8.  10.  12. 

5 - r/R(t) 

Figure  A-3 . Pressure  distributions,  $ - p/PQ-l  in  the  liquid  phase  vs. 

dimensionless  radial  distance  for  isothermal  buJcble  at 
several  stages  during  collapse. 
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APPKNPIX  B 

Approximate  Analys  la  of  Liquid  C jprssalblllty 

Hi*  object! ve  of  the  present  calculation  is  to  obtain  an  estimate 
of  the  effect  of  liquid  compressibility  upon  the  maximal  pressure  and 
temperature  of  a collapsing  gam  bubble.  To  carry  out  the  calculation 
without  solving  the  entire  liquid-bubble  dynamic  problem , certain  simpli- 
fications were  necessary.  The  strongest  assumption  made  herein  was  that 
the  pressure  distribution  in  the  liquid  surrounding  the  bubble  is  approxi- 
aately  equal  to  that  obtained  from  incompressible-liquid  analysis  of  a 
collapsing  bubble,  at  the  instant  of  full  compression,  !•••»  when 
0 » dR/dt  - 0. 

The  ensuing  calculation  is  based  on  the  compression  work  balance 
mentioned  in  section  7.3,  viz . , 


W m W + W 

STROKS  CL  CB 


(B-l) 


where  each  of  the  terms  is  evaluated  at  the  state  of  maximal  bubble 
compression . The  above  equation  will  be  shown  to  depend  upon  the 
final  state  R/R  ■ n , of  the  bubble,  and  will  have  to  be  solved 
iteratively.  ° 


The  results  of  this  approximation  analysis  are  presented  in  tabular 
farm  at  the  end  of  this  appendix,  as  well  as  in  Figures  B-2  and  B-3, 
where  comparisons  with  the  incompressible  case  (adiabatic  bubble)  are  made. 


B.l  The  pressure  Distribution 


FOr  purposes  of  the  present  calculation,  only  the  state  of  maxi me 1 
bubble  compression  is  of  interest , when  U ■ dR/dt  » 0.  The  pressure 
distribution  in  th*.  liquid  is  approximated  by  that  of  tJbm  incompressible 
liquid  case  (Rayleigh  solution).  At  tbs  instant  when  U * 0,  Bernoulli's 
equation  yields. 


Vn) 

c 


(B-2) 


where  p • p(C,n) , C * r/R(t)  end  n * R(t) /R.  p is  the  uniform  initial 

o o 

driving  pressure  in  the  liquid,  and  pg (n)  is  the  (uniform)  bubble  pressure. 

Ifcen  0 - °,  n - n^,  R ’ 'via  PB  " % 

XIX 

Considering  the  adiabatic  case. 


(B-3) 
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wham  y - C /C„  for  the  bubble  gas,  and  hence 
P v 


with  n ■ w 


(»-4) 


B.2  Hie  Caagreasion  Stroke  Work 

Mhen  liquid  coapressibility  ie  considered,  the  net  change  in  volune 
for  a full  coepression  stroke  (for  a single  bubble)  involves  the 
sue  of  bubble  and  liquid  c expression,  such  that  the  stroke  work  is  defined 


W - -pol4v»  * 4V' 


(»-S) 


The  bubble  volusw  change  is 


Av  - ir(RS  - Rs)  - i iru*  (n*  - 1) 

B O O JO 

where  n ■ nm1n  • Determination  of  the  liquid  vol 
■ore  involved,  and  undertaken  as  follows. 


(»-6) 

change  is  considerably 


The  variation  of  voluew  is 


<5vl  - 6 (M/p)  - 6m/p  ♦ MS  (1/p)  (b-7) 

where  N is  the  ease  of  a liquid  volusw  unit  undergoing  caepresoion,  which 
rwains  invariant.  Thus , - 0 and  the  voltssa  variation  i» 

«VL  - M6  (1/p)  (•-€) 

Considering  a differential  shell  of  radius  r end  thickness  dr, 

N » 4ipr2dr 


and  the  total  liquid  voluee  change  is  defined 


AV, 


I 


Pjr) 

4irr*drp  / d(l/p)  - 4w 


l 


r2dr(l  - p/p J 
o 


(■-9) 


The  integrand  in  the  last  equality  vanishes  identically  for  the  incoapressible 

jease,  m p » p is  unifoxu.  further,  as  r — ► m the  integrand  vanishes 
o 

aeyep  tot  icall  y p — ► oq,  but  a node  rate  driving  pressure  PQ  is  required, 

such  that  p/p  — ► 1 faster  then  1/r2  » 0. 

o 
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Uti lining  the  liquid  «— praaaibility  axpraaaion, 
C<jn  » p ♦ B 


(»-10) 


( for  water,  B ■ 3000  ata  and  n ■ 7 ia  proposed  aa  good  approximation 1 ') , 


" <rri) 

o 


^ _ l/h  1/h 

Bi  - (A  * i) 


aha  ra 


♦ - P/»,  - 1 

and  B'  ■ (1/P  ♦ 1) 

o 


Tr ana forming  to  dimanaianlaaa  coordinataa , 5 • r/R,  and  using  tha  praaaura 
distribution  given  by  tbs  incompressible  aolution, 


Avt  - «*y  jr  [ 1 - (t  ♦ 


(»-ii) 


where  a * 1^/1'  and  ♦ • tg  ia  defined  by  Eq.  (1-3)  > tha  notation 

r^irt  applies  aa  anil,  and  will  ba  oadttad  haraon.  Caana  of  intaraat 

involve  Pq  ~ 1000  ata,  and  if  tha  trand  ahoan  in  Appendix  A for  incompressible 

liquid  ia  correct,  ia  ej^ected  to  bn  large,  Thnmforn,  a » 1 can 

bo  anticipated.  To  avoid  tha  aingularity  at  infinity  involved  in  tha  last 
integral,  a convenient  diannaionlaaa  coapraaaihility  radios  is  dafinad  by 
a - Tjj/S* . sad  liquid  onymaaion  is  cons  ids  md  only  within  tha  bounds 

!<(<«.  With  tha  following  change  of  variable. 


x - C/» 


tha  problem  —cunts  to  calculating  tha  dafinita  intngral 


(1-12) 


1.  Tha  praaaura  distribution  and  region  of  integration  an 
itically  in  Figure  i-l.  Thus, 


= .■/[. 11  **> 
■U  »i/n 


1/n  1 

J 


x2dx 


Sines  x < 1,  tha  tan  in  brackets  can 

sarins t 


ba  evaluated  by  use  of  a binomial 


12a 


123 


(1  ♦ x)*  • ^T^(*)x*  • 1 + qx  ♦ x*  + OU1) 


taken  up  to  third  order  accuracy.  Hence, 


f r 2 2-1 /n  X3‘1/n  n(n  _1)  4-1/n  1 . 

X*a  J ^ -x  “ — -IT * J** 

l/a  . 

3 3 [l  £±  x1'1^  r £ ~ *>  2-1/nf 

“ LT  ~ 3-1/n  “ n(4-lA0  " 2(5-l/n)  * J1/a 

I - -0.041a*  (®-13) 

where  values  at  tha  lower  1'jd.t  of  Integration  vara  assumed  negligible, 
as  1/a*  « 1 is  axpactad.  Itw  atroka  work  can  now  ba  calculated,  with 
tha  evaluation  of  the  liquid  compression  tens: 


Avt  <•  - 4ir**n 

U O 

hence,  from  Sq.  ®»-5  ) , 


rvn)  1* 

■ « oW  - - if  *’n’  « o.u3 


■ "o'o  Jd-n’l  ♦ o.U3n’  | 


Vh*r*  n ’ Vn 


(B-14) 


B.3  Tha  Ccaaprassion  work  Terms 


Tha  bubble  compression  work  is  evaluated,  assuming  rapid  process, 
by  use  of  the  adiabatic  relationship: 


V ° 

•J  <r*  ro  ‘ty*  4™2"  " 

o n 


« (!s,  i-ji.-i’r1 

3 o o'*  Y - 1 

o 


**rm  n - n^. 


The  liquid  compression  work  is  given  by 


(1/P) 


T Vr> 

• - I 4irr2dro/  pd 


(■-15) 


(■-16) 
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The  internal  integral  can  ba  readily  evaluated,  using  tha  liquid 
coaqpression  isotherm 


Pn  1 


(»-17) 


Note  that  in  the  last  expressaxon. 


Cp  “P  + I ■ P B'i  so  that 
o o o 


p^pd(l/p)  - 


. ifftta  . . T .(£*j,l/nl 

n-1  P +B  lP  +B  I B lP  +T  | 

. ° o j L o J 


Therefore,  with  tha  Sana  change  of  variable  of  integration  as  in  Bq.  vb-11) 
for  AV_  , 

Is 


w - ♦ 4m 

CL  o 


on'Po  ] C*« 


•kb'  - i)  J [i  - <£  ♦ i)1/n  1 c*dc 

1 L J 


ia-18) 


Limiting  tha  integration  to  the  interval  1/a  < C < 1 and  further  tranafoi 
of  the  integration  variable  x = C/a » 


"cl  * 4’*in’ 


- dirm^’p 
o o 


po  / •'»**  * & 1 * 


l/« 


n*l  • 1 


’1/e 


where  I was  defined  by  Bq.  (B-13)  • Assuming  that  1/a*  « 1 , the  liquid 
compression  work  term  can  be  written 


" 5«*1n,p 
CL  3 o o 


1.51'  - In(B'-lHl)  0.123 


n-1 


(B-19) 


The  toecqy  Balance 

TtM  iMr^y  btlinc*  at  ms  latl  bubble  xmipressian , 
H_ 


- M (ADIABATIC)  ♦ W 
O CL 


can  be  now 


in  tanas  of  n - r^4|t  by  use  of  Bqs. 


(•-14) , 


(•-IS),  and  ft-19)  i after  substitution  and  dividing  through  by  4/J»R 9 , 

o o 


where 


a(n) 


n*  «•  o.uinV  - 


-3(^-1) 


■>] 


♦ i^nV 


(B-20) 


B’  - b/f  ♦ l 

O 


-A-^LaW^tU 

B n - 1 


Equation  (B-20)  can  be  solved  iteratively  to  determine  n ” n-4rt»  Hw 

method  chosen  for  this  purpose  is  Begula  Falsi  ( or  method  of  falsa 
position) . 

For  incompressible  liquid,  it  has  been  demonstrated  in  Appendix  A 
the  the  final  compression  ratio  n_<  ^ for  both  adiabatic  and  isothermal 

cases  depends  only  on  the  initial  pressure  p°/PQ  (and  y for  the  adiabatic 

bubble) . The  present  energy  balance  equation  for  the  compressible  liquid 
shows  that  depends  in  addition  upon  Pq,  (aa  well  as  B and  n which 

contain  the  liquid  compressibility  properties) , since  the  amount  of 
compression  work  dona  on  the  liquid  depends  on  the  actual,  value  of  the 
driving  pressure. 


How  that  n^4tt  is  t. reliable , the  normalised  (total)  bubble  compression 
work  can  be  calculated, 


WCB  " "cB^STBOKS 


(7.-n*40.123nJa*) 


(B-21) 


and  the  normalised  overall  liquid 


avion  work  is 


*CL  " *o/*S 


- 1 - w 


CS 


(B-22) 


the  maximal  bubble 
in  *q.  (B-3), 


is  defined  by  use  of  the  calculated 


pressure 


"o  - <vv*ej. 


and  the  final  bubble  taaparatura  ratio  .a  calcuiatad  throw? 
adiabatic  relationship, 


o 

(Pj/P,)  Y ■ 


i/n,(Y‘1) 


tha 


deploying  values  of  • » 3000  atm  and  a «_7  (for  water)  and  y - 1 
(for  air  bubble),  the  parameters  n-in*  ♦a>|U(>  and  Tj^r^/T^  ware 

calculated  for  (to  a ad  (P^  j p°/»()  data,  nurse  distinct  values  of 

**o  " 100  , 500,  and  1000  atm  ware  employed.  Itw  results  are  tabulated 

and  plots  of  and  versus  p°/PQ  are  shown  in  Figures  S-2  and 

S-3  respectively,  compared  in  each  case  with  tha  incampresaibla  liquid 
adiabatic  bubble  solution a . 
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4 

9 
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7 

* 

1# 

a 

a 

a 

14 

19 

U 

tr 

it 

it 

29 

21 

22 

23 

24 
IS 


-~s/ro 

.i »»»«*-*) 

aiMii-n 

• 24»ut-fa 

IIUIH] 

• 

• 4»m»-n 
tnttt-n 
Misn-fi 

• Hflll'll 

. 119999-91 

. tiniMi 
ntUMl 

• 1944 ll-#l 
.412*41-91 
.429399-91 
.•12929-91 
. 199999*99 
.139919*99 
.199339*99 
.299*39*99 
.344129*99 
•491239*99 
•••*499*99 


**»•*»  .A02au*rzc 


KM  <99/99*99) 
.199999*91 
. 199999*91 

•114999*91 

.12499%. 91 

.132991*91 

.149999*91 

.149991*91 

.194999*91 

.144999*91 

•173r»»*91 

• 1 *99917*91 

.109999*91 

.13*999*91 

.294999*91 

.242991*91 

.239991*91 

. 229999*91 

.33*999*91 

. 344991*91 

. 292999*91 

.2*9999*91 

.311999*91 

-7  7*99**91 

.304999*91 

.292999*91 


m/to-t 

.491*91*9* 

.221379*9* 
. 1 33299*9* 
.7*9111*99 
.442999*93 
.299449*99 
•147439*99 
.091919*94 
.492249*94 
70*099*94 
. *4 9199*94 
.999499*93 
.997499*93 
.329999*93 
.199999*93 
.111499*93 
.499199*92 
.394329*92 
.220139*92 
.134949*92 
.790999*91 
.493379*91 
.259321*91 
124139*91 
. 400999*99 


<r  - i.40) 


>391199*93 
.391379*93 
.191449*93 
139421*93 
.193739*93 
.•31779*93 
.•44919*92 
.934971*92 
. 429991*92 
.347949*92 
. 3794*9*92 
. 229991*92 
.177911*93 
.141919*92 
.113749*92 
.912919*91 
.731149*91 
.*••141*91 
. 449991*91 
.374791*91 
. 392991*91 
.242191*91 
.194991*91 
•199499*91 

. i24*4».jn 


conraMmj  uono.uutttic  mm 

a " 3000  «ta,  n - 7,  y m 1.40 


Pe  - 100  «ta 


a 

W/fO 

n/a* 

si/ro-i 

NCX./W4TM 

Nut.n/n* 

1 

1.999991-9) 

.31*191-91 

. 194411*94 

. 924991*99 

.424311*92 

2 

. 129191-92 

.342*71-91 

.179791*94 

.999211**9 

.572491*93 

3 

. 139491-92 

.349491-91 

.144991*94 

.194118*99 

.923121*93 

4 

. 199939*92 

. 39**9|-91 

.1399*1*94 

.474391*99 

.477948*92 

S 

.291191*92 

.439431-91 

.134911*94 

.199231*99 

.439998*92 

4 

.314231*92 

.449371-91 

. ! 244 31*94 

.•39331*9* 

. 394171**2 

7 

.399111*92 

.993471-91 

.112991*94 

.*99441*99 

.341111*92 

• 

.991191-92 

.449491-91 

.191211*94 

.744*31*99 

. 321921*92 

9 

.439*41-92 

.9919*1-91 

.993271*93 

.722131*99 

.237)21*92 

1# 

.794332-92 

.449971-91 

. 791411*93 

.479891*99 

.24*741*92 

13 

1.999992-92 

.792441-91 

.494771*93 

.499111*99 

.242941*92  • 

11 

.129991-91 

.749*91-91 

.397711*93 

.929391*99 

.214941*92 

13 

. 13S49B-91 

.•4*131-91 

.999721*93 

.448479*99 

.193131*92 

14 

.199921-91 

. 941*21-91 

.499741*93 

.349441*99 

.179311*92 

IS 

.291191-91 

.199741*99 

. *'13948*93 

.245191*99 

.14*231*92 

14 

. S1423S-91 

.129411*99 

.22*721*93 

.148391*99 

. 124*21*92 

17 

. 314111-41 

.139228*99 

. 194941*93 

.749171-91 

.194321*92 

IS 

.991191-91 

.142991*99 

.! 91991*93 

• 33979C-91 

.*•1*91*91 

19 

.439*41-91 

.191*21*99 

.43*421*92 

.141498-91 

.729321*91 

11 

.794331-91 

.229791*99 

.491391*92 

.S791*r-J2 

.594431*91 

21 

.199990*99 

.249921*99 

.234398*92 

.2492)1-92 

.492121*91 

12 

. 129*99*99 

.399731*99 

. 142921*92 

.191948-92 

.49*141*91 

23 

. 194491*99 

.349191*99 

.199991*92 

.444291-93 

.349341**1 

24 

.19933**99 

.414421*99 

. 48 *99l«f 1 

.194791-93 

. 2*S94aa*91 

IS 

.29119**99 

.479291*99 

491441*91 

. (49731-44 

.24179***1 

24 

.314231*99 

.341491*99 

.'•4241*91 

.394971-94 

.293431*91 

P0  ■ 500  ata 


1 

989/10 

n/M 

91/90-1 

NCX./WSTH 

mx.n/TM 

i 

1.9*9*91-93 

.992331-91 

.214441*93 

.997931*99 

.342991*92 

2 

.123*98-92 

.349341-91 

.243741*93 

.199391*99 

.331741*92 

3 

.131498-92 

.3*1421-91 

.243431*93 

.942491*99 

.393491*92 

4 

.199931-92 

.424241-91 

.224798*93 

.913111*99 

. * >7928*92 

5 

.23119I-S2 

.4Y4391-91 

.297191*93 

.922431*99 

.234238*92 

4 

.3142)1-93 

.724771-91 

. 1994*1*9! 

.999141*99 

. 232438*93 

7 

.391111-92 

.713471-91 

.174*41*93 

. 19 5141*99 

.212411*92 

• 

.991191-92 

.*49991-91 

.149191*93 

.•77171*99 

. 193978*92 

9 

.439941-92 

.912191-91 

. 144178*93 

.*37991*99 

. 174998*92 

19 

.794331-92 

.983278-91 

.132991*93 

.933571*99 

. 141341*92 

11 

i.99*num-*3 

. 194341*99 

. 129411*93 

.*93231*99 

. 144*91*92 

12 

.123*91-91 

.119348*99 

.191941*93 

.771421*99 

.133948*92 

13 

.131491-91 

.123941*19 

.971*31*92 

.731791*99 

.121111*92 

14 

.199931-91 

.139171**9 

.*42748*92 

.4*4221*99 

.199711*92 

13 

.291191-91 

.147918*99 

.757721*92 

.427131*99 

.999311*91 

l< 

.314331-91 

.141411*99 

.494241*92 

.941291*99 

.*99491*91 

17 

.391111-91 

.177411*99 

.997991*92 

.413471*99 

.7949*1*91 

11 

.391191-91 

. 199771*99 

.442828*92 

.399211*99 

. 797191*91 

19 

.439948-91 

.217431*99 

.371991*92 

.299991*99 

.423391*91 

29 

. 794)31*  91 

.244191*99 

.284198*92 

. 294941*99 

.942991*91 

21 

.199991*99 

.274821**9 

.219181*92 

.1299*1*99 

.447991*91 

22 

. 123*31*99 

. 314390*99 

.14772**92 

.4*1991-91 

.397941*91 

23 

.1914)1**9 

.343991*99 

.199431*92 

.33*931-91 

.334371*91 

24 

. 1999)1*99 

.418341**9 

.473*21*91 

.149991-91 

.2*4311*91 

13 

.291191*99 

.4*9271*99 

.444741*91 

.749341-92 

.241111*91 

!4 

.314231*99 

. 94**81*99 

.392*18*91 

.332271-92 

.299411*91 

COMPRESS I <L£  LIQUID,  ADIABATIC  BUBBLE  (con't 


pD  • 1000  atin 


f 

I 

i 


* 

1 

2 

s 

4 

3 

« 

7 

• 

9 

14 

11 

12 

13 

14 

15 
1« 
17 
«,! 
19 
II 
21 
7 » 

2 3 

24 

iC 


P19/PO 

l.lfllll-13 

. 123991-92 
.139491-92 
199331-92 
.231191-92 
.311231-92 
. 399119-92 
.391191-92 
.439942-92 
.794331-92 
1 .999999-92 
.123999-91 
.139499-91 
.199939-91 
. 2S119W-91 
.314231-91 
.399119-91 
.391199-91 
.439949-91 
.794331-91 
. 199999*99 
. 1 2 3UV 
. l5849C*tJ0 
.199339*40 
.251199*99 
.314239*99 


99/99 

.394232-91 
. 641379-91 

.499939-91 

.742342-91 

.799311-91 

.41*411-91 

.924939-91 

.999741-91 

. 197472*99 

.114132*99 

.123339*99 

.133499*99 

. 144422*99 

.139921*99 

.171911*99 

.144729*99 

.293339*99 

.222199*99 

. 243349*99 

.249479*99 

.247949*99 

. J 3 2->31>0-t 

. 37517E*00 

.42547 9*09 

.4841 JE*44 

.339899*99 


P9/PO-1 
.138949*93 
. 127879*43 
.114329*93 
. 199349*93 
. 191991*43 
.931479*42 
.437741*42 
.784479*42 
.723249*42 
.441739*42 
.493431*42 
. 5444-2*42 
.444439*42 
.444779*42 
.399339*92 
.333941*92 
. 3 14*18*42 
.249391*42 
.227911*92 
.194971*92 
.131911*92 
. 1 1 :>u2*iz 
. 5J7  3 30E*U1 
.622J5E**1 
. 42363E*01 
.286829*41 


MCL/WSTRK 
.943642*44 
. 964242*44 
.934429*44 
.946749*44 
.934341*44 
.923349*44 
.917132*44 
.943429*44 
.444232*44 
.474414*49 
.8416 11*44 
.823442*44 
.793742*49 
.734762*44 
.717391*49 
.664911*44 
.614941*44 
.343691*44 
.466241*94 
.379441*44 
.289412*49 
.1937 ?n*fff 
. l.!4  2 8E*0U 
.694452-41 
. 3 54  39E-41 
.168342-91 


! 

I 


MAX . TB/TBr 

.294792*42 

.274472*92 

.247371*42 

.226322*92 

.247372*42 

. 189772*92 

.173691*42 

. 134739*92 

.143441*42 

.132439*42 

.124991*92 

.114169*42 

. 149342*4: 

.911932*41 

.827499*91 

.749172*41 

.676342*41 

.649349*41 

.944392*41 

.494332*41 

.427749*41 

. 3 742jE*J1 

. 3:4  29E*vn 

.279842*41 

. .13881E  + 01 

.244512*41 


(10*R  . /R  ) 
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APPRWDIX  e 


rum  vtsoalizatiow  tbstrr  n&aMMUt  prootocrx 

The  following  part*  require  "0" -rings  prior  to  installation) 


PART 

Pneumatic  Piston 
Separator  Piston 
Projectile  Piston 
Septus  Access  Plug 


”0"  RING  SIZE 

♦2-012-H552-70 

I2-011-MS52-70 

#2-011-N552-70 

#2-006-KS52-70 


QT5f 

2 

2 

2 

1 


Install  the  following  items  into  the  pneumatic  loading  L.P.  cylinder 
and  visualisation  chamber  in  order  as  shown  (reference  Plow  visualisation 
Tester  Assembly  Drawing) . 


1.  Pneumatic  Piston  with  "o"  rings. 

2.  1/4  tube  to  3/8  pipe  fitting  (cover  pipe  threads  with 
teflon  tape  prior  to  installation) . 

3.  1/4  tube  to  1/8  pipe  fitting  (cover  pipe  threads  with 
teflon  tape  prior  to  installation) . 

Secure  pneumatic  load  L.P.  cylinder  to  tester  stand  (not  shown). 

4.  Silicon  rubber  septum. 

5.  Septum  retainer  (cover  with  teflon  tape  prior  to  installation) . 

6.  Septum  access  plug  (coat  "0"  ring  and  end  of  plug  with 
vacuum  grease) . 

7.  Poppet  valve  seat  with  "0"  ring.  Ensure  "oM  ring  is  properly 
seated  in  poppet  valve  seat.  (Coat  tcp  sealing  surface  of 
seat  with  vacuum  grease.) 

8.  Drop  poppet  valve  thru  poppet  valve  seat  and  pneumatic  load 
L.P.  cylinder. 

9.  Install  flow  guide  into  visualization  chamber  (coat  top  sealing 
surface  with  vacuum  grease  prior  to  installation) . 

10.  With  visualization  chamber  oriented  as  shown,  place  chamber 
on  pneumatic  L.P.  cylinder  and  secure  with  (4)  1/2  - 13  UNC  x 
2"  long  high  strength  bolts  and  washers.  Torque  bolts  to 

15  in-lb  using  alternate  bolt  tightening  method.  Visually 
watch  sealing  surface  of  poppet  valve  seat  to  ensure  correct 
sealing  in  visualization  chamber. 

11.  Coat  Parker  "0"  ring  P/N  2-006-N552-70  with  vacuum  grease  and 
push  onto  poppet  stem.  Tighten  poppet  packing  gland  into 
pneuswtic  L.P.  loading  cylinder. 

12.  Install  poppet  valve  spring  into  countersink  in  poppet  packing 
gland  and  hold  in  place  using  (2)  #5-40  washers  and  (2)  #5-40 
nuts. 

13.  Ieolaticn  piston  with  "0"  rings  — install  piston  into  visualization 
chamber  with  chamfered  end  in  first.  Care  should  be  taken  during 
installation  not  to  slice  "0"  rings  while  pushing  past  threads 

in  ch— her . (Also  a requirement  for  installing  projectile 
piston.) 


Projectile  piston  with  "O"  rings. 

Visualisation  chaster  plug  — screw  in  plug  until  threads  are 
flush  with  outs  ids  edge  of  visualisation  cheater. 

Pressure  transducer  — install  transducer  with  adapter 
and  copper  washer  until  snug. 

Bleed  valve  stew.  (Cover  thread  with  teflon  tape  prior  to 
installation.) 

1/4  tube  to  1/S  pipe  fitting  (cover  pipe  threads  with  teflon 
tape  prior  to  installation) . 

19.  L.P.  drain  line. 

20.  Install  d2  gas  line  and  L.P.  fill  lina. 


t 


tot 
be  taken 
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Wtm  VISUALIEKTTOS  TPTIH  FILLING  PKOCEDW 

filling  Flow  Visualisation  Taster,  the  following  stops  oust 
(roforonoo  Flow  Visualisation  Schematic  Drawing) j 


Kosovo  M2  gas  lino  fron  solo no id  valvo  and  onsuro  pnousatic 
piston  is  In  propor  location.  Reconnect  gas  lino. 

Ensuro  all  linos,  valves,  and  fittings  aro  properly  tightened. 

Close  all  valvos  an  tost  fixture  and  rosovo  air  compressor  lino. 


Open  regulator  valvo  on  N2  gas  cylinder  and  sot  gas  pressure  as 
per  tost  requirements . 

Open  blood  valvo  stem  to  L.P.  drain  lino. 

Using  positioning  rod,  sot  projectile  piston  at  starting 
location  (i.o. , roar  edge  of  piston  in  lino  with  edge  of  fill  hole) 

Sot  poppet  spring  compression  (distance)  as  follows: 

for  300  psi  tests  ■ 21/32  total  compression 
for  500  psi  tests  - 19/32  total  coapreecion 

Fill  L.P.  reservoir 

Connect  L.P.  fill  hose  to  L.P.  fill  valve. 

Open  L.P.  reservoir  stopcock. 

Open  L.P.  isolation  valvo. 

Manually  hold  poppet  valvo  open. 

Slowly  open  L.P.  fill  valve  until  a slow  flow  is  detected. 

Natch  L.P.  re so voir  and  refill  when  necessary,  whan  l.p, 
is  seen  coming  out  the  drain  line  and  no  air  bubbles  emerge, 
close  L.P.  fill  valve. 

Move  poppet  up  and  down  to  dislodge  any  air  bubbles. 

Close  bleed  valvo. 

Close  stopcock  on  L.P.  re so  voir. 

Remove  L.P.  fill  hosa  from  L.P.  fill  valvo  and  insert  hose  plug. 

ullage  addition (if  no  ullage  go  to  step  24) . 

Open  L.P.  fill  valvo. 

Kosovo  septum  access  plug. 

Load  syringe  with  required  amount  of  gas. 

Insert  syringe  thru  septus  and  inject  gas  into  chamber. 

Kosovo  syringe  and  install  septus  plug. 

Close  L.P.  fill  valvo. 


24.  Close  L.P.  isolation  valvo. 

25.  Chock  that  all  valvos  aro  closed. 

26.  Ensuro  contact  strobe  wire  in  chamber  1 

27.  initiate  Visualisation  Vests  Checklist. 


plug  is  in  proper  position. 


VISUALIZATION  TESTS  CHECKLIST 


8 8 8 8 ft 


i 


PRE-FIRING! 


□□□□□  1. 
□□□□□  2. 

□ □□□□  3. 

□ □□□□  4. 

□ □□□□  3. 

□ □□□□  6. 

□ □□□□  7. 

□ □□□□  8. 

□ □□□O'*. 

□□□□□ll. 

□ □□□□l2. 


After  filling,  ensure  all  valves  and  vent  plugs  ara  closad. 

Closa  LP  rasarvoir  atop,  disconnact  fill  hose,  and  ranova 
LP  rasarvoir  from  tast  araa. 

Dump  any  laftovar  LP  froa  drain  bottla  into  usad  LP  bottla 
and  put  drain  hose  back  into  bottla. 

Ensura  stroba  light  is  pluggad  in.  Thst  fira  sam. 

Unlock  camera  film  stop  and  advanca  film.  Ensura  Ians  covsr  is  off. 

Re -check  caaara  focus  using  focusing  rod  and  ra -check  "fH  stop. 

Closa  installation  covar  and  sat  charge  aof>  to  "OPE"  mode. 

Place  Blomation  in  ready  state.  Check  scale  (x:l)  and 
baseline  of  trace  lines. 

Turn  off  laboratory  lights. 

Open  caaara  shutter. 

Turn  on  lab  warning  light. 

Go  to  control  panel,  insert  key,  arm,  and  fire. 


POST  FIRING; 
Remove  key. 

D □ □ D D 14 . Close  c« 


ira  shutter. 

□ □□□□«.  Sat  charge  aap  on  ”GND". 

pneumatic  bleed  valve 
I — i p— 1 1 — 1 1 — (-—I  bleed  valve . 

LJ  LJ  L_I  LJ  LJ 17 . Disconnect  stroba  light. 

□.□□□□*•  Turn  off  lab  warning  light 


zasat  pneumatic  piston  and  open  chamber 


1.  Nashya  m*7. 

2.  Lanai  50mm 

3.  filat  Kodak  FXF-120  • a W,  ASA  125 

4.  Lana  at  "T"  tlaa  aatting  and  "f"  • top  sat  batwaan  fll  and  fl6. 

5.  Front  adga  of  Ians  placad  approsl  nataly  3"  front  can  tar  of 
Visualisation  Chaatbor  bora. 

Stroba  < 

1.  Vi vital  Auto  252 

2.  "Mad”  flash  aatting 

3.  Faoa  of  atroba  placad  approximately  6"  front  center  of  Visualisation 
Chanbar  bora.  Stroba  placad  at  45  da  gras  angle  to  front  faca 

of  visualisation  Chanbar 
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BIONATION*  Nodal  1015  wawa  for*  raoordar. 


IMPOTSs 

1.  Ch  1;  chaabar  praaaura  f 5 VFg 

2.  Ch  2 i flaah  discharge  voltage  • SO  VPS 

3.  Both  input*:  DC  nod* 
input  Nodus 

1.  Ch.  1 and  Ch.  2 (2044  word*  par  channal) 
Bnoord  Notes  Dalayod 


triggers 

1.  • 300  pal  - 200  words  dalayod 

2.  • 500  pal  ■ 370  worda  delayed 

3.  External  aourca 

4.  Single  awaap  aoda  1 

5.  D.C.  coupling 

6.  Poaitioa  alopa 

7.  trigger  level  at  9 o'clock 

Tlaa  Baaas 

1.  Dual 

2.  • 300  pal,  tiaabaaa  'A*  • 0.2  aa/vord 

"B"  • 0.05  as/word 

3.  • 500  pal,  tiaabase  "A"  • f 2 aa/word 

"B"  # C.02  aa/vord 


CHANGE  AMPLIFIERS 

I SPOT s 

1.  chaabar  praaaura  tram  "Kriatal"  praaaura  transducer,  typo 
601  A a/n  62094. 

2.  sensitivity  • 1.075  pc/pai 

3.  ran 9*  * 100  pai/V  for  300  psi  tost 

500  pai/V  for  500  psi  tast 


appctdix  d 


STAPfKP  OUUaOI  TKETKP  AJEEMBLT  TKXXCUT* 


The  following  porta  require  "O"- rings  prior  to  installation t 
PART  "0"-PIWG  SI A QTT 

STAkTXP  CHARGE  #2-026-KS52-90  1 

CHAMBER 

PRIMP  GAS  #2-020-N552-90  1 

XIRSGPATZNG 

CHAMBER 

Da  far  to  tha  Star  tar  Charge  Taatar  Aasaably  Drawing  for  tha  following 

assembly  procedure.  Nhara  loading  of  smokeless  powdar  and  wadding  ia 

mentioned,  rafar  to  tha  startar  charga  taat  matrix  for  tha  typa  and 

quantity  of  aach. 

1.  Install  alactric  primar  into  Priaar  Gas  Intagrating  Chaahar. 

Aftar  installation,  claan  contact  surfaca  of  priaar  and  faca  of 
Intagrating  Chamhar.  Using  a volt/oha  matar,  nota  rasistanca 
batwsan  priaar  and  Intagrating  Chaabar. 

2.  Load  powdar  and  wadding  into  In  tagrating  Chaabar  as  required. 

Ensura  powdar  doas  not  gat  inaida  thraads. 

3.  If,  in  stap  2,  tha  Intagrating  Chambar  was  fully  loadad  with 
smokeless  powdar  (i.a.,  2.60  grams),  covar  tha  four  holes  on  tha 
insida  faca  of  tha  Orifica  plug  with  callophana  tap*  (no  tapa  is 
raquirad  for  a partial  loading) . Covar  thraads  of  tha  Orifica 
Plug  with  graass  and  install  into  the  Intagrating  chamber.  Tighten 
plug  with  tha  special  spannar  tool.  Cover  outside  face  of  Orifice 
Plug  with  callophana  tapa. 

4.  Coat  all  outside  sealing  surfaces  of  Integrating  Chambar  with 

grease  and  drop  into  and  of  Starter  Charga  Chambar  as  shown.  Carefully 
push  Intagrating  Chamber  in  until  tha  "0"-ring  contacts  the  insida  bora. 
Using  a wooden  dowel,  carefully  tap  tha  Integrating  Chambar  in  until 
it  is  properly  seated  in  tha  starter  Charga  Chambar. 

5.  Claan  tha  contact  pin  and  and  Murface  of  tha  Electrical  Primer  Contact 
Assembly.  Coat  tha  thraads  of  the  Contact  Assembly  with  grease  and 
install  into  tha  Startar  Charga  Chamber.  Using  a volt/ohm  mater, 
measure  tha  rasistanca  between  canter  pin  of  electrical  connector 

and  casing.  If  rasistanca  has  increased  by  more  than  25,000  ohsw 
over  tha  reading  taken  in  stap  #1,  remove  tha  Electrical  Primer 
Contact  Assembly,  re-clean  it  and  re-clean  the  contact  surfaca  of  tha 
Integrating  Chaabar.  Re-install  tha  Contact  Assembly. 

6.  Turn  tha  Startar  Charga  Chamber  over  and  load  tha  chambar  with 
smokeless  powdar  and  wadding  as  required.  Ensure  chambar  doas  not 
rest  on  Electrical  Primer  Contact  Assembly. 

7.  Install  tha  Blast  Shi' Id  into  the  Blind  Ch amber,  oriented  as  shown, 
and  secure  with  (2)  # 6-32  UNO  x 3/4"  long  machine  screws  with  lock 
washers. 
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•.  Install  the  PCB  pressure  transducer  with  seal  ring  into  the 
Blind  Cbsdsr.  ‘targes  transducer  to  25  ft-lhs. 

non i Do  not  install  tawdar  Chamber  Want  Plug  at  this  tin*. 

9.  Grease  the  outside  coaling  surfaces  on  the  Starter  Charge  Chamtsr 
"neck"  and  *0--ring. 

10.  Carefully  push  the  Iliad  Ch— tar  down  onto  the  starter  Charge 
Chaster  "neck*  until  the  ends  of  both  parts  are  flush  against 
each  other. 

11.  Grease  all  asternal  threads  on  the  starter  Charge  Chaster  and  the 
Blind  Chaster. 


12. 


s 

E 

i 


Kith  Chanter  Ooupling  Collar  mounted  an  test  stand , put  tester 
asasatly  through  collar  in  direction  shown.  Screw  tester  assembly 
into  collar  until  collar  is  situated  evenly  over  both  parts  of 
tester  assembly. 

Grease  threads  on  Powder  chaster  vent  Plug  end  install  vent  plug 
in  sequence  as  show. 


t 

f 

p 


) 1 

t ' 

r * 


*4 


i 

k 


I 
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IS!!! 

i i & i £ sum  chMgjgan  cum  isi 

4l  4 4l  4l  41 

ll  l | l m=mm. 

[ H H ][  H ] 1.  Main  mire  that  the  safety  interlock  key  Is 

rawed  frost  tbs  remote  firing  panel. 

[ H H ][  )(  ] 2.  Short-out  the  canter  pin  of  tbs  pressure 

transducer  lead  to  the  connector  housing. 

f H ]f  ](  H ] 3-  Connect  the  pressure  transducer  lead,  asking 

certain  the  lead  is  out  of  the  path  of  the  sent  gas. 

the  firing  lead  to  the  Prismr  Electrical 


( H H H U ] 4. 
[ H H H H 1 s. 


Clone  the  safety  enclosure, 
insert  the  latch  bolt. 


the  handle  and 


( H H 

( )i  n 

l H U 
[ H H 
[ H H 

[ ][  H 
[ ][  ][ 


H H ] 4. 

][  H ] 7. 

H ][  ] s. 

][  )l  1 9. 

H H ho. 

H ][  in. 

If  H lia. 


Switch  charge 


to  "OPE" 


Ihke  oertain  that  channel  1 Record /Hold  button  on 
Slomstlcn  is  in  the  out  position  to  enable  the 
record 


Switch  Eioaation  to  "heady" 
Turn  on  lab  warning  light. 


Insert  safety  key  into  reaote  control  firing  panel 
to  activate  firing  circuit. 

Place  arming  switch  in  the  "on"  position. 

fire. 


PQSI-flRIftE 


I )[  ][  H H ]13. 

[ n n n u ]i4. 
t ][  n u )[  ns. 

[ H )[  ][  H Jit- 
r )[  )[  ][  H U7. 


he  turn  arming  switch  to  the  "off”  position, 
hemowa  safety  interlock  key  from  the  firing  panel. 
Push  in  channel  1 Record /Hold  button  on  llonatlon. 
Switch  charge  amp  to  "GND". 

Turn  off  lab  warning  light. 
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APPIMDIX  E 

press ion- Ignition  Sensitivity  Tester  Assam My  Procedure 


The  following  porta  require  "o"  ring*  prior  to  installation i 


Pneumatic  piston 
Separator  Piston 
Projectile  Piston 
Septun  Access  Plug 
Aland  Valve  stent 
Starter  Charge  Chamber 
Primer  Gas  Integrating 


"0”  Ring  Sise 


2 -012 -**74-70 
2-011-M352-90 
2-011-M352-90 
2-006 -*674- 70 
2-013-MS52-90 
2-026-1*552-90 
2-020-11352-90 


QTT. 


Refer  to  the  L.P.  Coapress ion- Ignition  Sensitivity  Taster  Aesenbly 
Drawing  for  the  following  assembly  procedure.  Install  the  following 


■ into  the  Pneussatic  Load  L.P.  Cylindar  in  order  as  shown. 

1.  Pneumatic  piston  with  "0"  rings. 

2.  1/4  tube  too  3/8  pipe  fitting  (cover  pipe  threads  with 
teflon  tape  prior  to  installation . 

3.  1/4  tube  to  1/8  pipe  fitting  (cover  pipe  threads  with 
teflon  tape  prior  to  installation) . 

4.  Silicon  rubber  septus. 

5-  Septum  retainer  (cover  with  teflon  tape  prior  to  installation). 
6.  Septun  access  plug  (coat  "0"  ring  and  end  of  plug  with  vacuus  gj 


7.  Poppet  valve  seat  with  "0"  ring.  Bn  sura  “0*  ring  is  properly  seated 
ia  poppet  valve  seat.  (Coat  top  eealing  surface  of  seat  with  vaouus 


8.  crop  poppet  valve  thru  poppet  valve 
L.P.  cylinder. 


9.  Coat  Parker  "O'”  ring  P/M  2-006-NS52-70  with  vacua  g 
onto  poppet  stem.  Tighten  poppet  peeking  gland  into 
L.P.  loading  cylinder. 


tic  load 


10.  Install  poppet  valve  spring  into  oountersink  in  poppet  peeking  gland 
and  hold  in  place  using  (2)  95-40  washers  wd  (2)  #5-40  nuts. 
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ll. Mt  poppet  ape  lag  oanpcMaiaa  (diataeoe)  as  follows i 

for  ISO  pai  toata  ■ 3/4  * total  ounpreesion 

for  300  pel  toata  - 21/32*  total  ooaprasaloa 

for  500  pai  toata  - 19/32-  total  oonpreeaion. 


Znatall  tha  following  itaaa  into  the  I a.  Conpressian  Chndnir  in 
ardtr  a a a town. 

12.  Neparatar  Platon  w/*0”  rings — inatal1  platan  with  etna  fared  and 
la  fir  at.  ON  a woodaa  Aowal  roJ  far  poshing  platan  down  chanhar 
bora,  an  soring  that  platen  ia  properly  aoatod. 

13.  Projectile  Platan — inatall  platan  with  special  depth  tool  for 
proper  platan  alignment. 

14.  Inatall  tha  blast  shield,  oriented  aa  shown,  and  secure  with 
(2)  #4-32  one  X 3/4-  long  naohine  screws  with  lock  washers. 


as  faro  installing  any  other  parte,  cheek  tha  chanter  bore,  light  aneor 
porta,  and  pressure  transducer  porta  and  clean  out  any  foreign  natter.  Next, 
attach  the  Micro switch  Ntg.  Precfcst  to  ths  Disc  detainer  and  secure  with 
(2)  #4-32  jmc  X 3/14"  long  aachins  screws  with  flat  washers.  Noant  the 
Wicroewitch  with  Push  Nod  onto  the  Mounting  Prwcket,  oriented  es  shown,  and 
secure  with  (2)  #4-40  U.<C  x 5/8"  long  screws  end  ruts.  Bold  tha  Shear  Disc 
against  tha  and  of  tha  Diac  Natainar  with  tha  Microswitch  Push  Mod  passing 
thru  the  Shear  Disc  hole.  Adjust  tha  Push  Nod  so  that  l/N"  protrudes  pest 
the  Shear  Disc  to  allow  for  proper  switch  actuation. 

15.  Shear  Disc  - coat  totn  aides  of  the  Disc  with  vacuun  gross*  and 
inatall  into  chaster  seat. 

16.  Disc  Natainar  with  Microsuitch  and  bracket  — grease  all  external 
threads  and  install  into  chanter.  Wrench  tighten. 

17.  Grease  all  external  threads  on  tha  L.P.  Oanpreaslon  Chanter. 

16.  With  rhnater  Coupling  Collar  aauntod  an  tost  stand,  screw  tha 
Oopnnioa  Chantoer  into  the  collar  until  (1)  ana  tKr^ir*  is 
visible  on  rh  aster 

19.  Notate  tha  Goapreuion  nianher  until  tha  (3)  light  censor  porta 

era  an  top.  Install  ths  (3)  light  sensors  end  *tti*  to  30  tt.-lte. 

20.  Cover  ths  pipe  threads  on  tha  1/4  tube  to  1/*  pipe  fitting  with 

teflon  tape  and  install  the  fitting  and  tha  L.P.  Drain  as  shown. 

21.  Notate  the  Oonprassioc  Chanhar  as  necessary  «d  install  ths  (3) 

pressure  transducer  a with  seal  rings  and  tha  transducer  port 
plug  as  show.  Torque  to  25  ft.-lta. 
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22.  fcjtate  ths  Compression  Chamber  so  that  tha  L.P-.  Pneumatic  load 
Cylinder  mounting  hole*  (4)  are  on  top. 

23.  Plow  Guide  — coat  top  sealing  surface  with  vacuum  grease  prior 
to  installation. 

24.  With  tha  L.P.  Pneumatic  Load  Cylinder  oriented  opposite  from 
shown,  place  cylinder  on  L.P.  Compression  Chamber  and  secure 
with  (4)  1/2-13  tINC  X 2"  long  high  strength  bolts  and  washers. 
Torque  bolts  to  SO  ft. -lbs.  using  alternate  bolt  tightening 
issthod. 

25.  Rotate  parts  so  orientation  is  as  shown. 

26.  Bleed  Valve  stem  — Coat  threads  with  grease  and  coat  "nose" 
with  vacuum  grease. 

27.  Install  the  L.P.  Pill  Line  Valves  to  L.P.  Reservoir  and  the  N2 
gas  line  from  solenoid  valve  to  Pneumatic  Cylinder. 


Where  loading  of  smokeless  powder  and  wadding  is  mentioned  in  the 
following  section,  refer  to  the  test  matrix  for  type  and  quantity  of  each. 

28.  Install  electric  primer  into  Primer  Gas  Integrating  Chamber. 

After  installation,  clean  contact  surface  of  primer  and  face  of 
Integrating  Chamber.  Using  a volt/ohu  meter,  note  resistance 
beeween  primer  and  Integrating  Chamber. 

29.  Load  powder  and  wadding  into  Integrating  Chamber  as  required. 

Ensure  powder  does  not  get  inside  threads. 

30.  If,  in  step  J,  tha  Integrating  Chamber  was  fully  loaded  with 
smokeless  powder  (i.e.  , 2.60  grans) , cover  the  four  holes  on  the 
inside  face  of  the  Orifice  Plug  with  cellophane  tape  (no  tape  is 
required  for  a partial  loading) . Cover  threads  of  the  Orifice 
Plug  with  grease  and  install  into  the  Intagrating  Chamber.  Tighten 
plug  with  the  special  spanner  tool.  Cover  outside  face  of  Orifice 
Plug  with  cellophane  tape. 

31.  Ooat  all  outside  sealing  surfacas  of  Integrating  Chamber  with  grease 
and  drop  into  end  of  starter  Charge  Chamber  as  shown.  Carefully 
push  Integrating  Chamber  in  until  the  "0"-ring  contacts  the  inside 
bore.  Using  a wooden  dowel,  carefully  tap  the  Integrating  Chamber 
in  until  it  is  properly  seated  in  the  starter  Charge  Chamber. 

32.  Clean  the  contact  pin  and  end  surface  of  the  Electrical  Primer 
Contact  Assembly . Coat  the  threads  of  the  Contact  Assembly  with 
grease  end  install  into  the  Starter  Charge  Chamber.  Using  a volt/ 
ohm  mater,  measure  the  resistance  between  center  pin  of  electrical 
connector  and  casing.  If  resistance  has  increased  by  more  than 


? 


25,000  otes  over  the  reading  taken  in  step  *1,  raaoea  the 
electrical  Priaer  Contact  Aeeeadoly,  r 'lean  it  mA  re -clean 
the  contact  surface  of  the  Xntegrct'  aeber . ka- install 
the  Contact  ksaenbly. 

33.  Turn  the  Starter  Charge  rhaafier  tear  and  load  the  dunbar  * ith 
aaokeless  powder  and  wadding  vr,  required,  ensure  chsubar  does 
not  rest  on  electrical  Prlwer  Contact  Assembly. 

34.  Grease  the  outside  sealing  surfaces  on  tbs  Starter  Charge 
Ohsnhar  ’neck"  and  "O' -ring. 

35.  Carefully  screw  the  Starter  Charge  Cheaher  into  the  Ctebtr 
Ooi^ling  Collar  until  the  "neck"  is  proparly  seated  in  the 
L.l*.  Cowprassion  Chaster. 


4.- 


i 

i 


\ 

\ 


1.  auon  all  11am,  valves,  and  fittings  arv  proparly  tightened. 

2.  Clone  all  valves  on  test  fixture  sad  r— nve  air  compressor  line. 

3.  Opm  regulator  valve  on  th  gas  cylinder  and  set  gas  pressure  as 
per  test  requirements. 

4.  Open  bland  valve  staa  to  L.P.  drain  line. 

5.  mi  I ..P.  reservoir. 

6.  Oonaact  L.P.  fill  hose  to  L.P.  fill  valve. 

7.  Open  L.P.  reservoir  stopcock. 

8.  Open  L.P.  isolation  valve. 

9.  Manually  hold  poppet  valve  open. 

10.  Slowly  open  L.P.  fill  valve  until  a slow  flow  is  detected.  Match 
L.P.  reservoir  and  refill  when  necessary.  Mien  L.P.  is  seen  coning 
out  the  drain  line  and  no  air  bubbles  merge,  close  L.P.  fill  valve. 

11.  Nova  poppet  up  and  down  to  dislodge  any  air  bubbles. 

12.  dose  bleed  valve. 

13.  Close  stopcock  an  L.P.  reservoir. 

14.  Mmove  L.P.  fill  hose  free  L.P.  fill  valve  and  insert  hose  plug. 


For  ullage  addition  (if  no  ullage  go  to  atop  21. 


15.  Opeu  L.P.  fill  valve. 

16.  —wove  septus  access  plug. 

17.  Load  syringe  with  required 

18.  insert  syringe  thru  septan 
20.  Close  L.P.  fill  valve. 


j 

i 
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nt  of  gas. 
inject  gas  into 
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8 8 8 8 8 

i i si  s 

□ □ □ □ □ 

□ □ □ □ □ 
□ □ □ □ □ 


w-mnci 

1.  After  fining,  aniur*  all  nlvtf  cad  vast  plug*  «•  clo— d. 

2.  Clo—  L.V.  re— roolr  ttogOBCfc  dleoonnact  fill  Inn,  and 

3.  Dog  any  leftover  L.P.  froo  drain  bottle  into  and  L.P. 
toottlo  and  rano—  bottla. 


□ □ □ □ □ «• 


~ □ »• 


Main  sura  that  aafaty  interlock  bays  an  namd  fWa  tha 
ranota  firing  panal  and  tha  aolanold  val—  analog  —Itch. 


Z □ □ 


6.  Naha  certain  that  the  X-docir  Uel  la  out  of  tha  paid*  of 

tha  rant  gaa  and  do—  aafaty  ancloaura.  Bnaoaa  handle  and 
lnaart  latch  bolt. 


u 


n 7.  switch  eharga 


to  "operate* . 


Naha  certain  that  all  hacocd/Nold  battens  an  Bloantlan  are 
oat  In  the  "hi coed"  Position. 


z c 

i z 


14. 


switch  slanatlon  to  'hody" . 

Torn  on  lab  Naming  Ligit. 

In—  rt  safety  hay  into  firing  papal  and  actnata  aolarold 
firing  circuit. 

Go  to  i — nla  control  firing  panel,  in— rt  safety  interlock 
hay  and  tarn  to  "0e". 

Switch  aiming  to  "on". 

Pin. 


